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EXECUTIVE  SUMMARY 


BACKGROUND 

Carbon  fibers  bonded  in  epoxy  constitute  one  of  today's  "miracle" 
structural  materials.  For  particular  applications  structural  components  can 
be  designed  with  great  strength,  while  being  considerably  lighter  than  the 
conventional  parts  they  replace.  Carbon  fiber  composite  material  has  been 
used  to  date  in  sporting  goods  and  military  aircraft  primarily.  Limited  near- 
term  use  in  civil  aircraft  is  expected  to  grow  considerably  in  the  future. 

The  material  is  also  expected  to  be  used  in  automobiles.  In  view  of  the  many 
potential  applications  of  carbon-fiber  composites  considerable  concern  was 
engendered  as  the  result  of  evidence  that  individual  fiber  segments  could 
cause  electrical  and  electronic  equipment  to  fail  under  certain  operating 
conditions.  Such  individual  fiber  segments  could  be  released,  for  example, 
in  a  fire  involving  the  composite  material.  As  a  result  of  this  concern  a 
national  multi-agency  program  was  established  under  the  aegis  of  the  Office 
of  Science  and  Technology  Policy.  In  this  program  the  National  Aeronautics 
and  Space  Administration  was  assigned  the  responsibility  of  examing  the  risk 
due  to  possible  accidents  and  fire  involving  civil  aircraft,  with  carbon  fiber 
composite  structural  components. 

The  scenario  envisaged  in  the  NASA-funded  ORI  investigation  is: 

t  A  commercial  jet  aircraft  in  an  accident  with  fire  leading 

to  ... 
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fe  Release  of  substantial  numbers  of  individual  carbon  fiber  seg¬ 
ments  which  . . . 

e  Are  carried  far  from  the  accident  scene  by  prevailing  winds 

and  . . . 

•  Enter  buildings  and  parked  aircraft  in  their  path  causing  ... 

•  Failures  of  electrical  and  electronic  equipment  resulting  in  .. 

•  Economic  impact 

In  order  to  examine  the  potential  magnitude  of  the  omerall  effect  ORI,  Inc 
developed  a  computer  simulation  model  which  replicates  these  events  as  far  as 
possible.  In  its  final  Phase  I  Report,  published  in  May  1979  ORI  described 
its  initial  modeling  effort  and  presented  the  results  obtained  from  a  large 
number  of  simulation  runs  for  nine  major  airports  in  the  United  States,  using 
the  best  estimates  available  for  all  required  input  data.  The  airport  results 
were  combined  statistically  to  obtain  an  estimate  of  the  total  national  risk. 

In  Phase  II  the  ORI  airport  risk  assessment  model  was  extended  in 
several  respects,  principally  to  increase  the  variance  --  i.e.  improve  the 
likelihood  that  extreme  values  would  be  generated.  Additional  experimental 
data  made  available  by  other  NASA-sponsored  efforts  were  used  in  the  new 
calculations.  A  new  national  risk  model  was  developed;  it  was  designed  to  be 
more  useful  for  estimating  the  statistical  confidence  that  can  be  assigned 
to  the  results.  An  essentially  independent  study  was  mounted  to  assess  the 
risk  to  the  electric  power  distribution  system. 

SINGLE  AIRPORT  RISK  ASSESSMENT 

Method 

The  method  used  to  estimate  the  risk  associated  with  accidents  in¬ 
volving  aircraft  with  carbon  fiber  composite  structural  components  is  essentially 
a  Monte  Carlo  simulation  model.  The  method  requires  the  generation  of  many 
aircraft  accidents,  with  variables  used  in  the  calculation  drawn  at  random 
from  defined  probability  instructions.  The  impact  of  each  accident  is  cal¬ 
culated  and  saved;  after  many  accidents  have  been  simulated  the  model  computes 
several  statistical  measures  from  the  results  of  all  the  accidents.  The 
method  is  illustrated  schematically  in  Figure  1, 


FIGURE  1.  SCHEMATIC  REPRESENTATION  OF  ORI  AIRPORT  CF  RISK  ASSESSMENT 
MODEL.  SHADED  BOXES  REPRESENT  SINGLE  ACCIDENT  SCENARIO. 
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The  principal  elements  or  submodels  in  the  improved  ORI  Phase  II 
Risk  Assessment  Model  are: 

•  Random  accident  generation.  For  each  replication,  comprising 

a  year's  accidents,  the  model  first  determines,  using  appropriate 
stochastic  methods,  the  actual  number  of  accidents  for  each 
size  aircraft.  The  mean  annual  fire-accident  rate  is  based  on 
the  analysis  of  historical  accident  data  conducted  by  the  major 
airframers.  This  submodel  makes  use  of  aircraft  manufacturers' 
projections  of  increased  use  of  carbon  fiber  composite  material 
in  individual  aircraft,  and  the  changes  in  the  fleet  mix  between 
now  and  1993,  the  target  year  for  the  risk  assessment.  The  air¬ 
framers  also  prepared  a  detailed  analysis  of  historical  jet  air¬ 
craft  accidents  involving  fires,  which  provided  critical  inputs 
describing  accident  impact.  This  submodel  randomly  selects  the 
aircraft  involved  in  each  simulated  accident,  based  on  the  pro¬ 
jected  fleet  mix  -  and  randomly  determines  the  extent  of  damage  to 
the  aircraft.  The  location  is  also  determined  by  drawing  a  random 
sample  from  a  distribution  obtained  from  the  historical  data. 


•  For  the  specific  accident  characteristic  the  resulting  fire 
plume  is  modelled  using  standard  methods.  This  determines  the 
height  to  which  the  fibers  are  carried. 

•  An  improved  Phase  II  transport  and  diffusion  calculation  deter¬ 
mines  the  concentration  of  individual  fibers  at  selected  repre¬ 
sentative  points  downwind  from  the  accident  scene.  Weather  vari¬ 
ables  used  in  the  calculation  are  selected  at  random  from 
historical  data  for  each  airport. 

•  County-based  data  are  used  to  describe  the  numbers  of  each  type 

of  vulnerable  business  and  industrial  facility,  and  housing  patterns 
within  a  50-mile  range  of  the  airport  in  all  directions. 

•  Each  household,  and  each  size  and  type  of  business  facility  is 
characterized  by  a  specific  type  of  building  or  other  structure. 

A  transfer  submodel  determines  the  fraction  of  fibers  outside 
the  building,  determined  in  the  transport  and  diffusion  calcula¬ 
tion,  that  gets  inside  the  building. 


•  For  each  type  and  size  of  business  at  risk  a  production  model 
describes  the  power  flow  for  that  type  of  business.  This  pro¬ 
duction  model  is  used  with  the  interior  dosage  and  equipment 
failure  inputs,  based  on  current  experimental  data,  to  deter¬ 
mine  which  equipments  fail,  and  the  resulting  impact  on  each 
business.  A  similar  calculation  is  made  for  each  class  of 
vulnerable  aircraft  at  the  airport  itself. 

•  Repair  costs  are  computed  for  each  piece  of  equipment  that  fails 
in  the  parked  aircraft,  and  in  each  business  and  industrial 
facility.  In  addition,  we  determine  on  a  stochastic  basis  those 
business  establishments  that  close  as  a  result  of  the  cumulative 
effects  of  individual  equipment  failures.  The  economic  impact 
of  such  closings  is  determined  by  allocating  to  each  business 
establishment  its  share,  based  on  payroll,  of  the  Gross  Domestic 
Product  for  the  particular  type  of  business  represented  by  that 
establishment. 

The  methods  briefly  described  above  are  applied  repetitively  to  a 
large  number  of  randomly  generated  accidents  at  one  airport.  The  result  is 
the  developm.ent  of  a  set  of  accident  impact  costs  for  many  replications  of  the 
year  1993  at  each  airport.  These  results  are  then  examined  to  provide  statistical 
measures  of  the  risk.  Sample  results  are  described  below. 

Results 

The  average  annual  impact  (in  1976  dollars)  and  the  average  impact 
per  accident  for  each  of  the  airports  analyzed  in  the  standard  1993  scenario, 
are  summarized  in  Table  1.  The  average  impact  over  all  simulated  accidents 
(about  2250)  at  all  airports  is  $5  for  household  equipment  damage,  $172  for 
business  and  industry  equipment  repair  and  business  dislocation,  and  less  than 
one  dollar  for  repair  of  damaged  avionics  equipment.  In  addition  to  the  mean 
values,  the  analysis  of  the  simulation  results  provides  considerable  statistical 
insight  into  the  results.  For  example,  the  most  costly  accidents  generated 
at  any  of  the  airports,  for  each  of  the  major  impact  categories  were: 

$2,665  in  household  damage:  Kennedy  Airport,  New  York 

$274,000  in  business/industrial  impact:  Logan  Airport,  Boston 
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TABLE  1 


SUMMARY  OF  INDIVIDUAL  AIRPORT  RISK  ASSESSMENT  RESULTS  -  1993 


Airport 

Mean  Impact  per 
Accident  ($) 

Mean  Annual 
Impact  ($) 

O' Hare/Chicago 

169 

17 

John  F.  Kennedy/New  York  City 

212 

15 

Washington  National  Airport/Washington,  D.C. 

315 

12 

Lambert/St.  Louis 

69 

3 

LaGuardia/New  York  City 

384 

24 

Logan/Boston 

153 

9 

Harts field/Atlanta 

73 

8 

Miami  International/Miami 

31 

2 

Philadelphia  International /Philadelphia. 

200 

8 

VI 


$3,910  in  avionics  equipment  damage:  Kennedy  Airport,  New  York. 

These  results  are  from  a  total  of  approximately  2500  simulated  accidents  at 
each  of  the  airports,  or  from  a  total  of  more  than  22,500  accidents  simulated 
at  all  airports. 

A  risk  profile  prepared  from  the  computer-generated  results  for 
Washington  National  Airport  appears  in  Figure  2.  It  shows  that,  although 
the  average  annual  cost  resulting  from  these  accidents  is  $12,  there  is 
some  chance  of  exceeding  this  figure  by  a  considerable  amount.  However,  the 
probability  that  the  annual  impact  will  exceed  $100,000  is  only  .00003  (3  in 
100,000).  Expressed  another  way  this  implies  that  costs  of  this  magnitude 
might  be  incurred  once  every  33,333  years,  on  the  average.*  The  computed 
statistical  confidence  limits  applied  to  these  results  indicate  that  we  can  be 
quite  confident  that  the  statistical  uncertainties  inherent  in  the  computer 
simulation  method  would  not  cause  us  to  raise  these  probabilities  significantly. 

In  addition  to  the  purely  statistical  uncertainty  involved,  the 
model  was  also  used  to  examine  the  sensitivity  of  the  results  to  possible  input 
data  errors.  This  can  be  done  relatively  easily  by  changing  particular  input 
data  elements  and  rerunning  the  model.  Several  such  sensitivity  tests  were 
conducted;  none  indicate  that  the  impact  of  input  errors  or  changes  in  assump¬ 
tions  would  require  us  to  significantly  change  our  conclusions  regarding  the 
nature  of  the  risk.  An  example  of  one  rather  drastic  sensitivity  calculation 
is  the  one  conducted  for  O'Hare  Airport.  In  this  case  we  changed  the  inputs 
to  reflect  an  assumption  that  all  aircraft  with  composite  operating  at  O' Hare 
were  loaded  —  that  is;  they  were  all  heavy  jets  with  15,619  kilograms  of  com¬ 
posite  material  onboard.  In  our  standard  1993  fleet  this  aircraft  only  com¬ 
prises  about  a  half  of  one  percent  of  the  aircraft  with  carbon  fiber  composite 
components.  The  average  amount  of  composite  on  the  'standard'  1993  aircraft 
with  carbon  fiber  composite  aboard  is  about  2800  kilograms.  This  worst-case 
O' Hare  Airport  risk  profile  is  compared  with  the  standard  case  in  Figure  3. 

The  result  shows  that,  even  in  this  "worst  case",  the  probability  of  exceeding 
an  annual  impact  of  $100,000  is  approximately  .003  (3  in  1000)  at  the  nation's 

*Since  the  likelihood  of  more  than  one  accident  at  the  airport  in  one  year  is 
very  smal 1 ,  this  may  be  safely  paraphrased  as:  "one  such  accident  might  occur 
every  33,333  years,  on  the  average." 


Probability  of  Exceeding 


100  1,000  10,000  100,00C 

Annual  Cost  in  Dollars 


FIGURE  2.  1993  WASHINGTON  NATIONAL  AIRPORT  ANNUAL  RISK  PROFILE 


viii 


1/10 


FIGURE  3 


c 

o 

s: 

(/) 

■M 

O 

CD 

Q. 


O) 

C 

'"O 

0 

0 

O 

X 

UJ 

H- 

o 


-Q 

0 

JD 

O 


1/10,000* - j ^ - 1 

100  1,000  10,000  100,000 


Annual  Impact  in  Dollars 


COMPARISON  OF  ANNUAL  1993  O' HARE  AIRPORT  RISK  PROFILES: 
STANDARD  AIRCRAFT  MIX  AND  ALL  CF  AIRCRAFT  WITH  MAXIMUM 
AMOUNT  OF  COMPOSITE 
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busiest  airport.  The  results  also  show  that  the  mean  annual  impact  increased 
by  a  factor  of  ten  from  $17  for  the  standard  aircraft  mix  to  $172  for  all  air¬ 
craft  "loaded"  with  CF. 

NATIONAL  RISK 

In  order  to  estimate  the  total  national  risk  due  to  release  of  CF 
in  an  aircraft  fire-accident,  and  the  subsequent  damage  to  electric  and 
electronic  equipment,  we  assumed  that  the  nine  airports  encompassed  all  of  the 
commercial  aircraft  activity  in  the  United  States.  This  greatly  overestimates 
the  risk  since  these  tend  to  be  busy  airports  with  considerable  surrounding  busi- 
business  and  industry.  The  model  generates  a  random  number  of  accidents  with 
fire  occuring  in  the  entire  United  States,  based  on  the  mean  values  determined 
in  the  previously-referred  to  airframers'  analysis  of  historical  accidents. 

The  individual  accidents  are  assigned  to  one  of  the  airports  previously  analyzed 
according  to  the  relative  traffic  level  at  that  airport.  The  impact  of  that 
accident  is  obtained  by  randomly  drawing  one  of  the  accidents  that  was  pre¬ 
viously  simulated  at  that  airport.  The  result  is  a  conservative  estimate  (that 
is  —  on  the  high  side)  of  the  national  risk,  in  that  accidents  will  not  be 
allocated  to  other  low-risk  airports.  The  resulting  national  annual  risk  pro¬ 
file  is  shown  in  Figure  4.  A  summary  of  the  results  also  appears  in  Table  2, 
where  we  have  separated  out  the  avionics  equipment  impact  because  of  special 
interest  in  that  aspect  of  the  risk  assessment.  The  risk  profile  may  be  inter¬ 
preted  as  showing  that  the  probability  of  exceeding  $100,000  in  annual  impact 
is  approximately  .00015;  the  estimated  probability  of  exceeding  an  economic 
impact  of  $1,000,000  is  less  that  than  .00001.  (one  year  out  of  a  hundred 
thousand)  The  tests  of  statistical  confidence  and  the  sensitivity  tests  con¬ 
ducted  during  the  study  indicate  that  we  may  be  confident  that  these  results 
are  statistically  valid  and  conservative,  in  that  they  tend  to  overestimate 
the  actual  risk. 

IMPACT  ON  ELECTRIC  UTILITY  SYSTEMS 

Parallel  efforts,  primarily  under  the  sponsorship  of  the  Department 
of  Energy,  investigated  the  vulnerability  of  electrical  transmission  equipment 
to  carbon  fiber  incursion.  These  indicated  some  vulnerability  of  individual 
components  in  the  electrical  distribution  system,  but  tend  to  show  that  high 
high  voltage  (above  38  kilovolts)  and  low  voltage  systems  (below  2.4  kilovolts) 
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Probability  of  Exceeding  Impact  Shown 


TABLE  2 


NATIONAL  RISK  IMPACT  RESULTS 


Type  of  Cost 

Mean 

Worst  Case 

Business,  Industry,  and 
Household  Impact 

$466 

$275,000 

Avionics  Equipment 

Failures 

$2 

$3,900 

Annual  Economic  Impact  (Dollars) 


FIGURE  4.  1993  ANNUAL  NATIONAL  RISK  PROFILE 
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are  both  essentially  invulnerable.  In  the  risk  assessment  model  described 
above,  the  vulnerability  of  equipment  in  the  intermediate  range  was  accounted 
for  in  specific  cases  by  examining  the  probability  of  failure  of  stepdown  trans¬ 
formers  for  specific  industrial  facilities.  At  NASA's  request  ORI  conducted  a 
separate  analysis  to  estimate  the  impact  of  possible  aircraft-fire  accident 
on  the  total  national  electric  power  distribution  system;  this  was  conducted 
in  parallel  with  a  historical  review  of  electrical  power  outage  data. 

In  order  to  conduct  the  analysis,  typical  and  actual  distribution 
systems  operating  in  the  vulnerable  voltage  range  were  defined.  Estimates 
of  the  effects  of  carbon  fibers  on  the  circuits  were  based  on  experimental 
data  made  available  to  ORI  by  the  Department  of  Energy  team  under  NASA  auspices. 
In  order  to  be  conservative  this  calculation  assumed  a  downwind  pattern  of 
carbon  fiber  exposure  based  on  a  combination  of  parameters  defined  as  the 
worst  possible  case.  Individual  insulator  and  bushing  failure  probabilities 
are  based  on  these  worst  case  exposure  values  and  experimental  failure  data 
for  these  components. 

The  results  were  obtained  for  several  sets  of  conditions.  The  cir¬ 
cuits  examined  include  a  typical  (textbook)  electrical  distribution  circuit 
operating  in  the  7.5  KV  range,  and  an  actual  suburban  system  provided  by  an 
operating  electric  utility  system  in  the  23  KV  range.  For  the  typical  cir¬ 
cuit,  failures  were  computed  based  on  published  component  reliability  data,  and 
the  worst-case  carbon  fiber  scenario.  For  the  23  KV  circuit  we  obtained  actual 
reported  outages,  and  also  computed  the  worst-case  carbon  fiber  scenario 
failures.  In  all  cases  the  failures  were  extrapolated  to  a  national  base, 
assuming  that  all  users  were  served  by  the  system  under  examination,  and  using 
the  national  average  annual  number  of  carbon  fiber  aircraft-fire  accidents  for 
1993.  The  results  are  summarized  in  Table  3.  As  a  further  comparision  we 
note  that  annual  bulk  outages  reported  to  the  Federal  Energy  Regulatory  Com- 
mision  indicate  that,  typically,  3,000,000  (3  x  10^)  customers  suffer  an 
outage  in  one  year  for  an  annual  average  of  about  0.05  outages  per  utility 
customer.  Such  bulk  outages  comprise  an  interruption  occuring  at  69KV  and 
above  and  resulting  in  a  loss  of  at  least  100  megawatts  for  at  least  15  minutes, 
or  a  loss  of  more  than  one  half  of  a  small  system's  annual  peak  load.  These 
outages  clearly  comprise  only  a  fraction  of  the  total  outages  in  the  nation. 


TABLE  3 

COMPARISON  OF  ANNUAL  ACTUAL  AND  "NORMAL"  ELECTRIC  POWER  DISTRIBUTION 
SYSTEM  OUTAGES  WITH  WORST-CASE  CARBON  FIBER  RELATED  OUTAGES 


Measure 

Distribution  Circuit 

Typical : 

7.5  Kv 

Actual :  23  Kv 

"Normal" 

Worst-Case  CF 

Reported 

Worst-Case  CF 

Total  National 
Outages  per  Year 

22  X  10® 

23 

140  X  10® 

3300 

Annual  Outages 
per  Customer 

0.32 

10-« 

2.1 

5  X  10"® 

All  of  the  carbon-fiber  related  outage  results  are  based  on  assump¬ 
tions  that  are  all  on  the  conservative  side;  that  is,  they  overestimate  the 
resulting  number  of  expected  power  outages.  The  result  is  that  we  can  expect 
less  than  one  civil  aircraft  fire-accident  carbon  fiber  release  related  power 
outage  for  about  every  200,000  to  million  outages  that  occur  for  a  variety  of 
other  reasons. 

Clearly,  then,  we  conclude  that  the  carbon  fibers  potentially  released 
in  a  civil  aircraft  accident  with  fire  represent  a  relatively  insignificant 
threat  to  the  electric  power  distribution  network. 
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I.  INTRODUCTION 


This  is  the  Final  ORI  Report  on  Phase  II  of  its  Civil  Aviation  Car¬ 
bon  Fiber  Risk  Assessment  study  performed  for  the  National  Aeronautics  and 
Space  Agency  under  Contract  No.  NASI -15379.  The  NASA-funded  effort  is  part 
of  a  major  national  program  directed  toward  estimating  the  potential  risk 
of  increased  use  of  carbon  fiber  composite  material  in  a  variety  of  applica¬ 
tions.  This  program  was  initiated  as  the  result  of  evidence  that  electrical 
and  electronic  equipment  may  fail  as  the  result  of  the  deposit  of  carbon  fibers 
released  by  burning  of  the  composite  material.  Carbon  fiber  -  expoy  composite 
materials  offer  considerable  advantage  over  more  conventional  material  due  to 
the  ability  to  engineer  in  superior  strength  while  achieving  a  considerable 
weight  saving.  The  national  program  was  established  to  investigate  the  nature 
of  this  potential  hazard  in  the  light  of  projected  increased  use  of  these 
materials.  The  NASA  program,  of  which  the  ORI  investigation  reported  here  is 
a  small  part,  is  directed  particularly  at  the  possible  risk  associated  with  the 
use  of  carbon  fiber  composite  materials  in  civil  aviation. 

In  Phase  I  of  its  investigation  ORI  developed  a  computer  simulation 
model  that  was  used  to  generate  risk  statistics  for  accidents  at  several  air¬ 
ports,  which  were  later  combined  to  estimate  the  national  risk.  The  Phase  I 
model,  although  using  many  Monte  Carlo  -  or  stochastic  -  submodels,  did  compute 
the  business-industry  impact  on  an  expected  value  basis.  This  may  have  tended 
to  limit  the  variance  of  the  final  results  and  thus  reduced  the  likelihood  of 
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generating  extreme  values  on  the  "tails"  of  the  accident  cost  distribution. 

In  the  Phase  II  model  this  computation  has  been  made  stochastic,  and  several 
other  subroutines  in  the  complete  model  have  similarly  been  made  to  operate  in 
a  random  rather  than  an  expected  value  mode. 

The  calculation  is  essentially  "input  driven";  that  is,  the  results 
depend  on  a  host  of  input  data,  many  of  which  are  from  sources  not  directly 
linked  to  the  problem  at  hand  (for  example,  national  economic  data).  Other  input 
data  elements,  such  as  the  amount  of  carbon  fiber  composite  on  an  airplane,  and 
the  fraction  of  carbon  that  would  be  released  if  that  airplane  were  to  crash 
and  burn,  are  documented  for  the  Phase  I  calculation.  For  Phase  II  additional 
experimental  data  were  available,  and  were  used.  Thus  the  results  presented 
in  this  report  are  based  on  the  use  of  a  computer  model  that  is  more  sophis¬ 
ticated  than  the  one  used  in  Phase  I,  operating  on  a  much  more  solid  data  base. 

In  addition  ORI  was  requested  to  investigate  the  potential  risk  to  the  electric 
utility  subdistribution  system. 

The  basic  technical  approach  to  the  risk  assessment  problem  is  to 
simulate  many  aircraft  accidents  with  fire,  each  one  characterized  by  many 
random  variables,  and  then  compile  statistics  based  on  the  analysis  of  the 
computed  impact  of  the  series  of  accidents.  The  availability  of  high-speed 
digital  computing  techniques  makes  this  approach  feasible.  This  is  an  applica¬ 
tion  of  the  so-called  Monte  Carlo  simulation  technique.  The  principal  elements 
in  the  scenario  that  is  simulated  are: 

•  Aircraft  accident  with  fire 

•  Release  of  carbon  fiber  material 

•  Entrainment  of  the  carbon  fibers  in  a  smoke  plume 

e  Transport  of  the  carbon  fiber  material  downwind 

•  Transfer  of  some  of  the  fibers  into  the  interior  of  buildings 

§  Failures  of  electrical  and  electronic  equipment 

•  Economic  impact  of  these  failures 

In  addition  to  the  simulation  of  these  events  £er  ^  the  complete  model  must 
perform  many  other  functions  related  to  the  selection  of  appropriate  random 
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variables,  as  well  as  what  may  be  termed  "housekeeping"  functions.  A  simpli¬ 
fied  form  of  the  complete  airport  risk  assessment  model  is  illustrated  in  the 
flowchart  appearing  in  Figures  1.1a  and  1.1b. 

The  simulation  of  one  accident  requires  the  random  selection  of  the 
accident  location.  This  selection  depends  on  input  data  that  is  the  result  of 
a  detailed  analysis  of  all  historical  jet  aircraft  accidents  in  which  fires 
were  involved,  performed  under  NASA  sponsorship  by  the  principal  airframe 
manufacturers.  The  estimate  of  the  amount  of  fibermaterial  released  depends 
on  inputs  that  define  the  mix  of  aircraft  in  the  fleet  for  the  target  year  of 
1993.  The  fraction  of  the  material  that  is  released  as  fibers  in  the  size 
range  of  interest  is  based  on  recent  experimental  results.  The  computer  pro¬ 
gram  models  the  behavior  of  the  resulting  fire  plume  that  carries  the  released 
fibers  aloft.  The  downwind  transport  and  diffusion  processes  are  then  modelled 
using  methods  that  are  somewhat  more  general  and  refined  compared  to  those 
used  in  Phase  I.  The  necessary  meteorological  inputs  for  this  calculation 
are  drawn  at  random  from  a  body  of  data  for  each  of  the  airports  for  which 
the  calculations  are  m.ade. 

The  transport  and  diffusion  calculation  provides  the  exposure  or 
dosage  at  predefined  points.  These  points  are  selected  in  advance  to  be  repre¬ 
sentative  of  the  area  at  risk  surrounding  the  airport.  In  all  cases  these 
extend  to  a  range  of  at  least  50  miles  from  the  airport  at  which  the  simulated 
accidents  occur.  The  points  are  selected  to  represent  concentrations  of 
businesses,  industry,  and  private  residences  in  individual  counties,  in  order 
to  make  use  of  readily  available  county-based  economic  and  other  census  data. 
All  types  of  vulnerable  businesses  and  industry,  as  well  as  households,  are 
characterized  by  particular  types  of  buildings  at  the  representative  points. 

The  definition  of  the  building  types  includes  a  set  of  ventilation  parameters, 
obtained  from  standard  engineering  sources  modified  by  particular  experimental 
data  appropriate  to  the  carbon  fibers.  These  parameters  are  used  in  the  cal¬ 
culating  the  fraction  of  the  fibers  outside  each  building  that  would  get  in¬ 
side,  termed  the  transfer  function.  In  this  way  the  risk  assessment  model 
determines  the  exposure  or  dosage  to  which  vulnerable  equipment  is  subjected. 
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SET  AIRPORT 


FIGURE  1.1  a.  ORI  RISK  ASSESSMENT  MODEL  FLOWCHART 

(continued) 
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Si 


FIGURE  1.1  b.  ORI  RISK  ASSESSMENT  MODEL  FLOWCHART 

(concluded) 


With  the  interior  exposure  available  we  can  then  compute  the  pro¬ 
ability  of  failure  for  specific  classes  of  vulnerable  equipment,  previously 
defined  for  each  type  of  business,  industry,  and  household  in  the  airport 
environs.  The  failures  of  individual  equipments  are  then  determined  using 
randomized  or  Monte  Carlo  methods.  The  impact  of  these  failures  is  then 
assessed  in  several  ways.  First  all  equipment  failures  are  totalled  to  esti- 
m.ate  the  required  cost  to  repair  the  damaged  equipment.  The  likelihood  that 
the  equipment  failure  would  be  severe  enough  to  cause  a  place  of  business  to 
close  is  then  examined.  This  is  based  on  the  assumption  that  the  business 
facility  would  close  if  electric  power  were  lost,  if  its  principal  control 
systems  were  knocked  out,  or  if  half  of  its  production  equipment  fails.  These 
events  are  also  determined  on  a  random  basis  for  essentially  each  place  of 
business  in  the  downwind  path  of  the  plume  resulting  from  the  and  accident 
and  fire.  The  financial  impact  of  such  closings  is  estimated  by  allocating 
to  each  place  of  business  its  estimated  share  of  the  Gross  Domestic  Product 
for  that  class  of  business  or  industry.  If  a  place  of  business  is  determined 
to  have  closed  as  a  result  of  the  carbon  fiber  release  incident,  a  clean  up 
cost  is  also  assessed.  The  household  equipment  impact,  because  of  the  large 
number  of  essentially  identical  equipments  at  risk  in  very  similar  environ¬ 
ments  is  treated  on  an  expected- val ue  basis. 

The  vulnerability  of  avionics  equipment  aboard  aircraft  parked  at 
the  airport  is  also  examined  in  the  risk  assessment  calculation.  The  number 
of  aircraft  in  a  potentially  vulnerable  state  at  the  airport  is  determined 
from  data  provided  by  the  aircraft  manufacturers  via  NASA.  For  each  air¬ 
craft  the  calculation  determines  the  number  of  each  class  of  equipment  in 
each  of  several  at-risk  states.  The  number  of  failures  is  then  determined  on 
a  random  basis,  and  the  input  repair  costs  are  used  to  determine  the  total 
impact  of  such  failures. 

After  all  business  facilities  and  households  at  risk,  and  all  parked 
aircraft  have  been  examined,  the  m.odel  has  generated  an  estimate  of  the  total 
impact  of  one  accident,  in  1976  dollars.  These  results  are  available  for  the 
three  principal  impact  categories:  business/industry,  household,  and  avionics 
equipment  in  parked  aircraft.  The  computer  program  management  module  then 
returns  to  the  "front  end"  to  generate  the  details  associated  with  the  next 
accident  in  the  sample  year  being  replicated,  if  there  are  any  more.  Once  all 
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accidents  in  a  sample  year  have  been  simulated  the  model  generates  another 
sample  year,  with  its  aircraft  accidents,  simulates  the  details  of  each  acci¬ 
dent,  then  does  the  next  sample.  This  process  continues  until  the  preset 
number  of  annual  replications  has  been  completed.  At  this  point  we  have  the 
impact  for  many  sample  years  and  compute  the  statistics  over  all  samples. 

These  include  such  measures  as  the  mean  annual  impact,  and  the  mean  impact 
per  accident,  as  well  as  the  risk  profiles.  The  risk  profile  shows  graphically 
the  probability  that  the  annual  impact  will  exceed  any  value. 

The  method  described  briefly  above  has  been  applied  to  a  sample  of 
nine  airports.  These  were  selected  as  reasonably  representative  of  U.S.  air¬ 
ports  in  Phase  I,  although  they  were  purposely  chosen  to  be  a  conservative 
group,  in  the  sense  of  representing  the  "high  side"  of  the  risk.  In  order  to 
compute  the  total  risk  at  the  national  level,  this  set  of  airports  was  assumed 
to  represent  the  entire  United  States.  The  national  risk  model  generates  ran¬ 
dom  accidents  during  a  sample  year  and  determines,  on  a  random  basis,  which 
of  the  nine  airports  that  each  accident  would  have  taken  place.  A  pre¬ 
viously  simulated  accident  at  that  airport  is  selected  at  random,  and  its  com¬ 
puted  impact  added  to  the  running  total  for  the  current  replication's  national 
impact.  By  repeating  this  process  many  times  the  national  model  generates  the 
statistics  necessary  to  produce  the  national  risk  profile. 

The  methods  outlined  here  are  described  in  detail  in  the  remainder 
of  this  report.  The  accident  details,  including  extrapolation  of  the  1993  com¬ 
mercial  aircraft  fleet  mix,  and  other  necessary  inputs  are  described  in  Section 
II  which  follows  immediately.  The  fire  plume  calculation  is  described  in 
detail  in  Section  III. 

The  improved  Phase  II  methods  for  computing  the  downwind  transport 
and  diffusion  of  the  material  contained  in  the  plume  are  described  in  Section 
IV.  The  methods  used  to  compute  the  transfer  of  the  diffused  material  into 
the  interior  of  buildings  and  other  structures,  including  the  use  of  new  Phase 
II  data,  appears  in  Section  V.  The  following  part  of  the  report.  Section  VI, 
discusses  the  methods  of  treating  equipment  failures.  At  this  point  in  the 
logic  flow  of  the  simulation  we  are  ready  to  calculate  failures  and  need  to 
convert  those  failures  into  dollar  measures  of  impact.  The  required  methods 
are  described  in  detail  in  Section  VII;  this  part  of  the  methodology  comprises 
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the  major  Phase  II  improvement  over  Phase  I.  The  required  data  bases  for  the 
economic  inputs  are  also  detailed  in  that  section.  Results  of  the  single  air¬ 
port  simulations  are  presented  in  Section  VIII,  including  several  sensitivity 
tests,  which  examine  the  impact  on  the  risk  results  of  significant  changes  in 
input  data  and  associated  assumptions.  The  national  model  and  the  results  it 
generated  are  presented  in  Section  IX.  The  ORI  conclusions  appear  in  Section 
X.  The  analysis  of  failures  in  the  electrical  power  distribution  system  is 
described  in  Appendix  A. 
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II.  ACCIDENT/RELEASE 


In  this  section  of  the  report  we  describe  the  method  used  to  "gen¬ 
erate"  an  accident  and  determine  the  amount  of  fiber  released  in  the  accident. 

The  model  is  applied  to  a  single  airport  at  a  time.  The  generation  of  an  acci¬ 
dent  with  fire  involving  an  aircraft  with  CF  composite  in  its  structure  is  the 
first  step  in  the  scenario  simulated  by  the  ORI  risk  assessment  model. 

1993  FLEET  COMPOSITION 

In  order  to  estimate  the  amount  of  carbon  fiber  that  might  be  re¬ 
leased  in  an  accident  with  fire  it  was  first  necessary  to  estimate  the  amount 
of  carbon  fiber  that  would  be  on  particular  aircraft,  as  well  as  the  mix  of 
aircraft  in  the  1993  fleet.  The  principal  aircraft  manufacturers,  working  with 
NASA,  and  in  consultation  with  ORI,  prepared  descriptions  of  the  different  air¬ 
craft  configurations  to  be  introduced  from  now  until  1993.  These  are  defined  by 
the  amount  of  composite  material  in  all  structural  components;  each  of  the  air¬ 
craft  types  is  defined  in  Table  2.1.  Several  aircraft,  defined  early  in  Phase  II, 
were  later  dropped  when  it  was  determined  that  it  was  unlikely  that  they  would 
be  in  the  1993  fleet.  For  this  reason  no  aircraft  of  types  3  and  4  appear  in 
Table  2.1.  Several  different  aircraft  defined  by  the  airframe  manufacturers 
that  were  essentially  identical  from  the  composite  distribution  viewpoint  were 
combined  in  preparing  the  table.  Retirement  schedules  were  developed,  and 
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TABLE  2.1 

DEFINITION  OF  TYPES  OF  AIRCRAFT  EXPECTED  TO  BE  IN  THE  1993 
COMMERCIAL  FLEET  BY  SIZE  AND  AMOUNT  OF  COMPOSITE  (KILOGRAMS) 


17.  Floor  Beams  0  0  0  0  0  558  0  0  0  909  909  909  0  711  711  711  711 

18.  Cockpit  Window  Sup.  00000000000000000 

_ Composite  65  192  367  582  1505  3786  155  286  1234  2018  5125  6071  17R9  R4t;i  15  619 

Note;  Types  3  and  4  defined  early  in  Phase  II  were  later  dropped.  Types  20,  21,  22  not  shown  in  table  are  aircraft  with  no  CF. 


introduction  of  new  aircraft  "played"  for  each  year.  This  straightforward  cal¬ 
culation  led  to  the  development  of  an  estimate  of  the  1993  fleet  mix.  The  re¬ 
sult  is  shown  in  Table  2.2. 

ACCIDENT  RATE 

In  Phase  I,  ORI  conducted  a  limited  analysis  of  individual  aircraft 
accident  reports  and  summary  data  available  through  the  National  Transportation 
Safety  Board.  In  Phase  II,  under  NASA  auspices,  the  major  aircraft  manufac¬ 
turers  completed  detailed  analyses  of  approximately  100  jet  aircraft  accidents 
in  which  fire  played  a  part.  These  analyses  provided  estimates  of  the  damage 
to  each  major  aircraft  structural  component.  Based  on  this  data  base,  it  was 
determined  that  the  annual  fire-accident  rate  pertinent  to  the  risk  assessment 
was  3.8  per  year;  this  has  been  accepted  as  the  best  estimate  available  for 
the  1993  scenario.  For  the  risk  assessment  calculation  we  are  only  concerned, 
with  aircraft  containing  composite  material,  estimated  to  be  about  70  percent 
(cf.  Tables  2.1,  2.2)  of  the  1993  fleet,  for  a  resulting  national  mean  number 
of  carbon-fiber  aircraft  accidents  with  fire  of  2.6  per  year. 

The  simulation  model  treats  one  aircraft  size  at  a  time.  Accordingly, 
for  airport  A  and  aircraft  of  size  S,  we  estimate  the  annual  accident-with- 
fire  rate  by: 

X(A,S)  =  '^A  S  ^ 

^  ^  N 

S 

A  S 

where  N^  ^  is  the  number  of  operations  of  aircraft  of  size  S  at  airport  A;  thus 
the  sum 


comprises  all  operations  in  the  U.S.  In  any  one  replication  (a  random  year) 
the  number  of  accidents  is  assumed  to  fit  a  Poisson  distribution.  The 
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TABLE  2.2 

1993  COMMERCIAL  AIRCRAFT  FLEET  PROJECTION 


Aircraft 

Size 

ORI  Identification 
Number 

Number  in 
Fleet 

Composite  per 
aircraft 
(kilograms) 

1 

71 

65 

Smal  1 

2 

80 

192 

20 

409 

0 

5 

125 

367 

6 

37 

582 

Medium 

7 

255 

8 

329 

3,786 

21 

34 

0 

9 

7 

155 

10 

5 

286 

11 

326 

1 ,234 

12 

54 

2,018 

13 

80 

5,125 

Large 

14 

192 

6,071 

15 

53 

1 ,282 

16 

53 

3,038 

17 

143 

4,633 

18 

79 

5,413 

19 

11 

15,619 

22 

U— _ 

396 

0 

2~ 


probability  of  exactly  n  accidents  with  fire  involving  aircraft  of  size  S  at 
airport  A  is  given  by  ^ ^ : 

P(n;  X(A,S))  =  e"  ^ 

n! 

The  Monte  Carlo  simulation  model  makes  a  random  draw  from  this  distribution 
during  each  replication  for  one  airport,  and  one  size  aircraft.  Because  the 
number  of  accidents  in  any  one  replication  is  very  small  the  computer  model 
actually  uses  a  double  precision  sampling  technique. 

AMOUNT  OF  CARBON  FIBER  RELEASED 

The  airframer  accident  analysis  generated  estimates  of  the  amount 
of  each  major  aircraft  component  that  was  involved  in  each  accident.  These 
output  results  were  combined  with  the  characterization  of  each  aircraft--amount 
of  composite  in  each  component--to  provide  an  estimate  of  the  amount  of  com¬ 
posite  that  would  have  been  involved  in  each  of  the  historical  fire  accidents. 

For  each  projected  1993  aircraft  type  identified  in  Table  2.1  the  ORI  risk 
assessment  team  computed  the  sum: 

Composite  Consumed  =  S  (Fraction  Consumed)  x  (Amount  of  Composite) 

c  c  c 

for  all  accidents  in  the  airframers'  analysis,  where  the  index  c  refers  to  an 
aircraft  component.  Thus,  for  one  aircraft  type,  defined  by  a  distribution 
of  composite  material,  we  estimated  the  total  amount  of  composite  material  that 
would  have  been  consumed  in  each  of  the  analyzed  historical  accidents.  The 
results,  a  sample  of  which  are  shown  in  Table  2.3,  comprise  one  of  the  major  in¬ 
put  data  sets  for  the  risk  calculation.  In  each  simulated  accident  involving 
an  aircraft  in  a  particular  size  category,  the  specific  aircraft  type  is  deter¬ 
mined  in  a  random  draw.  The  probability  that  the  aircraft  is  of  type  k  is 
determined  by  the  ratio  of  the  number  of  aircraft  of  type  k  in  the  fleet  to  the 
total  number  of  aircraft  in  the  size  category.  The  simulation  model  then  de¬ 
termines  the  amount  of  composite  material  involved  in  the  fire  by  randomly 
selecting,  on  an  equally  likely  basis,  one  of  the  accidents,  i.e.,  an  amount 
of  composite  from  the  appropriate  column  of  the  complete  form  of  Table  2.3. 


l)w.  Feller,  1950.  An  Introduction  to  Probability  Theory  and  Its  Applications, 
Vol  I,  John  Wiley,  New  York.  Page  158  et  seq. 
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1993  COMflERCIAL  AIRCRAFT  IN  SELECTFD  ACCIDENTS 


To  actually  determine  the  amount  of  carbon  fiber  released,  it  is 
assumed  that  one  percent  of  the  carbon  fiber  involved  is  actually  released  as 
3-mm  single  fibers.  Analysis  of  the  accident  data  indicated  that  approximately 
1  out  of  30  accidents  during  landings  and  takeoffs  involved  explosions.  The 
model  uses  this  value  as  the  appropriate  probability  to  randomly  determine 
whether  an  explosion  occurred  in  the  simulated  accident.  If  an  explosion  did 
occur,  an  additional  two-and-a-half  percent  of  the  carbon  fiber  is  released 
due  to  the  agitation  of  the  composite  material.  This  input  is  based  on  experi¬ 
mental  evidence  obtained  after  completion  of  Phase  I,  and  is  in  marked  contrast 
to  the  input  used  then.  In  Phase  I  the  fraction  of  carbon  fiber  released  as 
single  fibers  was  assumed  to  be  0.20. 

Each  accident  in  the  historical  file  is  also  characterized  by  the 
operational  phase  during  which  the  accident  occurred  and  the  degree  of  severity. 
The  accidents  were  analyzed  to  obtain  a  distribution  of  locations  for  landing 
and  takeoff  accidents.  The  generalized  distribution  was  applied  to  each  of 
the  airports  for  which  the  risk  calculations  were  made.  In  each  simulated 
accident  the  location  distribution  is  sampled  to  draw  an  actual  location  to  be 
used  in  the  calculation.  The  severity  measure  associated  with  each  accident 
is  used  later  in  determining  the  plume  height.  It  should  be  noted  that  the 
methods  described  here  represent  a  great  increase  in  the  amount  of  variance 
permitted  in  the  calculation  over  those  employed  in  the  Phase  I  risk  assess¬ 
ment. 
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III.  PLUME  HEIGHT  CALCULATION 


The  simulated  release  of  graphite  fibers  starts  with  the  aircraft 
accident  and  resulting  fire.  In  the  preceding  section  we  described  the 
methods  used  to  estimate  the  fraction  of  the  aircraft  consumed  in  the  fire 
and  the  calculation  of  the  amount  of  fiber  released.  The  next  step  in  the 
simulation  is  described  in  this  section. 

As  a  result  of  the  fire  a  hot  buoyant  plume  is  formed  that  rises  to 
a  "stabilization"  height  which  is  a  function  of  the  energy  available,  the  wind 
speed,  and  the  atmospheric  stability.  The  graphite  fibers  enter  the  buoyant 
plume  and  are  lifted  to  the  stabilization  height. 

GENERAL  METHOD 

As  in  Phase  I,  calculation  of  the  plume  rise  (or  elevation),  H,  at 
stabilization  from  an  open  fire  follows  the  work  of  Briggs^ since  no  improved 
approach  has  been  located  in  Phase  II.  In  the  Briggs  model,  as  adapted,  the 
height  of  the  plume,  in  meters,  is  given  by: 

H  =  2.9  (F/us)^/^ 


—  G.A.  Briggs:  "Some  Recent  Analyses  of  Plume  Rise  Observations."  Paper 
presented  at  the  1970  International  Air  Pollution  Conference  of  the  Inter¬ 
national  Union  of  Air  Pollution  Prevention  Associations. 
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for  stable  conditions,  and 

H  =  when  x  <3.5x*  (3.1  b) 

H  -  1.6F^/\'^(3.5X*)^/^,  when  x<3.5x*  (3.1  c) 

for  neutral  or  unstable  conditions,  where  u  is  the  mean  wind  speed  in  meters 
per  second  and: 

X*  =  14F^/®,  when  F<55  (3.2  a) 

X*  =  34F^/^,  when  F>55  (3.2  b) 

The  buoyancy  flux  parameter,  F,  appearing  in  the  above  equation,  is  given  by 

F  =  9Qr 

where: 

2 

G  =  acceleration  of  gravity,  9.8  m/sec 
=  heat  emission  rate,  kcal/sec 
Cp  =  specific  heat  of  air  at  constant  pressure, 

.2391  kcal/kg°K. 

O 

p  =  atmospheric  density,  1.23^  kg/m 
T  =  ambient  temperature,  °K. 

The  atmospheric  stability  parameter,  s,  is  defined  by: 

c  =  990 
^  TT 

where : 

^  A  p. 

=  gradient  of  potential  temperature,  0.35  /km 
for  stable  conditions. 

HEAT  EMISSION  RATE 

In  order  to  use  the  Briggs  formulas,  it  is  necessary  to  specify  Qp, 
the  heat  emission  rate  for  a  burning  aircraft;  this  is,  in  turn,  the  product 
of  the  rate  measured  in  gallons  per  unit  time,  and  the  fuel  heat  content  per 
gallon.  In  Phase  I  a  standard  burn  rate  was  used,  based  on  the  experimental 
data  available  at  that  time.  In  Phase  II  we  were  able  to  turn  to  the  detailed 
fire-accident  analysis  previously  referred  to.  In  this  case,  it  vias  possible 
to  estimate  the  fuel  burn  rate  for  accidents  occuring  during  different  opera¬ 
tional  phases,  as  well  as  accidents  of  different  severity.  The  reported 
accidents  involved  small  jet  aircraft  almost  exclusively,  so  a  scaling  factor 


3-2 


proportional  to  the  relative  volume  of  the  aircraft  fuel  tanks,  as  reported 
in  Janes'—'*^  was  used  to  estimate  the  burn  rates  for  other  size  aircraft.  The 
results  are  summarized  in  Table  3.1. 

Another  major  input,  or  modelling  assumption,  concerns  the  behavior 
of  the  plume  at  an  inversion.  In  the  ORI  Phase  I  Final  Report  this  matter 
was  discussed  at  some  length.  On  the  basis  of  the  evidence  available  then, 
and  not  significantly  increased  during  Phase  II,  we  continued  to  model  the 
plume  so  that  it  does  not  penetrate  the  inversion.  In  subsequent  sections  of 
this  report  we  examine  the  impact  of  this  assumption  on  the  final  results. 

With  the  inputs  described  here,  and  the  above  decision  regarding 
behavior  of  the  plume  at  an  inversion,  the  computer  model  implementing  Equa¬ 
tions  (3.1)  and  (3.2),  can  determine  the  stabilization  height  for  the  plume 
resulting  from  the  simulated  accident  involving  any  of  the  projected  1993 
aircraft  for  any  combination  of  wind  speed  and  stability  conditions. 


Janes'  All  The  Worlds  Aircraft-1977-78,  J.  Taylor. 
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TABLE  3.1 

ESTIMATED  FUEL  BURN  RATES  (Liters/Minute) 


Operational  Phase 

Damage  Severity 

Aircraft  Size 

Small 

Medium 

Large 

Takeoff 

Mi  nor 

238 

397 

794 

Substantial 

1590 

2650 

5300 

Landing 

Minor 

719 

1192 

2385 

Substantial 

1590 

2650 

5300 

Static 

- 

19 

19 

19 
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IV.  DOWNWIND  TRANSPORT  AND  DIFFUSION  OF  FIBERS 


After  a  simulated  accident,  fire,  and  release  of  carbon  fibers,  a 
buoyant  plume  carries  the  fibers  aloft,  as  described  in  the  preceding  section 
of  this  report.  The  plume  at  its  stabilization  height  may  be  considered  a 
point  source;  meteorological  transport  and  diffusion  methods  are  then  applied 
to  determine  the  downwind  dosage  (or  exposure)  at  points  of  interest. 

BASIC  CONCEPTS 

In  Phase  I,  ORI  adapted  an  essentially  standard  EPA  Gaussian  plume 
transport  and  diffusion  model  to  the  needs  of  the  risk  assessment  study.  The 
model  provides  for  downwind  transport  and  diffusion  of  material  in  the  form 
of  a  plume  that  diffuses  simultaneously  in  the  crosswind  and  vertical  directions. 
The  emitting  source  can  be  elevated  at  any  specified  height.  The  atmosphere 
is  characterized  as  being  in  one  of  several  stability  classes.  Dispersion  para¬ 
meters  that  govern  the  rate  of  crosswind  and  downwind  diffusion  are  associated 
with  each  stability  class.  The  plume  rise  calculations,  described  previously, 
give  the  source  height  which  is  used  explicitly  in  the  transport  and  diffusion 
model . 

In  Phase  II  further  extensions  were  made  to  the  ORI  transport  and 
diffusion  model.  These  allow  for  multiple  reflections  of  the  diffusing 
particles  and  provide  an  improved  mechanism  for  accounting  for  particle  fallout 
at  downwind  distances  that  are  so  large  that  the  cloud  is  uniformly  dispersed 
in  the  vertical. 
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The  wind  speed  at  plume  height  is  treated  as  representative  of  the 
layer  in  which  the  carbon  fibers  are  dispersing.  The  standard  power  law  for 
the  variation  of  wind  speed  with  height  may  be  written: 

U  =  Uq  (H/7)P  (4.1) 

where  H  is  the  height  in  meters.  In  the  cases  presented  in  this  report,  then, 

H  is  typically  assigned  a  value  equal  to  the  stabilization  height  of  the  plume 
resulting  from  the  fire  following  the  aircraft  accident.  The  exponent  "p"  is 
assigned  specific  values  for  different  atmospheric  stability  conditions,  as 
shown  i n  Tab! e  4. 1 . 

In  most  cases  rather  stringent  physical  conditions  must  be  met  for 
the  plume  to  "punch  through"  an  inversion.  Observations  indicate  that  this 
typically  does  not  occur.  It  was  therefore  considered  reasonable  to  assumie  that 
if  the  com.puted  plume  height  is  greater  than  the  height  of  the  inversion,  it 
can  be  set  equal  to  the  inversion  height.  The  impact  of  relaxing  this 
condition  and  permitting  the  plume  to  penetrate  the  inversion  is  examined 
in  Section  VII  below. 

ORI  TRANSPORT  AND  DIFFUSION  EQUATIONS 
General  Case 

The  Phase  I  meteorological  transport  and  diffusion  equations  were 
modified  to  include  the  direct  component  and  five  reflected  components  at  down¬ 
wind  locations.  The  method  of  treating  the  multiple  reflections  follows  that 
presented  by  Cramer,  et  al ^  generalized  for  additional  reflections  in 

2/ 

accordance  with  the  concepts  presented  by  the  Environmental  Protection  Agency—. 
The  general  result  for  the  dosage  (exposure)  at  a  point  on  the  surface  at  a 
location  (x,y)  in  units  of  particle -seconds  per  cubic  meter  is: 

n  1  ? 

E  (x,y,0,H' )  =  2 _  exp  [  -  -^(y  )j 

Trow  a^U  a 

^  y 


I  exp  [  -  ■^(  H_  )]  -F 
(  az 


-  H.E.  Cramer,  et  1972:  Development  of  Dosage  Models  and  Concepts,  U.S. 
Army  Dugaway  Proving  Ground,  Dugway,  Utah.  AD893  341  L. 

2/ 

User's  Manual  for  Single-Source  (CRSTER)  Model,  EPA,  July  1977,  EPA-450/ 
2-77-013. 
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TABLE  4.1 


WIND  PROFILE  EXPONENT 


Pasqui 11 -Gifford 

Stability  Class 

Exponent,  p 

A  -  Most  Unstable 

0.10 

B 

0.15 

C 

0.20 

D  -  Neutral 

0.25 

E 

0.30 

F  -  Most  Stable 

0.30 

-3 


(4.2) 


+  rexp  L  -  Y  ^  m)  1  +  [exp  -  ^  ^m) 

g  L  0  _ 

2  Z 

.  r2  exp  I  -  ^  ^  rexp  [-  ^  ('H'  ^  ^  ] 

a  a 

z  2 

+  exp  [  -  ].  (-H'  .  6  I 

z 

where: 

X  =  downwind  distance  from  source  to  receptor, 
y  =  crosswind  distance  from  source  to  receptor, 
u  =  mean  wind  speed,  m/sec, 

Q  =  number  of  particles  released 

=  Standard  deviation  of  the  wind  speed  in  the  crosswind  direction, 
as  a  function  of  x  and  the  stability  class 
=  standard  deviation  of  the  wind  speed  in  the  vertical,  as  a  func¬ 
tion  of  X  and  the  stability  class 
r  =  reflection  coefficient,  the  fraction  of  particles  that  are 

reflected  from  the  ground  surface.  The  corresponding  coefficient 
for  reflections  from  the  base  of  the  inversion  is  assumed  to  be 
un i ty . 

In  order  to  incorporate  the  effect  of  particle  fallout  into  our  cal¬ 
culations  we  adopted  the  tilted-plume  method  presented  by  Van  der  Hoven,— and 
also  used  by  Cramer—^.  Equation  (4.2)  makes  use  of  the  effective  plume  height, 
H' ,  given  by: 

H‘  =  H  -  (Vg/u)  X,  (4.3) 

where  v^  is  the  particle  fall  rate.  This  is  essentially  the  method  previously 
used  in  Phase  I  to  account  for  particle’s  falling  out  of  the  cloud. 


3/ 

Meteorology  and  Atomic  Energy  1968,  David  H.  Slade,  Editor,  AEC,  July 
1968. 

i/ 

Cramer,  cit. 
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Modification  for  Large  Distance  Downwind 

When  the  vertical  range  over  which  the  plume  is  mixed  becomes  equal 
to  the  depth  of  the  mixed  layer  (below  the  inversion),  we  can  assume  that  a 
relatively  uniform  distribution  of  particles  in  the  vertical  exists.  The 
model  therefore  makes  the  distribution  of  graphite  fibers  uniform  in  the  ver¬ 
tical,  from  the  ground  surface  to  the  base  of  the  inversion,  when  becomes 
larger  than  1.6  and  where  is  the  height  of  the  base  of  the  inversion. 

At  distances  far  enough  downwind  [o.  >1.6  H„)  that  mixing  results 

z  m 

in  an  essentially  uniform  distribution  of  the  fibers  in  the  vertical,  we  there¬ 
fore  use: 


E(x.y,0,H')=  Q 

_  Fexp  -  ■  j  (^  fl 

exp 

r  -  ''s  X  1 

2.5066a  H  u 
y  m 

L  ^  y-l 

The  general  form  of  this  expression  follows  Turner  —  t  except  for  the  final 
term,  which  accounts  for  the  fallout  of  the  particles  due  to  gravitational 
settling.  This  result  may  be  derived  by  considering  the  change  in  the  number 
of  particles  in  a  uniformly  distributed  layer  during  a  small  time  interval  of 
length  dt: 

dN  =  -(N/H^) 

where  N  is  the  number  of  particles,  and  the  other  variables  have  been  defined 
previously.  Upon  integration  we  obtain  _ 

N  =  N(0)  exp  (1-r)  tj  (4.5) 

where  N(0)  is  the  number  of  particles  present  at  time  t=0.  Since  t  may  be 
estimated  by  the  ratio  x/u,  we  obtain  the  final  term  appearing  in  Equation 
(4.4). 

INPUTS  TO  TRANSPORT  CALCULATIONS 

Actual  mixing  height  values  were  developed,  as  in  Phase  I,  from 
climatological  mean  values  reported  by  Holzworth— {  modified  for  different 


— ^D.  Bruce  Turner,  Workbook  of  Atmospheric  Dispersion  Estimates,  EPA,  1970. 
Publication  No.  AP-26. 

^olzworth,  Mixing  Heights,  Wind  Speeds,  and  Potential  for  Urban  Air  Pollution 
Throughout  the  Contiguous  United  States,  EPA,  January  1972. 


stability  conditions  as  suggested  by  Calder  — Sensitivity  tests  to  determine 
the  impact  of  changes  in  mixing  height  values  are  presented  later  in  this 
report. 

In  many  diffusion  problems  it  is  customary  to  determine  the  location 
of  an  upwind  virtual  point  source  from  which  a  diffusing  plume  would  have  grown 
to  the  size  computed  at  plume  stabilization.  In  view  of  the  large  uncertainties 
in  other  phases  of  the  risk  calculation,  and  our  concern  with  effects  some 
miles  downwind  from  the  accident  site,  we  have  set  the  virtual  point  source 
directly  over  the  accident  -  fire  site. 

The  reflection  coefficient  has  been  set  equal  to  1  at  the  inversion 
and  to  0.7  at  the  ground.  These  values  were  developed  in  consultation  with 
Messrs.  Cramer  and  Tretheway  at  a  meeting  convened  by  the  NASA  Graphite  Fiber 
Risk  Assessment  Program  Office.  The  general  association  between  the  2  centi- 
meter-per-second  fall  rate  of  the  fiber  particles  and  the  reflection  coefficient 
at  the  ground  has  been  demonstrated  by  Dumbauld,  Rafferty,  and  Cramer. 

The  diffusion  calculation  requires  input  values  of  the  dispersion 

parameters,  Oy  and  cr^,  as  functions  of  the  downwind  distance,  x,  and  the 

prevailing  stability  conditions.  The  standard  in  this  case  is  provided  by 

the  well-known  Pasqui 11 -Gifford  curves  —  .  Several  invesigators  have  questioned 

their  universal  applicability;  the  reader  is  referred  to  Pasquill's  recent  work 
9  / 

on  this  subject  —  .  In  view  of  the  fact  that  no  generally  accepted  modifica¬ 
tion  of  the  Pasqui 11 -Gifford  curves  exists,  we  adopted  these  curves  for  the 
Phase  I  calculations  and  continued  to  use  them  in  Phase  II.  For  present  pur¬ 
poses  there  were  most  conveniently  used  in  the  form  of  a  computer  program 
made  available  by  EPA. 

The  basic  weather  inputs  required,  surface  wind  speed  and  direction, 
and  stability  class,  are  drawn  from  historical  data.  These  data  were  obtained 


—  Y..  L.  Calder,  "A  Climatological  Model  for  Multiple  Source  Urban  Air  Pollu¬ 
tion,"  Appendix  D  to  A.  D.  Buse  and  J.  R.  Zimmerman,  User's  Guide  for  the 
Climatological  Dispersion  Model,  EPA-71-024,  December  1973. 


— ^See  Turner,  op  cit,  for  example. 

9/ 

—  F.  Pasqui 1 1 ,  Atmospheric  Dispersion  Parameters  in  Gaussian  Plume  Modeling, 

Part  II,  "Possible  Requirements  for  Chanae  in  the  Turner  Workbook  Values," 
EPA-600/4-76-0306. 
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from  the  National  Weather  Records  Center  for  the  airports  we  studied;  the  data 
provide  the  frequency  for  each  combination  of  the  three  weather  parameters. 

The  simulation  model  makes  a  random  draw  of  one  of  these  combinations  weighted 
by  the  input  frequency. 

TRANSPORT  AND  DIFFUSION  MODEL  SENSITIVITY  TESTS 

ORI,  Inc.  was  required  by  the  Phase  II  contract  to  test  the  model 
results'  sensitivity  to  different  particle  sizes.  For  this  purpose  the  trans¬ 
port  and  diffusion  calculations  were  performed  independently  of  the  complete 
risk  assessment  model.  Figures  4.1  to  4.4  compare  the  downwind  "foot-prints" 
for  200  kilograms  of  carbon  released  as  single  fibers  and  as  brush-clumps  in 
different  meteorological  conditions.  The  input  conditions  for  these  results 
are  summarized  in  Table  4.2. 

In  order  to  normalize  the  comparison,  the  amount  of  carbon  release  in 
the  two  forms  is  kept  constant  in  the  different  calculations.  The  combina¬ 
tion  of  the  reduced  number  of  particles  —  10^  clumps  versus  10^^  singles  —  and 
the  higher  fall  rate  of  the  clumps  results  in  a  greatly  reduced  footprint  for 
the  clumps.  The  maximum  exposures  for  the  clump  calculations  are  lower  than 
for  the  singles  by  at  least  two  orders  of  magnitude.  The  dominant  factor  in 
these  comparisons  is  the  difference  in  the  number  of  particles  per  kilogram. 

GEOGRAPHICAL  INPUTS 

The  methods  described  to  this  point  permit  the  computer  model  to  gene¬ 
rate  an  accident,  determine  a  release  amount  of  CF,  the  height  to  which  the 
carbon  fibers  are  lofted,  and  the  downwind  transport  end  diffusion  of  these 
fibers.  All  of  these  events  are  randomized  using  appropriate  Monte  Carlo 
methods.  The  transport  and  diffusion  calculation  provides  the  dosage  or 
exposure  at  particular  points  defined  by  their  by  x,  y  coordinates.  It  is 
appropriate  to  define  these  points  here,  although  much  of  the  underlying 
motivation  comes  from  the  methods  used  in  the  cost  calculations,  described 
later.  Briefly,  it  is  pointed  out  here  that  much  of  the  required  economic 
data  is  county  based.  For  this  reason  the  focus  of  our  interest  is  on 
counties  surrounding  the  airport  at  or  near  which  simulated  accidents  may 
occur. 
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TABLE  4.2 


SUMMARY  OF  INPUT  CONDITIONS 
AND  RESULTS  OBTAINED 
IN  TEST  OF  TRANSPORT/DIFFUSION 
MODEL  FOR  DIFFERENT  FIBER 
PARTICLE  SIZES 


(Release  =  200  kilograms) 


Meteorology 

Output 

Particle  Size 

Single 

Clump 

No. 

Fall  Rate 
(M/sec) 

No. 

Fall  Rate 
(M/sec) 

10^^ 

.02 

108 

1 .0 

Stability  Class  6 
Plume  Height  100  m 
Mean  Wind  5-5 
m/sec 

Fig.  No. 

4.1 

4.2 

Max  Exp. 

( FSec/m^) 

2x10^ 

2x10^ 

Stability  Class  1 
Plume  Height  1900m 

Mean  wind  3.5 

M/sec 

Fig.  No. 

4.3 

4.4 

Max  Exp. 

(FSec/m^) 

1 .6x105 

20 
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100 


100 


10  8642024  68  10 

km 


FIGURE  4.1.  DOWNWIND  EXPOSURE  PATTERN  (Fiber-seconds/cubic  meter); 
SINGLES,  STABILITY  CLASS  6. 
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FIGURE  4.2.  DOWNWIND  EXPOSURE  PATTERN  (Fiber-seconds/cubic  meter)* 
CLUMPS,  STABILITY  CLASS  6. 
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100 


As  in  Phase  I,  then,  county-based  economic  data  were  adopted  for 
computer  input;  in  many  cases  counties  were  divided  into  smaller,  homogeneous 
geographical  units.  In  each  case  the  center  of  the  county  or  sub-county  geo¬ 
graphical  unit  was  selected  and  a  representative  circle  inscribed  within  that 
area.  The  input  data  set  includes  the  coordinates  of  the  center  and  the 
associated  radius.  The  exposure  and  resulting  impact  calculations  are  made 
at  the  center  and  points  a  distance  equal  to  two-thirds  of  the  radius  to  the 
east,  west,  north,  and  south  of  the  center. 

Figure  4.5  shows  this  geometrical  pattern  schematically.  This 
method  was  adopted  to  provide  area-sensitivity  in  the  resulting  impact  cal¬ 
culation.  The  use  of  the  two- thirds  radius  mesh  interval  was  selected  so  that 
representative  points  selected  in  neighboring  circles  could  not  be  colocated. 
The  resulting  mesh,  if  all  circles  were  equal  in  size,  would  be  square  with 
all  points  equidistant  from  one  another. 

In  each  case  the  county-based  business/industry  sites  are  uniformly 
distributed  over  these  five  points.  The  concept  is  illustrated  in  Figure  4.6, 
as  it  was  applied  to  one  county  for  the  Washington  National  Airport  risk  cal¬ 
culations.  In  all  cases  this  method  was  applied  to  the  area  around  each  air¬ 
port  to  a  distance  of  50  miles  or  more. 
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FIGURE  4.6. 


Definition  of  Areas  at  Risk  for  Washington  National 
Airport.  Howard  County,  Maryland,  Outlined  in  Upper 
Map,  Shown  in  Detail  in  Lower  Map.  Circles  Represent  Con¬ 
centrations  of  Business,  Industry,  and  Residences.  Method 
was  Applied  to  all  Counties  Within  50  Miles  of  Airport. 
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V.  TRANSFER  OF  FIBERS  INTO  INTERIOR  OF  STRUCTURES 


In  computing  the  impact  on  electrical  and  electronic  equipment  of 
exposure  to  carbon  fibers  we  are  principally  concerned  with  equipment  inside 
buildings.  It  is  therefore  necessary  to  compute  the  exposure  inside  a  building 
resulting  from  a  known  value  of  the  exposure  outside.  This  section  of  the 
report  describes  the  methods  developed  by  ORI  to  accomplish  this.  In  the 
logic  flow  of  the  risk  assessment  simulation  model  the  transfer  calculation 
follows  the  computation  of  exterior  exposure  values. 

METHOD 

When  a  building  is  impinged  on  by  a  plume  of  carbon  fibers,  some 
of  the  fibers  may  enter  the  building  through  air  conditioning  or  other  ventila¬ 
tion  systems  and  by  various  leakage  paths.  Once  inside  the  building  or 
enclosure,  fibers  will  be  removed  by  fallout  and  through  leakage  paths  back 
to  the  outside.  If  inside  air  is  recirculated  and  filtered,  additional 
fibers  will  be  removed.  The  concentration  of  fibers  that  produce  failure 
stresses  on  equipments  in  a  building  or  enclosure  at  any  time  may  be  deter¬ 
mined  from  equations  describing  the  net  flow.  These  have  been  developed  in 
a  relatively  simple  form  by  Slade. 


—  David  H.  Slade,  Editor,  Meteorology  and  Atomic  Energy,  AEC,  July  1968. 
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In  Phase  I,  ORI  was  able  to  show  that  the  "transfer  function"  or 
ratio  of  interior  to  exterior  exposure  can  be  expressed  as: 


I  -Jh _ 

0 


(5.1) 


where: 


v.j  =  rate  at  which  fiber-borne  air  enters  the  building,  or  enclosure 
through  the  air  conditioning  system  and  all  leakage  paths 
Vp  =  rate  at  which  fiber-borne  air  leaves  the  building,  including 
that  removed  by  recirculation 
Vg  =  fall  rate  of  carbon  fibers 

v^  =  rate  at  which  fibers  are  removed  by  recirculation  filtering 
s  =  volume  of  building  or  enclosure 
a  =  area  of  space  subject  to  fallout. 


IMPLEMENTATION 


As  in  the  earlier  Phase  I  effort.  Equation  (5.1)  provided  the  basis 
for  calculating  interior  exposure  values.  In  Phase  I,  ORI,  Inc.  defined 
several  types  of  buildings  and  other  enclosures;  each  was  characterized  by 
size,  types  of  doors  and  windows  and  ventilating  equipment.  These  basic 
enclosure  types  were  used  with  some  minor  revisions  in  Phase  II.  It  was  no 
longer  necessary  to  treat  equipment  enclosures  explicitly  since  all  equip¬ 
ment  failure  tests  (see  Section  VI  below)  included  the  effect  of  typical 
enclosures.  The  following  principal  building/enclosure  categories  were 
defined  in  Phase  II: 

1.  Small  Equipment  Building  or  Van 

2.  Medium  Equipment  Building 

3.  Large  Equipment  Building  or  Factory 

4.  Equipment  Room  in  Building 

5.  Utility  Room 

a)  filtered 

b)  unfiltered 

6.  Residence 

a)  air  conditioned 

b)  not  air  conditioned 

7.  Retail/Wholesale  Establishments. 
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Generalized  design  factors  are  associated  with  each  of  these  building/enclo¬ 
sure  types  in  Table  5.1. 

It  was  shown  in  Phase  I,  that,  as  long  as  basic  architectural  rela¬ 
tionships  are  maintained,  the  ventilation  mode  of  a  building  is  essentially 
independent  of  the  actual  size  of  the  building.  These  design  factors  are 
used  to  determine  the  air  conditioning  flow  rates,  filter  efficiencies,  and 

air  leakage  rates  used  in  Equation  (5.1).  Ventilation  rates  were  based  on 

2  3/ 

published  industry  standards-^—  .  The  values  of  filter  efficiency  used  in 
Phase  I  were  changed  to  incorporate  new  experimental  results.  It  was  also 
shown  in  Phase  I  that,  although  ventilation  rates  are  a  function  of  wind 
speed,  the  "fallout  term"  in  Equation  (5.1)  tends  to  be  dominant.  Accordingly, 
transfer  functions  were  computed  for  a  nominal  10  meter-per-second  wind 
speed.  The  resulting  transfer  functions,  shown  in  Table  5.2,  were  used  in 
all  Phase  II  calculations.  Specific  building  types  were  associated  with 
different  categories  of  business  and  industry,  as  described  in  Section  VI, 
below. 


2/ 

—  Carrier  Air  Conditioning  C.,  Handbook  of  Air  Conditioning  System  Design, 
McGraw-Hill  Book  Co..  1965  - 

3/  Baumeister  &  Marks,  Standard  Handbook  for  Mechanical  Engineers,  McGraw- 
Hill  Book  Co.,  1967.  - - 
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TABLE  5.2  -  TRANSFER  FUNCTIONS  FOR  STANDARD  ENCLOSURES 


Enclosure 

Category 

Transfer 

Function 

.012 

2.  Medium  Equipment  Building 

.010 

3.  Large  Equipment  Building  or  Factory 

Building  (per  floor) 

.004 

4.  Equipment  Room  in  Building  (one  exterior  wall) 

.010 

5.  Utility  Room 

F i 1 tered 

.023 

Non-Filtered 

.094 

6.  Residence 

Air  Conditioned 

.058 

Not  air  Conditioned 

.004 

7.  Retail /Wholesale  Establishment 

.004 
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VI.  EQUIPMENT  FAILURES 


FAILURE  MODEL 

The  probability  of  failure  of  equipment  which  is  exposed  to  car¬ 
bon  fibers  is  obtained  from  the  exponential  expression: 

Pp  =  1  -  exp  (-E/E)  (6.1) 

where: 

Pp  =  probability  of  equipment  failure 

E  =  exposure  level  in  the  immediate  vicinity  of  the 
vulnerable  equipment,  in  fiber-seconds  per  cubic 
meter 

F  =  average  exposure  causing  a  failure. 

During  Phase  I,  the  U.S.  Army  Ballistics  Research  Laboratory  (BRL)  at  Aber¬ 
deen,  Maryland,  determined  that  experimental  failure  data  for  many  classes  of 
equipment  fit  an  exponential  failure  Later,  it  was  shown  that 

certain  failures  were  mutl iple-fiber  events.  It  appeared  that  the  generalized 
Weibull  distribution  provided  a  better  fit  to  failure  data  for  those  equip¬ 
ments.  In  Phase  II  it  has  been  shov/n  that,  even  for  those  equipments  whose 


Shelton  and  Moore,  Have  Name  Vulnerability  of  the  Improved  Hawk  System,  BRL 
Report  No.  1964,  February  1977. 

2_/ 

ORI  discussions  with  BRL,  August  15,  1978. 
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failures  do  not  obey  the  exponential  law,  it  is  conservative  to  use  the 
exponential  law  in  estimating  failures.  The  exponential  relationship  gives 
a  higher  value  of  the  failure  probability  for  low  values  of  the  exposure  than 
the  Weibull  distribution,  thus  overestimating  failures,  and  providing  the 
desired  conservatism  in  estimating  the  overall  risk.  Typical  values  of  the 
exponential  failure  parameter  for  generic  equipment  types  are  shown  below  in 
Table  6.1.  It  should  be  noted  that  the  failure  concepts  developed  here  apply 
only  to  equipment  when  it  is  energized. 

The  exposure  used  in  Equation  (6.1)  is  that  directly  impinging  on 
the  vulnerable  equipment.  When  this  equipment  is  inside  a  building,  the 
interior  exposure  may  be  obtained  from  the  exterior  exposure  by  multiplying 
the  exterior  exposure  by  the  appropriate  transfer  function  (TF),  as  described 
in  Section  V,  above.  Since  the  transfer  function  and  the  mean  exposure  to 
failure,  F,  are  constants  for  a  particular  piece  of  equipment  in  a  particular 
building,  we  define  a  failure  parameter: 

K,J  =  (TF)yr,  (6.2) 

where: 

K.jj  =  overall  failure  parameter  for  equipment  of  type  i 
in  a  building  of  type  j 

(TF)j  =  penetration  factor  (transfer  function)  for  a  building  of 
type  j 

E.j  =  mean  exposure  to  failure  for  equipment  of  type  i. 

In  subsequent  applications,  the  parameter  is  substituted  into  Equation 
(6.1)  to  give  the  probability  of  failure  for  equipment  of  type  i  in  a  building 
of  type  j  for  any  exterior  exposure: 

Pp_,j  =  1  -  exp  (-K,jE^)  (6.3) 

Thus,  although  the  exterior-to-interior  transfer  process  has  been  discussed 
as  a  separate  entity  in  the  preceding  section,  we  were  able  to  combine  the 
failure  and  transfer  calculations  in  one  procedural  step  by  defining  specific 
equipments  in  specific  types  of  buildings.  These  methods  are  described  in  more 
detail  below. 
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EQUIPMENT  CONFIGURATIONS 

In  treating  typical  equipment  configurations  it  is  convenient  to 
develop  expressions  for  the  collective  probability  of  failure  of  the  complete 
configuration.  In  particular,  if  n  identical  equipments  are  in  series  so  that 
a  failure  of  one  causes  the  entire  "line"  to  fail,  the  probability  that  the 
line  fails  is: 

Pp  (LINE)  =  1- 

-  (6.4) 

Similarly  if  n  like  equipments  are  in  parallel,  so  that  the  operation  fails 
only  if  all  equipments  fail,  the  aggregate  probability  of  failure  is: 

Pp(Operatlon)  =  P^^pj  (P  5, 

The  computer  program  that  determines  the  impact  of  each  simulated  aircraft 
accident  and  associated  release  of  graphite  fibers  uses  Equations  (6.3)-(6.5) 
to  estimate  the  probability  that  each  business  or  industry  in  the  geographical 
area  of  interest  is  affected. 

One  of  the  major  efforts  in  Phase  I  was  the  characterization  of  each 
business-industry  sector,  defined  by  an  SIC  (Standard  Industrial  Classifica¬ 
tion)  number,  by  a  specific  set  of  equipments  installed  in  a  specific  type  of 
building.  This  effort  was  extended  and  made  more  detailed  in  Phase  II.  The 
generalized  business/industry  equipment  configuration  showing  the  electric 
power  flow  appears  in  Figure  6.1;  in  any  one  class  of  business  or  industry 
portions  of  this  configuration  may  not  be  present.  Typical  individual  equip¬ 
ments  in  each  of  the  modules  shown  in  Figure  6.1  are  defined  in  Table  6.1,  with 
their  estimated  values  of  mean  exposure  to  failure. 

The  equipment  configuration  was  made  specific  to  plants  of  different 
size  (small,  medium,  and  large)  in  each  pertinent  SIC  -  number  category.  An 
example  will  illustrate  the  method.  A  large  plant  in  Category  28A  (comprising 
all  3-digit  SIC  code  numbers  under  28,  basically  chemical  and  allied  products) 
has  an  internal  power  interface  characterized  by  one  set  of  input  power  service 
equipment,  one  distribution  panel,  and  an  auxiliary  generator.  Its  common 
module  consists  of  two  computers  in  parallel  and  two  keyboard  display  units 
in  parallel.  The  plant  has  25  lines  in  its  distributed  module.  Each  line 
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POWER  MODULE 


a.  CL  LL  3 
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TABLE  6.1 


GENERIC  BUSINESS/INDUSTRY  EQUIPMENTS  WITH  MEAN  EXPOSURE  TO  FAILURE 
VALUES  (E  IN  FIBER  SECONDS/METER-^) 


Equipment 

Failure 

Module 

Code 

Definition 

rcirame  L6r 

(E) 

Power 

SW 

' 

Input  power  service  equipment  -  trans-  ^ 

formers,  breakers,  switchgear 

10^ 

DIST 

Power  distribution  buses  and  panels 

10^ 

AUX 

Auxiliary  power  supply  in  parallel 
with  power  input 

10^ 

Common 

COMP 

Standard-size  computer  used  as  a  cen¬ 
tral  facility  controller 

10^ 

K/D 

Keyboard-display  unit 

10^ 

Distri¬ 

buted 

PS 

High-voltage  power  supply  at  a  machine 
station 

10^ 

I  NT 

Interface  unit  used  to  buffer  cen¬ 
tral  computers  to  line  controllers 

10« 

MC 

Manual  controller,  associated  with  each 
electrically-operated  machine 

10^ 

MPC 

Mini -computer  used  as  a  programmable 
controller 

10^ 

liPC 

Microprocessor  used  as  a  controller 

10^ 

MM 

High-voltage  motor  controller 

10^ 

MS 

Machine  station  servo-mechanism 

10® 

MH 

Heater  or  oven  control 

10® 

SENSOR 

Device  to  measure  temperature, 
thickness,  weight,  position, 
motion,  etc. 

10^ 
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consists  of: 


•  5  high-voltage  power  supply  units 

•  5  interface  units 

•  5  manual  controllers 

0  5  minicomputers,  used  as  controllers 

•  2  high-voltage  motor  controllers 

«  2  machine  station  servo-mechanisms 

•  1  heater  control  unit 

•  5  sensor  units. 


Similar  configurations  were  defined  for  all  vulnerable  categories  of  business 
and  industry.  The  data  was  developed  as  a  result  of  an  extensive  literature 
search,  augmented  by  site  visits  during  Phases  I  and  II. 

The  data  collection  effort  during  ORI's  risk  assessment  contract 
included  visits  to  one  or  more  plants  in  each  of  the  following  major  categories 


•  2011  -  Meatpacking 

•  2331  -  Womens  Blouses 

•  262  -  Paper  Mills 

t  2721  -  Periodicals 

•  2732  -  Book  Printing 

•  3519  -  Internal  Combustion  Engine  Manufacturing 

•  3661  -  Telephone  and  Telegraph  Equipment 

•  3662  -  Electronic  Equipment 

•  458  -  Air  Transportation  Services 

9  481  -  Radio  and  Television  Broadcasting 

•  491  -  Electric  Services 

e  806  -  Hospitals 


The  results  of  these  site  visits,  conferences  with  NASA  personnel,  and  the 
earlier  literature  surveys  are  summarized  in  Table  6.2.  Typically,  a  large 
factory  has  more  than  250  employees,  a  medium  size  factory  50  to  249  employees, 
and  a  small  factory  20  to  49  employees. 
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TABLE  6.2 

EQUIPMENT  CONFIGURATIONS  FOR  MANUFACTURING 
FACILITIES  BY  SIC  GROUP  AND  SIZE 
(NO.  OF  EQUIPMENTS) 


SiC 

GROUP/ 

SIZE 

POWER 

MODULE 

COMMON 

MODULE 

DISTRIBUTED  MODULE 

SW 

DIST 

K/D 

COMP 

K/D 

NO.  OF  LINES 

PS 

INT 

MC 

MPC 

^PC 

K/D 

MM 

MS 

MH 

SENSOR 

20A 

L 

1 

1 

1 

0 

0 

25 

0 

0 

5 

0 

0 

0 

0 

0 

0 

0 

M 

1 

1 

0 

0 

0 

10 

0 

0 

5 

0 

0 

0* 

0 

0 

0 

0 

S 

0 

1 

0 

0 

0 

1 

0 

0 

5 

0 

0 

0 

0 

0 

0 

0 

20B 

L 

1 

1 

1 

(2) 

5 

1 

12 

12 

12 

0 

0 

0 

0 

0 

0 

M 

1 

1 

0 

■■ 

0 

3 

1 

1 

12 

1 

12 

1 

0 

5 

1 

12 

S 

0 

1 

0 

0 

0 

1 

0 

0 

12 

0 

0 

0 

0 

0 

0 

0 

21A 

L 

1 

1 

1 

0 

0 

25 

0 

0 

5 

0 

1 

0 

0 

1 

0 

1 

M 

1 

1 

0 

0 

0 

10 

0 

0 

5 

0 

0 

0 

0 

0 

0 

0 

S 

0 

1 

0 

0 

0 

1 

0 

0 

5 

0 

0 

0 

0 

0 

0 

0 

22A 

L 

1 

1 

1 

(2) 

(2) 

50 

1 

6 

6 

6 

0 

1 

3 

0 

6 

M 

1 

1 

0 

0 

0 

10 

1 

1 

6 

1 

6 

1 

1 

3 

0 

6 

S 

0 

1 

0 

0 

0 

5 

0 

0 

6 

0 

0 

0 

1 

0 

0 

0 

23A 

L 

1 

1 

1 

0 

0 

75 

0 

0 

4 

0 

2 

0 

0 

0 

0 

2 

M 

0 

1 

0 

0 

0 

25 

0 

0 

4 

0 

2 

0 

0 

0 

0 

2 

S 

0 

1 

0 

0 

0 

5 

0 

0 

4 

0 

0 

0 

0 

0 

0 

0 

24A 

L 

1 

1 

1 

0 

0 

10 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

M 

1 

1 

0 

0 

0 

5 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

S 

0 

1 

0 

0 

0 

1 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

25A 

L 

1 

1 

1 

0 

0 

10 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

M 

1 

1 

0 

0 

0 

5 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

S 

0 

1 

0 

0 

0 

1 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

26A 

L 

1 

1 

1 

(2) 

(2) 

3 

4 

8 

8 

8 

3 

0 

4 

3 

1 

4 

M 

1 

1 

0 

0 

0 

3 

2 

0 

6 

1 

5 

1 

3 

3 

0 

2 

28A 

L 

1 

1 

1 

(2) 

(2) 

25 

5 

5 

5 

5 

0 

2 

2 

1 

5 

M 

1 

1 

0 

0 

0 

10 

5 

1 

5 

1 

5 

1 

2 

2 

1 

5 

S 

0 

1 

0 

0 

0 

1 

0 

0 

5 

0 

0 

0 

2 

0 

0 

0 

29A 

L 

1 

1 

1 

(2) 

(2) 

10 

5 

1 

8 

1 

8 

0 

8 

0 

0 

5 

30A 

L 

■■ 

(2) 

(2) 

25 

5 

5 

5 

5 

0 

0 

2 

2 

1 

5 

M 

■■ 

0 

0 

10 

5 

1 

5 

1 

5 

1 

2 

2 

1 

5 

S 

0 

0 

1 

0 

0 

5 

0 

0 

0 

2 

0 

0 

0 

32A 

L 

1 

1 

1 

0 

0 

5 

2 

1 

8 

6 

0 

1 

0 

5 

2 

5 

33A 

L 

1 

1 

1 

(2) 

(2) 

5 

2 

1 

8 

1 

6 

1 

5 

0 

0 

5 

35A 

L 

1 

1 

1 

50 

1 

1 

5 

1 

5 

1 

0 

5 

0 

2 

!VI 

1 

1 

0 

20 

1 

0 

5 

0 

5 

0 

0 

5 

0 

2 

S 

0 

1 

0 

0 

0 

5 

0 

0 

5 

0 

0 

0 

0 

0 

0 

0 

35B 

L 

1 

1 

0 

0 

0 

5 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

36A 

L 

1 

1 

1 

mam 

(2) 

50 

1 

1 

5 

1 

5 

2 

2 

0 

2 

M 

1 

1 

0 

■■ 

0 

20 

1 

0 

5 

0 

5 

0 

2 

2 

0 

2 

S 

0 

1 

0 

0 

0 

5 

0 

0 

5 

0 

0 

0 

2 

0 

0 

0 

1 

1 

1 

(2) 

(2) 

50 

mm 

1 

5 

0 

0 

3 

1 

2 

1 

1 

0 

0 

0 

20 

0 

5 

0 

0 

3 

1 

2 

0 

1 

0 

0 

0 

5 

0 

0 

0 

0 

0 

0 

0 

36C 

L 

1 

1 

0 

0 

0 

6 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

37A 

L 

1 

1 

1 

(2) 

(2) 

50 

2 

1 

8 

1 

8 

1 

2 

4 

1 

5 

1 

1 

1 

0 

0 

25 

1 

1 

5 

1 

5 

1 

0 

4 

1 

5 

1 

1 

0 

0 

0 

10 

1 

0 

5 

0 

5 

0 

0 

4 

1 

5 

0 

1 

0 

0 

0 

5 

0 

0 

5 

0 

0 

0 

0 

0 

0 

0 

(  )  DENOTES  EQUIPMENTS  IN  PARALLEL 

SIC  GROUPS  ARE  DEFINED  IN  NOTES  FOLLOWING 
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Components  and  Assessories 

36C  3662  Radio  and  TV  Communication  Equipment 

37  A  37  X  Transportation  Equipment 

38  A  38  X  Instruments  and  Related  Products 


The  description  above  provides  the  linkage  from  SIC  number  and 
size  to  equipment  configuration,  then  to  specific  equipments  and  their 
associated  failure  parameters.  It  was  also  necessary,  as  described  above  in 
Section  V,  to  relate  specific  building  types  to  each  vulnerable  class  of 
business  and  industry.  These  results  are  summarized  in  Table  6.3.  The 
table  associates  the  different  building  types  defined  in  Table  5.1  (with 
transfer  functions  in  Table  5.2)  with  each  of  the  major  sections  (modules) 
of  plants  of  different  sizes  in  different  SIC  groups. 

COMPUTER  IMPLEMENTATION 

The  mean  exposure-to-failure  values  for  the  generic  equipments 
defined  above  were  summarized  in  Table  6.1.  In  using  these  inputs  the  equip¬ 
ment-specific  value  of  F  was  combined  with  the  building-specific  transfer 
function,  in  accordance  with  Equation  (6.2).  In  order  to  estimate  the  impact 
on  specific  business  and  industrial  complexes  it  was  assumed  that  the  plant 
is  down  if  electric  power  is  lost  inside  the  plant,  if  the  common  module 
fails,  or  if  more  than  one  half  of  the  ‘‘lines"  in  the  distributed  module  fail. 
The  implementation  of  these  modeling  concepts  is  described  in  more  detail  in 
the  following  section  of  the  report. 

Phase  II  results  reported  by  other  investigators  indicated  that  the 
high-voltage  power  supply  system  is  essentially  invulnerable;  it  was  assumed 
that  an  equivalent  piece  of  equipment  representing  the  bushings  and 
bus  of  a  step-down  transformer  could  be  used  to  represent  the  possibility  of 
an  exterior  power  supply  failure. 
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TABLE  6.3 

ENCLOSURE  TYPES  BY  SIC/SIZE  CATEGORY 
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VII.  COSTS  DUE  TO  EQUIPMENT  FAILURES 


This  section  of  the  report  presents  ORI's  Phase  II  methodology  for 
determining  the  costs  associated  with  equipment  failures.  The  most  significant 
changes  to  the  Phase  I  methodology  were  introduced  in  this  part  of  the  risk 
assessment  calculation.  Three  categories  of  cost  were  considered  for  business 
and  industry  impacts: 

•  Repair  of  damaged  electrical  equipment 

•  Facility  cleanup 

•  Business/industry  disruption. 

In  the  Phase  I  risk  assessment,  attention  was  focussed  on  the  latter  cost 
category  using  an  expected  value  technique.  In  Phase  II  the  model  has  been 
expanded  to  treat  all  the  above  categories  explicitly,  while  disruption  costs 
are  now  computed  using  a  Monte  Carlo  random  process.  Household  equipment 
failures  are  treated  as  in  Phase  I,  using  an  expected-value  algorithm.  A 
completely  new  submodel  has  been  developed  to  compute  the  cost  incurred  as 
a  result  of  failures  of  avionics  equipment  aboard  commercial  aircraft  on  the 
ground. 
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BUSINESS/INDUSTRY  REPAIR  COSTS 


For  each  of  the  generic  types  of  equipment  defined  previously  (cf. 
Table  6J)  a  repair  cost  was  estimatedj  using  data  provided  by  the  Ballistics 
Research  Laboratoryj  information  gained  on  the  industrial  site  visits,  and 
other  sources.  Equipment  repair  cost  inputs  are  summarized  in  Table  7.1. 

In  many  cases  it  v/ould  be  expected  that  repairs  could  be  effected  by  the 
simple  act  of  vacuum  cleaning  the  equipment  that  failed.  It  was  assumed, 
hov/fiver,  that  a  minimum  repair  cost  would  still  be  incurred  to  cover  trouble¬ 
shooting  and  repair  time  on  the  basis  that  equipment  users  would  not  usually 
ba  av/are  of  this  fact. 

As  shown  in  the  preceding  section  each  business  or  industrial  facil- 
'ity  defined  by  SIC  number  and  size  has  a 'defined  equipment  "suit"  (cf.  Table 
b.?).  The  computer  model  treats  all  the  equipments  of  one  type  at  one  geo- 
g':"&ph1ce,l  location  collectively.  It  first  computes 


N(i)  2  (Equipments  of  Type  i)  SIC  S 


SIC  S 


(7.1) 


to  obtain  the  total  number  of  equipment  of  type  1  at  the  location.  At  its 
most  straightforward  the  simulation  would  have  been  written  to  test  each  of 
the  M(1)  equipments,  and  determine  whether  each  failed  using  a  procedure  that 
ccmpares  a  random  number  with  the  computed  failure  probability  Pp(i).  This 
prccedure  Is  easy  to  program,  but  is  somev/hat  inefficient  and  wastes  computer 
tlnia  if  thare  are  many  pieces  of  equipment.  Since  each  piece  of  equipment 
either  fails  or  not  the  process  is  an  example  of  a  Bernoulli  trial.  The  pro- 
ability  that  exactly  k  equipments  fail  is  given  by 


bik;  N(i)s  Pp(i)] 


(7.2) 


Further.,  if  H(i)  is  relatively  large,  and  Pp(i)  is  small,  which  is  true  for 
the  cases  of  interest  here,  and  we  define 


Np(i)  =  N(1)P 


F' 


(i): 


(7.3) 


number  of  failures  of  equipment  of  type  i  at  the  particular  loca¬ 


tion,  then  is--^  : 


b  \  ki  N(1),  Pp(i) 


%  0“^p( f( i ) 


k! 


(7.4) 


Feller,  od. 


cit. ,  p .  1 48  et  seg. 
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Table  7.1 

INPUT  REPAIR  COSTS  FOR  ORI  STANDARD  EQUIPMENTS 


Equipment 

CODE 

Definition 

Repair  Cost  ($) 

POW 

Exterior  Step-Down  Transformer 

300 

SW 

Input  power  service  equipment  -  trans¬ 
formers,  breakers,  switchgear 

3,000 

DIST 

Power  distribution  buses  and  panels 

2,600 

AUX 

Auxiliary  power  supply  in  parallel 
with  power  input 

5,000 

COMP 

Standard-size  computer  used  as  a  central 
facility  controller 

50,000 

K/D 

Keyboard-display  unit 

3,000 

PS 

High-voltage  power  supply  at  a  machine 
station 

2,000 

INT 

Interface  unit  used  to  buffer  central  com¬ 
puters  to  line  controllers 

600 

MC 

Manual  controller,  associated  with  each 
electrically-operated  machine 

2,500 

MPC 

Mini -computer  used  as  a  programmable  con¬ 
troller 

10,000 

juPC 

Microprocessor  used  as  a  controller 

7,000 

MM 

High-voltage  motor  controller 

5,600 

MS 

Machine  station  servo-mechanism 

1,000 

MH 

Heater  or  oven  control 

1,000 

SENSOR 

Device  to  measure  temperature,  thickness, 
weight,  position,  motion,  etc. 

6,000 

-3 


Equation  (7.4)  Is  the  Poisson  distribution  with  mean  equal  to  the  expected  num¬ 
ber  of  failures.  In  performing  the  simulation  the  model  computes  the  expected 
number  of  failures  for  each  class  of  equipment  in  turn.  The  number  of  failures 
is  then  obtained  by  drawing  a  random  sample  from  the  appropriate  Poisson  dis¬ 
tribution.  This  method  is  essentially  equivalent  to  "playing"  the  failure  of 
each  equipment  individually  but  is  much  more  economical. 

The  procedure  described  above  is  used  in  several  places  in  the  cal¬ 
culation,  because  of  the  simplification  and  economy  it  introduces  into  the 
calculation,  with  only  very  little  loss  in  generality.  In  cases  where  the 
same  class  of  equipment  is  located  in  facilities  with  different  transfer  func¬ 
tions  they  are  treated  as  different  equipments  types  for  computational  pur¬ 
poses.  Once  the  number  of  failures,  Np(i),  is  obtained  by  sampling  the  Poisson 
distribution,  the  total  repair  cost  for  that  equipment  type  is  the  product  of 
the  repair  cost  per  equipment  (Table  7.1)  and  Np(i).  This  is  repeated  for  all 
types  of  equipment  at  a  given  location  in  the  downwind  path  of  the  plume.  The 
computer  program  logic  is  illustrated  in  Figure  7.1,  which  is  a  schematic  flow 
chart  for  this  calculation. 

FACILITY  CLEAN-UP  COSTS 

Estimates  of  facility  cleanup  costs  were  made  for  different  businesses 
and  industries  on  the  basis  of  type  of  business  and  size  of  plant.  Using  infor¬ 
mation  gained  during  the  Phase  II  site  visits  it  was  estimated  that  the  decision 
to  institute  a  special  plant-wide  cleanup  would  be  made  on  the  basis  of  evidence 
of  major  impact  of  the  presence  of  carbon  fibers.  Accordingly,  it  is  assumed 
that  an  intensive  plant  cleanup  is  implemented  whenever  the  plant  is  shut  down 
due  to  equipment  failures,  as  described  below.  For  each  plant  or  other  facility 
that  is  shut  down  the  model  looks  up  the  input  cleanup  cost  for  a  plant  of 
that  SIC  number-size  combination.  The  calculation  of  plant  shut  down  is 
described  below. 

DISLOCATION  COST 

It  was  assumed  that  a  plant  or  place  of  business  would  be  shut  down 
if  power  were  lost,  the  common  module  failed,  or  more  than  half  of  the  produc¬ 
tion  lines  failed.  Figure  7.2  illustrates  this  concept  in  a  decision  tree 
formulation.  The  computation  is  done  for  all  plants  in  one  SIC-code  number 
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FIGURE  7. 


1  FLOWCHART  FOR  COMPUTING  REPAIR  COSTS  DUE  TO  EQUIPMENT  FAILURE 
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group  at  one  location.  In  contrast  to  Phase  I,  then,  we  determined  plant 
closings  on  a  stochastic  basis,  rather  than  employing  an  expected-value 
algorithm.  The  shutdown  calculation  proceeds  through  each  SIC  number  and 
each  size  group  with  that  SIC  number  at  each  geographical  location. 

Individual  Module  Failures 

The  probability  that  the  power  module  fails  is  determined  by  com¬ 
puting  the  probability  of  failure  of  the  primary  power  input,  the  switch  gear, 
the  distribution  panel,  and  auxiliary  power,  if  present.  The  probability  of 
a  power  failure  ahead  of  the  distribution  panel  may  be  expressed  as 

Pp(Power  In)  =  {l-{l-Pp(P0W)}  {l-Pp(SW)}}  Pp(AL)X)  (7.5) 

where  Pp(P0W)5  Pp(SW),  and  Pp(AUX)  are  the  computed  failure  probabilities  for 
the  primary  pov/er  source,  the  switch  gear,  and  the  auxiliary  power  system, 
respectively.  The  probability  that  the  plant  is  without  power  is  then  estimated 

by: 

Pp(Power)  =  Pp(Power  In)  +  {l-Pp(Power  In)}  Pp(DIST)  (7.6) 

where  Pp(DIST)  is  the  probability  of  a  failure  at  the  distribution  panel. 

The  probability  that  the  common  module  fails  is  estimated  by: 

Pp(Common)  =  1  -  {l-Pp(COMP)}  {l-Pp(K/D)}  (7.7) 

where  Pp(COMP)  and  Pp(K/D)  are  the  failure  probabilities  for  the  computer  and 
keyboard  displays  respectively,  and  n  and  m  are  the  numbers  of  each  in  parallel. 

The  probability  that  one  line  in  the  distributed  module  fails  is  given 

by: 

Pp(Line)  =  1  -  {1-Pp(i )  }"^^' ^  {1-Pp(  j  ^  .  (7.8) 

where  Pp(i)  is  the  probability  of  failure  for  equipment  of  type  i,  and  n(i) 
is  the  number  of  units  of  type  i  in  series  in  the  line.  Equation  (7.8)  indi¬ 
cates  that  the  line  fails  if  at  least  one  unit  in  series  in  the  line  fails;  in 
the  equation  we  have  indicated  that  there  are  n  units  of  type  i  and  m  units  of 
type  j  in  the  line. 
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DOES  POWER  MODULE  FAIL? 
YES  NO 

I 

PLANT 
FAILS 


DOES  COMMON  MODULE  FAIL? 
YES  NO 

PLANT 

FAILS 


DO  MORE  THAN  HALF  OF 
LINES  IN  DISTRIBUTED 
MODULE  FAIL? 

YES  NO 

I 

PLANT 
FAILS 

PLANT 

SURVIVES 


FIGURE  7.2.  SCHEMATIC  PLANT  FAILURE  DECISION  TREE 


Computational  Method 


Rather  than  examine  each  plant  individually,  the  model  examines  the 
group  of  plants  in  the  SIC  number-size  group.  The  logic  is  illustrated  sche¬ 
matically  in  a  flow  chart  appearing  as  Figure  7.3.  We  first  determine  the 
expected  number  of  power  module  failures  by  multiplying  the  number  of  plants, 
N(SIC,  Size),  by  Pp(Power).  The  actual  number  of  power  module  failures  is 
drawn  from  a  Poisson  distribution  with  this  mean  value.  This  method  is  entirely 
analagous  to  that  derived  above  for  the  equipment  failures. 

Next,  this  submodel  treats  the  surviving  plants,  those  of  the  orignal 
N(SIC,Size)  that  did  not  suffer  power  module  failures.  The  model  samples  a 
Poisson  distribution  with  mean  equal  to  the  product  of  the  number  of  survivors 
and  Pp(Common)  to  determine  the  number  of  facilities  that  fail  due  to  failures 
of  the  common  module. 

Those  plants  that  survive  the  power  and  common  module  "cuts"  are 
then  examined  one  by  one.  For  each  of  these  plants  we  determine,  again  by 
sampling  a  Poisson  distribution,  the  number  of  lines  that  fail.  The  expected 
number  (or  mean  of  the  distribution)  is  the  product  of  the  number  of  lines  in 
the  distributed  module  and  Pp(Line).  For  each  plant  the  computer  program  de¬ 
termines  whether  the  randomly  generated  number  of  lines  that  fail  is  equal  to 
half  or  more  of  all  the  lines  in  the  plant.  If  so,  the  plant  is  counted  as 
"failed."  This  is  repeated  for  all  the  survivors  to  determine  the  distributed 
module  "cut." 

The  sum  of  all  plants  that  failed  due  to  power  module  failures,  com¬ 
mon  module  failures,  and  failures  of  more  than  half  the  lines  in  the  distributed 
module  yields  the  number  of  plants  shut  down  due  to  carbon  fiber  impact.  The 
computer  program  then  turn  to  plants  of  the  next  SIC  number.  These  methods 
constitute  a  calculation  that  generates  a  considerable  increase  in  the  variance 
relative  to  the  Phase  I  methodology. 

Cost  Impact  of  Business/Industry  Closings 

The  Monte  Carlo  submodel  described  above  yields  the  number  of  plants 
in  each  SIC,  size  group  that  are  shut  down  as  the  result  of  each  simulated 
accident.  In  order  to  compute  the  impact  of  those  plant  closings  in  dollar 
terms  we  first  estimate  the  fraction  of  the  industry  shutdown  at  the  location 
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FIGURE  7.3  CONCLUDED 
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that  is  affected.  This  is  done  by  using  an  employee  weighted  fraction  of  pro¬ 
duction  lost.  The  expression  used  in  the  computer  model  is  equivalent  to: 


(F.C.) 


SIC 


(Employees)  (No.  of  Plants  Shut) 

Size _ SIC,  Size _ SIC, Size  (7.9) 

^  >  (Employees)  (No.  of  Plants) 

Size  SIC, Size  SIC, Size 


Equation  (7.9)  provides  an  estimate  of  the  loss  in  capacity  or  output  in  an 
industry  identified  by  one  SIC  number.  The  numerator  is  the  number  of  employees 
in  those  facilities  that  are  shut  down  in  one  SIC  catetory;  the  denominator  is 
the  total  number  of  employees  in  the  same  SIC  category  at  the  same  location. 


The  ORI  risk  assessment  model  estimates  the  impact  of  plant  closings 
due  to  carbon  fiber-related  equipment  failures  by  using  the  Gross  Domestic 
Product  allocated  to  a  particular  business-industry  segment.  The  Gross  Domestic 
Product  (GDP)  is  equal  to  the  Gross  National  Product  (GNP)  debited  by  the  value 
of  foreign  production  of  American  companies  and  credited  with  the  value  of  pro¬ 
duction  by  foreign  companies  in  the  United  States.  In  this  sense  it  measures 
the  value  of  goods  and  services  associated  with  particular  business  and  in¬ 
dustrial  sectors  in  the  United  States.  The  GDP  measures  more  than  value  of 
production  alone  and  is  therefore  the  most  useful  readily  available  economic 
indicator  for  our  use.  GDP  estimates  are  published  at  the  3-digit  SIC  code 
level  on  a  national  basis  by  the  Department  of  Commerce.  In  order  to  allocate 
the  GDP  to  the  local  level  we  used  county-based  payroll  data,  published  in  the 
County  Business  Patterns  for  individual  SIC  numbers.  The  ORI  model  tacitly 
assumes  that  local  productivity  is  essentially  equal  to  the  national  average 
productivity,  on  an  industry-by-industry  basis.  The  GDP  allocable  to  one 
industry  in  a  particular  county  is  estimated  by: 


(County  Payroll  )2jj. 
(National  Payroll) 

SIC 


(GDP) 

SIC 


(7.10) 


The  economic  impact  (not  the  GDP  lost)  is  estimated  by  the  product  of  the 
expression  (7.10)  and  defined  by  Equation  (7.9).  One  further  adjust¬ 

ment  is  required.  Since  the  GDP  data  are  usually  annualized  and  the  payroll 
data  used  for  a  common  time  interval  it  is  necessary  to  multiply  a  factor 
that  is  the  ratio  of  the  length  of  time  (say  in  number  of  days)  a  plant  or 


business  facility  is  shut  down  to  the  number  of  business  days  in  a  year  (de¬ 
fined  as  K).  We  generally  assume  that  a  closing  of  the  type  contemplated  here 
would  last  one  day.  We  therefore  have 


Cost  =  K 


(Local  Payroll 
(National  Payroll) 

SIC 


(GDP)  (F.C.) 

SIC  SIC 


(7.11) 


As  indicated  above,  county  data  appearing  in  County  Business  Patterns  were  used 
for  the  employee  and  payroll  information.  National  data  from  the  Department 
of  Commerce  provided  the  required  GDP  inputs. 

HOUSEHOLD  IMPACT 


The  method  used  in  Phase  II  is  essentially  the  same  as  that  employed 
in  Phase  I.  We  estimated  the  fraction  of  households  in  an  area  that  are  air 
conditioned  (FAC)  and  use  the  methods  previously  described  to  estimate  the 
failure  probability  of  vulnerable  equipment  in  air  conditioned  and  non-air 
conditioned  households.  The  latter  calculation  includes  both  the  failure  and 
ventilation  parameters.  Transfer  functions  for  households  appeared  in  Table 
5.2.  If  the  fraction  of  time  that  a  unit  is  operating  is  T,  then  the  number 
of  failures  of  one  type  of  equipment  is 

HH  X  (No  Equip/HH)  x  T  x  FAC  x  Pp- 

+  HH  X  (No  Equip/HH)  x  T  x  (1-FA)  x  Pp^ 

where  HH  is  the  number  of  households  and  Pp  Pp  are  the  failure 

probabilities  for  the  equipment  in  air  conditioned  and  non-air  conditioned 
households  respectively.  If  the  repair  cost  for  this  equipment  is  RC  dollars, 
then  the  total  estimated  cost  to  repair  all  damaged  equipments  of  a  particular 
class  at  all  households  at  a  location  characterized  by  a  single  exterior  ex- 
plosure  value,  is  given  by: 

RC  X  HH  X  (No  Equip/HH)  x  T  x 

(Pp,  AC  ™  "  Pf,  NAC  (7-13) 

The  locations  and  numbers  of  residential  units  were  obtained  from  the  Bureau 
of  Census  publication.  County  and  City  Data  Book.  Based  on  the  latest  ex¬ 
perimental  evidence  our  attention  was  limited  to  household  television  and 
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high  fidelity  equipment.  Failure  parameters  (E)  for  both  were  set  equal  to 
10®  fiber  seconds  per  cubic  meter;  the  repair  costs  were  estimated  to  be  $50 
per  television  set  and  $100  per  high  fidelity  set.  Updated  Phase  II  ventila¬ 
tion  data  were  incorporated  in  the  calculation  (cf  Table  5.2).  It  was  further 
assumed  that  each  of  these  equipments  would  be  operated  about  half  of  the  time. 
The  equipment  failure  parameters  and  repair  costs  may  be  considered  typical  of 
a  wide  range  of  household  appliances,  so  that  any  two  appliances  may  be  con¬ 
sidered  treated,  rather  than  the  television  and  high  fidelity  sets. 

AIRCRAFT  VULNERABILITY 

Problem  Definition 

In  Phase  I  it  was  concluded  that  key  airport  operations  were  essenti¬ 
ally  invulnerable  to  carbon  fiber  incidents  due  to  the  many  designed  redun- 
diancies  in  the  system.  The  Phase  I  analysis  did  not,  however,  cover  the  risk 
to  aircraft  on  the  ground  at  the  time  of  the  accident.  Because  of  safety-of- 
flight,  as  well  as  other  factors,  it  was  decided  that  an  investigation  should 
be  made  of  the  risk  to  aircraft  on  the  ground,  at  passenger  gates  and  main¬ 
tenance  locations.  This  was  initiated  in  Phase  II,  and  focussed  on  failures 
of  avionics  equipment. 

In  a  cooperative  effort  the  aircraft  manufacturers  analyzed  data 
to  determine  the  number  of  aircraft  expected  to  be  at  passenger  boarding  gates 
and  at  maintenance  locations  on  the  airport  by  day  and  night.  This  was  done 
for  the  nine  airports  previously  selected  to  represent  the  entire  United  States 
(accounting  for  about  one  third  of  U.S.  operations  with  a  bias  toward  the 
larger  airports).  The  results  of  the  airframer  data  collection  effort,  based 
on  current  operations,  were  extrapolated  to  the  1993  time  frame  for  the  ORI 
risk  assessment.  The  results  are  shown  in  Table  7.2. 

For  the  principal  aircraft  types  the  airframers  reviewed  all  onboard 
electrical  and  electronic  equipment.  For  the  L-1011  Tristar,  for  example,  600 
types  of  equipment  were  surveyed,  and  258  components  and  assemblies  were  iden¬ 
tified  for  detailed  vulnerability  review.  After  examination  of  all  pertinent 
characteristics,  84  types  of  equipment  were  identified  as  susceptible  to  CF- 
induced  damage.  All  of  these  types  of  equipment  were  assigned  failure  para¬ 
meters  based  on  available  experimental  data,  extrapolated  where  necessary.  A 
few  examples  from  the  L-1011  are  shown  in  Table  7.3.  The  table  reveals  another 
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TABLE  7.2 


ESTIMATED  NUMBER  OF  AIRCRAFT  EXPOSED  ON  THE  GROUND  -  1993 


AIRPORT 

A/C 

SIZE 

DAY 

NIGHT 

GATE 

MAINT 

GATE 

MAINT 

O' Hare/ 

SMALL 

12 

3 

5 

6 

Chicago 

MED 

36 

0 

18 

11 

LARGE 

39 

8 

18 

20 

Kennedy/ 

SMALL 

12 

2 

3 

5 

New  York 

MED 

49 

0 

21 

65 

LARGE 

41 

8 

23 

32 

Lambert/ 

SMALL 

5 

1 

2 

1 

St.  Louis 

MED 

8 

0 

1 

0 

LARGE 

2 

0 

2 

0 

La  Guardi  a/ 

5 

0 

4 

6 

New  York 

0 

0 

0 

13 

3 

0 

4 

0 

Logun/ 

SMALL 

4 

1 

5 

4 

Boston 

MED 

12 

0 

4 

0 

LARGE 

9 

0 

5 

4 

Phi  la. 

SMALL 

2 

0 

3 

0 

Inf  1 ./ 

MED 

6 

0 

1 

0 

Philadelphia 

LARGE 

6 

0 

8 

0 

Washington 

SMALL 

2 

0 

3 

0 

National/ 

MED 

6 

0 

1 

0 

Washington,  D.C. 

LARGE 

6 

0 

8 

0 

Hartsfield/ 

SMALL 

9 

3 

8 

5 

Atlanta 

MED 

32 

11 

23 

13 

LARGE 

16 

8 

16 

12 

TABLE  7.2 

ESTIMATED  NUMBER  OF  AIRCRAFT  EXPOSED  ON  THE  GROUND  -  1993  (Continued) 


AIRPORT 

A/C 

SIZE 

DAY 

NIGHT 

GATE 

MAINT 

GATE 

MAINT 

Miami 

SMALL 

6 

1 

9 

2 

International/ 

MED 

34 

0. 

43 

26 

Miami 

LARGE 

24 

12 

20 

24 
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TABLE  7.3 

EXAMPLES  OF  VULNERABLE  EQUIPMENT  ABOARD  L-1011  TRISTAR 


EQUIPMENT 
I.D.  NO. 

USE 

NO.  PER 
AIRCRAFT 

LOCATION  ON 
AIRCRAFT 

Failure 

Parameter  ^ 

(E:  Fiber  sec/m'^) 

LI  2 

Radio  Communication 

2 

Avionics  Center 

10^ 

LI  3 

Radio  Communication 

3 

Avionics  Center 

10® 

L32 

Electric  Power 

1 

Flight  Station 

10® 

L65 

Navigation 

2 

Avionics  Center 

1.5  X  10^ 

L69 

Navigation 

1 

Flight  Station 

10® 

L78 

Airborne  Auxiliary 
Power 

1 

Passenger  Cabin 

10® 
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important  factor  in  the  analysis,  to  be  discussed  later:  similar  equipments 
may  be  installed  in  different  parts  of  the  aircraft.  There  are  several  diffe¬ 
rent  possible  paths  for  fiber-laden  air  from  the  exterior  to  the  different 
onboard  equipment  locations. 

In  order  to  consolidate  the  equipment  data  the  airframers  defined  an 
avionics  "suit"  for  a  typical  aircraft  in  each  size  category.  Equipment  type 
classifications  were  made  on  the  basis  of  failure  parameter  and  repair  cost 
primarily,  so  that  it  was  possible  to  reduce  the  total  number  of  types  of 
equipment  considerably  by  appropriate  aggregation.  In  this  way  generic  types 
of  avionics  equipment  were  identified  with  onboard  locations  indicated  for  each 
one,  as  well  as  mean-exposure-to-failure  values,  and  repair  costs.  Table  7.4 
summarizes  the  equipment  input  data  prepared  by  the  airframers:  the  number  of 
the  aircraft,  the  failure  parameter  (T),  and  the  repair  costs. 

As  indicated  above,  avionics  equipment  is  operated  in  several  loca¬ 
tions  on  the  aircraft;  this  factor  together  with  the  possibility  of  various 
doors  and  hatches  being  open  or  shut  resulted  in  the  definition  of  different 
ventilation  modes  for  each  equipment-aircraft  combination.  Here  again,  an 
independent  analysis  by  the  airframe  manufacturers  provided  values  of  the 
different  transfer  functions  and  the  fraction  of  time  (during  day  and  night) 
that  each  would  be  expected  to  prevail.  These  results  are  summarized  in 
Table  7.5. 

Computer  Methods 

Figure  7.4  is  a  flow  chart  illustrating  the  computer  submodel  that 
computes  avionics  failures  and  resulting  costs.  The  first  step  in  the  com¬ 
putation  is  the  determination,  on  a  stochastic  basis,  of  whether  the  simulated 
accident  took  place  during  the  day  or  night.  The  conditional  probability  for 
this  event  is  based  on  the  analysis  of  the  airframer  aircraft  accident  data 
base.  The  calculation  proceeds  through  each  aircraft  size  in  turn.  For  each 
size  aircraft  the  program  "looks  up"  the  number  of  aircraft  at  the  predefined 
gate  and  maintenance  locations  (for  the  airport  being  simulated).  At  each  loca¬ 
tion  in-turn  the  number  of  each  type  of  equipment  in  aggregated;  on  a  random 
basis  the  model  determines  the  number  in  each  of  the  predefined  ventilation 
modes  (i.e.,  finds  the  applicable  transfer  function).  With  the  value  of  the 
transfer  function  and  the  exterior  exposure  for  the  particular  location  the 
model  computes  the  failure  probability  for  the  equipment,  using  the  input 
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TABLE  7.4 

AIRCRAFT  AVIONICS  EQUIPMENT  CONFIGURATIONS  WITH  FAILURE  AND  COST  INPUTS 


Aircraft 

Size 

Avionics 

Equipment 

ORI  I.D.  No. 

Number  on 

Ai rcraf t 

E 

(Failure 

Parameter*) 

Repair  Cost 
(  $  ) 

1 

38 

10® 

100 

2 

7 

1.5  X  10^ 

100 

Small 

3 

6 

10® 

450 

4 

2 

1.5  X  10^ 

450 

5 

1 

10® 

300 

6 

18 

10® 

50 

7 

26 

1.5  X  10^ 

215 

8 

24 

10® 

220 

Medium 

9 

153 

10® 

175 

& 

10 

4 

10® 

250 

Large 

11 

22 

10® 

210 

12 

43 

10® 

385 

13 

3 

10® 

530 

14 

2 

10® 

1295 

15 

4 

10® 

1665 

*In  fiber-seconds  per  cubic  meter. 
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TABLE  7.5 

VENTILATION  FACTORS  (T.F.)  AND  ASSOCIATED  PROBABILITIES  FOR  AVIONICS 
EQUIPMENT  ABOARD  PARKED  AIRCRAFT  BY  LOCATION  AND  TIME  OF  DAY 


Gate 

Maintenance 

Equip. 

No. 

Day 

Night 

Day  &  Night 

Prob. 

T.F. 

Prob. 

T.F. 

Prob. 

T.F. 

1  - 

.99 

.70 

.70 

.70 

.23 

.70 

4 

.01 

1.0 

.30 

1.0 

.77 

1.0 

5, 

6 

.99 

.70 

1.0 

.70 

.96 

.70 

.01 

.0025 

.04 

.0025 

7 

.95 

.0025 

.20 

.0025 

.79 

1.0 

.01 

1.0 

.50 

.01 

.14 

.01 

.04 

.01 

.30 

1.0 

.07 

.0025 

8 

1.0 

.01 

1.0 

.01 

1.0 

.01 

9. 

10 

.99 

.01 

1.0 

.01 

.96 

.01 

.01 

.0025 

.04 

.0025 

11  - 

.95 

.0025 

.20 

.0025 

.78 

1.0 

15 

.01 

1.0 

.50 

.01 

.14 

.01 

.04 

.01 

.30 

1.0 

.08 

.0025 
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SET  GATE  LOCATION 


FIGURE  7.4  FLOWCHART  FOR  COMPUTING  AVIONICS  FAILURES  AND  ASSOCIATED  COSTS 
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value  of  E.  Using  a  Poisson  distribution  the  actual  number  of  failures  is  ob¬ 
tained  in  a  random  draw.  Using  the  input  repair  cost  for  the  equipment  the 
computer  model  then  determines  the  total  repair  cost.  This  procedure  is  re¬ 
peated  for  each  type  of  equipment  on  one  size  aircraft,  then  done  for  the 
other  size  aircraft,  then  repeated  at  the  next  location. 

COSTING  SUMMARY 

The  input  requirements  for  the  business/industry  impact  cost  model 
are  summarized  in  Table  7.6.  All  counties  within  50  miles  of  the  airport  are 
defined  by  a  set  of  geographical  coordinates.  At  one  geographical  location 
the  model  computes  business-industry  impact  as  the  sum  of  costs  of  equipment 
repair,  facility  cleanup,  and  business  disruption.  At  those  locations  defined 
as  residential  centers  the  model  computes  the  total  cost  due  to  household  equip¬ 
ment  failures.  At  the  airport  itself  the  model  computes  costs  required  to 
repair  failed  avionics  equipment.  Summary  results  for  each  simulated  accident 
present  the  total  of  costs  in  each  of  these  three  major  categories,  obtained 
by  adding  the  costs  over  all  geographical  locations  affected  by  the  simulated 
accident. 
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TABLE  7.6 
SUMMARY 

BUSINESS/ INDUSTRY  IMPACT  COST  MODEL 
INPUT  REQUIREMENTS 


Level 

Descriptor 

Input  Definition 

National 

SIC  Number 

Payroll 

Gross  Domestic  Product 

County 

SIC  Number 

Local  Payroll 

No.  of  Establishments  by  Size 

Facility 

SIC  Number 

&  Size 

No.  of  Equipments  by  Type 

Plant  Configuration 

Cleanup  Cost 

Equipment 

Standard  Type 

Repair  Cost 

Mean  Exposure-to-Failure 
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VIII.  INDIVIDUAL  AIRPORT  RESULTS 


The  simulation  model  was  run  for  a  large  sample  of  accidents  at  the 
following  nine  airports  previously  examined  in  Phase  I: 

O' Hare/Chicago 

John  F.  Kennedy/New  York  City 

Washington  National  Airport/Washington,  D.C. 

Lambert/St.  Louis 

LaGuardia/New  York  City 

Logan/Boston 

Hartsf i el d/Atl anta 

Miami  International/Miami 

Philadelphia  International /Philadelphia . 

SAMPLE  ACCIDENT 

To  set  the  stage  for  interpreting  the  airport  simulation  statistics 
we  present  detailed  output  for  one  random  accident  generated  at  Kennedy  Air¬ 
port.  The  basic  geometry  is  shown  in  Figure  8.1;  the  airport  and  accident 
location  are  indicated.  Randomly  generated  weather  data  are:  wind  from  the 
south  at  2  meters  per  second  and  stability  class  6.  The  accident,  based  on 
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Figure  8.1  Basic  Geography  Associated  With  One  Random  Accident  At 
Kennedy  Airport  in  New  York.  Airplane  shows  airport 
location;  asterisk  indicates  accident  location.  The 
circle  represents  Queens  County  for  modelling  purposes. 
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randomly-selected  factors  involved  a  large  jet  (type  16)  with  3000  kilograms 
of  composite  aboard;  1.3  x  10^^  fibers  were  released  in  a  take-off  accident, 
with  the  plume  formed  during  the  fire  limited  to  a  height  of  100  meters  by 

4 

an  inversion.  As  a  result  the  mean  exposure  in  Queens  County  is  3.5  x  10 
fiber-seconds  per  cubic  meter;  no  other  neighboring  counties  were  significantly 
affected. 

The  results  for  this  one  accident  include  the  computed  cost  of  re¬ 
pairs  to  household  equipment  in  Queens  of  $533.  The  impact  on  business  com¬ 
prises  $5,600  in  equipment  repair  costs  and  $31,577  due  to  business  closings 
and  cleanup  as  the  result  of  equipment  failure. 

SIMULATION  OUTPUTS 

A  typical  set  of  runs  includes  approximately  2500  simulated  accidents, 
each  generating  data  of  the  type  presented  above  for  one  accident.  The  com¬ 
puter  program  summarizes  the  data  from  all  simulated  accidents  to  provide  the 
following  outputs: 

•  Characteristics  of  the  ten  most  costly  accidents 

•  Probability  distribution  of  annual  costs  for  household,  in¬ 
dustrial,  and  avionics,  as  well  as  the  mean,  standard  devia¬ 
tion,  and  risk  profile 

•  Probability  distribution  of  costs  per  accident  for  household, 
industrial,  and  avionics  impacts  to  provide:  mean,  standard 
deviation,  and  risk  profile 

•  Distribution  of  number  of  accidents  per  year  (replication). 

Available  options  permit  printing  out  the  details  associated  with  each  accident, 
as  presented  above.  A  standard  printout  for  a  sample  (annual  replications)  re¬ 
sults  for  34,000  replications  at  Kennedy  Airport  appears  in  Figure  8.2.  The 
results  indicate  that,  for  example,  in  7  runs  the  total  cost  was  greater  than 
$100,  and  less  than  $178,  as  shown  in  the  column  headed  "TOTAL  COUNTS."  The 
class  intervals  were  selected  to  be  equal  on  a  logarithmic  scale  to  facilitate 
computing  the  risk  profiles,  essentially  cumulative  probability  distributions. 
The  risk  profiles  for  annual  household  avionics  failures,  and  business-industry 
impact  plotted  from  these  output  data  appear  in  Figure  8.3. 
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FIGURE  8.2  Annual  Result  Output,  for  34,000  Annual  Replications  -  Kennedy  Airport  -  1993 
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Figure  8.3  Annual  Risk  Profiles  for  Kennedy  Airport  -  1993:  Business/ 

Industry,  Household  Equipment,  and  Avionics  Equipment  Impacts; 
and  Total  Impact. 
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COMPARISON  OF  DIFFERENT  AIRPORTS 


In  this  section  representative  results  are  summarized  for  the  nine 
airports  listed  above.  For  each  airport  the  number  of  samples  (replications) 
was  selected  so  that  at  least  2500  accidents  were  simulated.  Computer  time 
is  directly  proportional  to  the  number  of  accidents  being  simulated,  and  com¬ 
parability  of  results  for  different  airports  required  a  common  basis.  Re¬ 
sults  presented  below  show  the  impact  of  changing  this  number. 

Table  8.1  summarizes  the  results  of  all  the  single  airport  accident 
simulations.  The  table  presents  accident  data  as  contrasted  with  annual  data. 
The  probability  of  more  than  one  accident  per  year  is  so  small,  however,  that 
the  average  of  results  over  all  accidents  is  essentially  equal  to  the  average 
of  results  over  all  simulated  years.  In  each  set  of  runs  the  program  presents 
detailed  information  about  the  ten  worst  (highest  cost)  accidents.  The  aver¬ 
age  of  the  ten  highest-cost  accidents  at  each  of  the  airports  also  appears  in 
the  table.  At  each  airport  these  ten  accidents  (0.4  percent  of  all  simulated 
accidents)  comprise  the  highest  cost. 

The  results  in  Table  8.1  show  that,  typically,  the  costs  resulting 
from  business  and  industrial  impact  are  considerably  greater  than  the  house¬ 
hold  impact  costs  and  the  avionics  failure  costs.  All  mean  costs  appear  rela¬ 
tively  small.  To  present  some  idea  of  the  range  of  these  results  we  note  here 
that  the  maximum  cost  in  each  of  the  three  categories  and  the  airport  at  which 
it  occurred  are: 

•  Households:  $2,665  at  Kennedy  Airport,  New  York 

t  Business-Industry:  $274,000  at  Logan  Airport,  Boston 

•  Avionics:  $3,910  at  Kennedy  Airport,  New  York. 

As  indicated,  these  represent  results  from  2500  accident  simulations  at  each 
airport,  so  that  these  extreme  values  were  experienced  with  an  empirical 
frequency  of  4  in  ten  thousand.  The  extreme  values  quoted  here  are  actually 
the  maxima  from  a  sample  of  approximately  9x2500  or  22,500  accidents.  The 
likelihood  of  an  accident  with  fire  in  any  year  Is  quite  low;  the  extreme  values 
reported  for  Kennedy  airport  and  Logan  airport  would  have  occurred  only  once 
in  34,000  years  and  once  in  67,000  years,  respectively. 
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993  RESULTS  FOR  SIMULATED  ACCIDENTS  AT  NINE  AIRPORTS  (1976  DOLLARS) 
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The  computer  simulation  model  also  generates  all  the  results  needed 
to  plot  the  risk  profiles,  as  shown  previously  in  the  sample  printout  appearing 
as  Figure  8.2.  The  risk  profile  is  typically  presented  on  an  annual  basis. 
Figure  8.4  shows  the  annual  risk  profiles  for  O'Hare  Airport,  Chicago,  Lambert 
Airport,  St.  Louis,  and  Hartsfield  Airport,  Atlanta.  O'Hare/Chicago,  the 
Nation's  busiest  airport,  has  a  risk  profile  that  shows  that  the  probability  of 
exceeding  $10,000  per  year  in  total  CF-related  impact  is  approximately  .0004. 
For  St.  Louis,  the  corresponding  probability  is  approximately  .0001  (one  in 
10,000).  These  three  airports  constitute  a  sample  of  different  combinations 
of  annual  commercial  operations  and  surrounding  population,  as  summarized 
here: 

O' Hare  (ORD):  High  population,  heavy  air  traffic 

hartsfield  (ATL) :  Low  population;  high  traffic  level 

Lambert  (STL):  Low  population,  low  traffic  level. 

This  stratification  is  reflected  in  the  annual  risk  profiles.  The  O'Hare  risk 
is  highest,  St.  Louis  lowest. 

The  Phase  II  computer  program  was  modified  to  generate  statistics 
on  a  per-accident  basis  as  well  as  the  customary  per-year  basis.  Figure  8.5 
shows  the  accident  risk  profiles  for  the  same  three  airports.  The  risk  is 
greatest  for  O’Hare,  due  to  the  relatively  high  concentration  of  business 
and  industry;  the  St.  Louis  and  Atlanta  risk  profiles  are  quite  similar,  in¬ 
dicating  that  their  separation  in  Figure  8.4  was  due  to  the  difference  in 
accident  incidence  (i.e.  in  our  model,  the  difference  in  number  of  operations). 
The  accident  risk  profiles  may  be  considered  conditional  probabilities.  For 
example,  given  that  a  CF-built  aircraft  crashes  and  burns  at  O'Hare  Airport 
the  probability  is  1  in  a  hundred  (.01)  that  the  impact  will  exceed  one  thous¬ 
and  dollars  ($1,000);  for  $100,000  the  probability  is  2  in  one  hundred  thous¬ 
and  (.00002),  by  extrapolation  from  Figure  8.5. 

STATISTICAL  CONFIDENCE  LIMITS 

In  Phase  I  ORI  shows  that  the  simulation  runs  for  one  airport  may  be 
considered  a  set  of  Bernoulli  trials.  As  a  result  we  derived  the  following 
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Figure  8.5  Accident  Risk  Profiles  -  1993  Scenario  -  for  O' Hare/Chicago  (ORD), 
Hartsfield/Atlanta  (ATL)  and  Lambert/St.  Louis  (STL). 
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expression  for  the  95%  statistical  confidence  limits: 

P  t  2  ^p(l-p.).. 

where  p  is  the  computed  exceedance  probability  after  simulating  n  samples. 

Figure  8.6  shows  the  Washington  National  Airport  risk  profile  with  the  95% 
confidence  limits.  The  confidence  limits  apply  to  the  purely  statistical 
nature  of  the  simulation,  and  not  to  the  impact  of  errors  in  input  data.  The 
results  do  show  that  conclusions  based  on  the  risk  profiles  need  not  be 
altered  because  of  inherent  statistical  uncertainty.  The  confidence  limit 
bars  shown  on  the  graph  appear  to  not  be  of  equal  size  above  and  below  the 
curve  due  to  fact  that  the  results  are  plotted  on  a  logarithumic  scale. 

ADEQUACY  OF  SAMPLE  SIZE 

In  a  major  simulation  modelling  effort  of  the  type  reported  here  one 
of  the  important  questions  is  whether  enough  runs  have  been  made.  This  is 
related  to  the  stability  of  the  model  and  variance  in  the  input  data.  Rather 
than  invoke  sophisticated  statistical  arguments  it  is  more  convenient  to  let 
the  results  "speak  for  themselves."  In  effect  we  compared  the  results  for 
two  different  numbers  of  simulation  runs. 

The  O'Hare  Airport/Chicago  simulation  was  run  for  22,000  and  44,000 
annual  samples,  resulting  in  2537  and  5038  accidents  respectively.  It  is  not 
possible  to  compare  the  two  risk  profiles  on  a  graph  using  a  scale  convenient 
for  this  report,  since  they  would  be  too  close  to  one  another.  We  are  there¬ 
fore  limited  to  the  results  summarized  in  Table  8.2.  A  significantly  larger- 
cost  accident  occurred  in  the  44,000-sample  run  than  in  the  22,000-sample  run 
which  is  typical  of  extreme-value  statistics.  In  this  case  the  contribution 
of  the  larger  accident  results  in  the  mean  values  being  somewhat  different. 

The  risk  probabilities  are,  however,  quite  similar.  It  is  also  interesting  to 
note  that  five  of  the  ten  highest-cost  accidents  in  the  44,000-sample  run  oc¬ 
curred  in  "second  half",  i.e.,  in  samples  after  number  22,000. 

SENSITIVITY  TESTS 

To  demonstrate  the  flexibility  of  the  ORI  Carbon  Fiber  Risk  Assess¬ 
ment  Model,  as  well  as  to  provide  insight  into  the  physical  mechanisms  at  work, 
several  input  parameters  were  varied  and  the  impact  of  the  variation  on  the  out¬ 
put  examined.  These  results  are  described  in  this  section  of  the  report. 


EXCEEDANCE  PROBABILITY 


TABLE  8.2 

1993  CHICAGO/O'HARE  COMPARISON  OF  DIFFERENT  SETS  OF  SIMULATIONS 


Measure 

22,000  Samples 

44,000  Samples 

No.  of  Accidents 

2537 

5038 

Mean  Accident 

$147 

$166 

Worst  Accident 

$54,000 

$110,299 

P  (Annual  Cost  >$1000) 

.000955 

.00111 

P  (Annual  Cost>$l 0,000) 

.000545 

.000545 

Inversion  "Punch  Through" 


It  is  assumed  in  the  results  reported  so  far  that  the  plume  does  not 
penetrate  the  inversion.  For  one  set  of  Washington  National  Airport  simulations 
we  permitted  the  plume  to  "punch  through"  the  inversion.  This  was  a  rela¬ 
tively  simple  program  modification  to  introduce,  since  the  standard  model  in¬ 
cludes  a  test  comparing  the  computed  plume  height  with  the  inversion  height. 

If  the  initially-computed  plume  height  is  greater  it  is  reduced  and  set  equal 
to  the  inversion  height  for  the  remainder  of  the  calculation.  In  the  sensitivity 
test  reported  here  this  comparison  was  bypassed  and  the  computed  plume  height 
was  used,  regardless  of  its  magnitude  relative  to  the  height  of  the  inversion. 

The  results  for  Washington  National  Airport  indicate  risks  so  low 
compared  to  the  base  case  that  no  risk  profiles  were  drawn.  It  is  only  neces¬ 
sary  to  cite  a  few  values  to  make  the  point.  In  the  test  case,  with  punching 
through  permitted,  the  probability  of  exceeding  $100  per  year  in  CF  -  impact 
costs  was  .000045  compared  to  the  1993  Washington  National  Airport  base  case 
(no  punch-through)  value  of  .003,  roughly  a  difference  by  a  factor  of  twenty. 

For  exceeding  $1000  per  year  the  corresponding  probabilities  are  .000015 
(punch  through)  and  .00048  (no  punch  through).  In  the  test  cases  the  mean 
annual  accident  cost  is  reported  as  zero  (actually  less  than  50d)  compared  to 
$12  in  the  standard  case.  Another  way  of  reporting  these  results  is  that,  of 
2590  random  accidents  generated  in  the  "punch-through"  runs,  only  three  had 
associated  CF-related  costs  of  more  than  $100. 

One  interesting  result  relates  to  the  stability  class  associated 
with  the  ten  most  costly  accidents.  In  the  base  case  these  are  all  class  6 

or  class  5,  -  the  most  stable  atmospheric  conditions;  these  stability  classes 

are  characterized  by  a  100-meter  inversion  height.  In  the  standard  simulations 
this  was  also  the  height  at  which  the  plume  was  stopped.  For  the  "punch- 

through"  runs  the  ten  most  costly  accidents  are  associated  with  stability  class 

6  although  the  average  plume  height  was  437  meters;  the  average  cost  incurred 
in  these  ten  costliest  accidents  was  $571,  compared  to  $62,497  (cf.  Table  8.1) 
in  the  base  case. 

Random  Inversion  Height 

In  other  results  presented  in  this  report  the  inversion  height  is 
linked  to  the  stability  class  by  a  one-to-one  relationship  for  each  airport. 
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This  is,  in  each  of  the  airport  accident  simulations  the  random  selection  of 
the  atmospheric  stability  category  automatically  determines  the  height  of  the 
inversion  (cf.  Section  IV).  As  called  for  in  the  Phase  II  contract  ORI  tested 
the  effect  of  this  approach  by  devising  a  methodology  to  provide  for  a  random 
selection  of  the  inversion  height.  This  was  done  by  first  associating  the 
different  stability  classes  with  the  period  of  the  day  during  which  each  is 
most  likely  to  occur.  For  example,  stability  classes  4,  5,  and  6  usually  pre¬ 
vail  during  the  night.  The  previously  developed  values  of  the  inversion 
height  were  selected  as  the  values  prevailing  at  the  midpoints  of  each  of  the 
appropriate  time  intervals  and  were  then  connected  by  straight  line  segments. 

The  actual  inversion  height  was  then  determined  by  a  random  selection  from 
that  part  of  the  continuous  inversion  height-time  relationship  appropriate  to 
the  randomly-selected  stability  class. 

The  resulting  risk  profiles  are  compared  in  Figure  8.7  for  the  1993 
Washington  National  Airport  scenario.  The  standard,  fixed  inversion  height 
per  stability  class  case  is  characterized  by  a  somewhat  higher  risk,  showing 
that  the  method  previously  used  is  relatively  conservative.  The  annual  mean 
impact  is  $12  for  the  fixed  inversion  height  case  and  drops  to  $5  in  the  vari¬ 
able  inversion  height  case.  The  corresponding  average  impacts  of  the  ten 
worst  accidents  are  $62,497  and  $28,994,  respectively.  It  also  turned  out 
that,  in  the  variable  inversion  height  runs,  the  ten  worst  accidents  were 
associated  with  stability  classes  5  and  6,  the  most  stable,  although  the  in¬ 
version  height  was  not  always  set  at  100  meters  as  it  is  in  the  standard  runs 
for  these  stability  classes. 

0*Hare  Airport  “Worst  Case" 

In  order  to  examine  the  impact  of  a  drastic  change  in  the  underlying 
assumptions  the  O'Hare  Airport  1993  scenario  simulations  were  run  with  the  in¬ 
puts  changed  so  that  all  aircraft  operations  involved  aircraft  "loaded"  with 
CF.  In  order  to  do  this  conveniently  the  inputs  were  adjusted  so  that  all  air¬ 
craft  operating  at  O'Hare  in  1993  with  CF  in  their  structure  were  our  previously 
defined  type  number  19  (cf.  Table  2.2).  This  is  a  large  jet  with  the  most  CF 
composite  of  any  plane  expected  be  in  the  1993  fleet;  the  total  onboard  is 
15,600  kilograms.  The  average  amount  of  composite  aboard  all  1993  aircraft  is 
only  about  2,800  kilograms  per  aircraft  with  composite  while  the  range  is  from 
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Figure  8.7  1993  Washington  National  Airport  Annual  Risk  Profiles  Showing 

Effect  of  Random  Selection  of  Inversion  Height. 
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65  to  15,600,  for  the  aircraft  with  composite.  The  worst-case  scenario  effec¬ 
tively  increases  the  amount  of  carbon  fiber  liberated  in  an  accident  by  a  fac¬ 
tor  of  approximately  five,  on  the  average.  Several  of  the  output  results  from 
the  two  sets  of  simulation  runs  are  compared  in  Table  8.3.  The  interesting 
result  is  that  the  mean  annual  impact  and  the  mean  accident  cost  each  increase 
by  a  factor  of  approximately  ten  due  to  the  average  increase  in  CF-release 
of  a  factor  of  about  five.  In  the  case  of  the  ten  most  costly  accidents  the 
ratio  of  the  mean  impacts  is  approximately  2.5.  The  average  amount  of  com¬ 
posite  released  in  the  ten  most  costly  accidents  for  the  worst-case  scenario 
is  about  5.8  times  the  average  for  the  corresponding  best-estimate  cases. 

The  results  for  the  probabilities  shown  in  Table  8.3  indicate  that 
the  effect  of  the  increase  in  the  average  amount  of  fiber  per  aircraft  is  to 
shift  the  peak  of  the  frequency  distribution  of  accident  impact  to  the  right; 
the  most  trival  accidents  have  been  eliminated  and  this  causes  the  larger 
shift  in  the  mean  values.  In  the  standard  runs  the  standard  deviation  of  annual 
impact  is  equal  to  about  40  times  the  annual  mean  impact;  in  the  worst  case  this 
ratio  is  about  15.  The  best-estimate  1993  O'Hare  annual  risk  profile  is  com¬ 
pared  with  worst  case  in  Figure  8.8.  The  comparison  shows  the  significant 
impact  of  the  increased  carbon  fiber.  However,  the  probability  of  exceeding 
$10,000  in  annual  damages  is  only  about  .005  (five  in  a  thousand),  which  is 
also  indicated  in  Table  8.3. 
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TABLE  8.3 

COMPARISON  OF  SELECTED  OUTPUTS:  O' HARE  AIRPORT.  1993,  BEST-ESTIMATE 
AIRCRAFT  MIX  AMD  ALL  AIRCRAFT  "LOADED" 
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EXCEEDANCE  PROBABILITY 


Figure  8.8  1993  Chicago/0 'Hare  Annual  Risk  Profile  Compared  to  "Worst 
Case." 
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IX.  NATIONAL  RISK 


In  Phase  I  ORI  used  a  mathematical  technique  called  a  convolution 
to  generate  the  national  risk  profile.  This  method  used  as  its  input  the  final 
results  of  the  airport  simulations:  the  risk  profiles,  or  probability  dis¬ 
tributions.  In  Phase  II  a  somewhat  more  straightforward  approach  was  used, 
which  takes  as  its  input  the  individual  accident  results  generated  for  each 
of  the  airports.  The  method  is  described  in  more  detail  below,  followed  by 
a  presentation  of  the  results. 

METHOD 

The  airports  previously  treated,  and  listed  in  Section  VIII,  account 
for  approximately  one-third  of  the  Nation's  commercial  air  traffic.  To  compute 
the  national  risk  it  is  assumed  that  these  airports  can  be  used  with  suitable 
adjustment  to  represent  the  entire  United  States,  at  least  as  far  as  commer¬ 
cial  aircraft-related  CF  risk  is  concerned.  Since  these  are  predominantly 
large,  busy  airports  near  large  metropolitan  areas,  this  method  may  be 
expected  to  overestimate  the  national  risk;  for  our  purposes  this  constitutes 
a  conservative,  and  therefore  desirable,  result. 

We  have  already  discussed  our  method  of  estimating  the  average 
number  of  accidents  in  the  U.S.  in  one  year  involving  commercial  aircraft 
with  CF  aboard  in  a  fire  (cf  section  II).  In  conducting  a  national  simulation 
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we  assume  that  the  actual  number  of  accidents  in  a  year  obeys  a  Poisson  distri¬ 
bution  with  mean  equal  to  2.6.  The  ORI  national  risk  model  then  draws  a  number 
of  accidents  at  random  from  this  distribution.  The  method  is  illustrated  in 
Figure  9.1.  Each  of  the  accidents  in  a  year  (one  replication)  is  assigned  to 
one  of  the  airports  previously  examined.  The  probability  that  the  accident 
takes  place  at  one  airport  is  simply  the  ratio  of  the  number  of  operations  at 
that  airport  to  the  total  for  all  nine  airports.  This  again  biases  the  results 
in  favor  of  the  relatively  busy  airports  we  examined  previously,  and  this  errs 
in  the  conservative  direction.  The  important  concept  in  this  method  is  that 
the  details  of  the  individual  aircraft  accidents  previously  simulated  at 
all  nine  airports  are  saved  and  used  here.  Figure  9.2  is  a  conceptual  repre¬ 
sentation  of  the  file  for  one  airport's  simulated  accident.  In  the  national 
calculation  we  simply  draw  one  of  these  accidents  at  random  off  the  list  for 
the  airport  to  which  the  national  model  has  allocated  the  accident.  This  pro¬ 
cess  is  repeated  for  all  accidents  in  the  replication  (i.e.,  simulated  year). 

The  sum  of  costs  over  all  accidents  in  one  year's  sample  is  a  conservative 
estimate  of  the  total  national  impact. 

RESULTS 

The  results  of  the  calculation,  using  outputs  from  the  individual 
airports,  and  the  weighting  factors  described  above,  indicate  a  maximum  annual 
impact  on  business  and  industry  of  $274,000,  with  a  mean  of  $466.  For  avionics 
impact  the  results  are  $3,900  and  $2,  respectively.  The  method  is  entirely 
analogous  to  that  used  for  the  single  airports  in  the  sense  of  being  Bernoulli 
trials,  and,  therefore,  our  previous  method  of  computing  confidence  limits  can 
be  used.  In  this  case  we  typically  simulated  about  10,000  "national  years." 

The  national  risk  profile  with  the  95%  statistical  confidence  limits 
is  shown  in  Figure  9-3.  The  risk  profile  indicates  that  the  probability  of 
exceeding  $10,000  per  year  in  CF-related  impact  is  about  1  in  a  hundred  (.01). 
The  statistical  confidence  limits  indicate  that  the  probability  is  .95  that  the 
actual  probability  in  this  case  is  between  .0133  and  .00867.  The  sensitivity 
tests  reported  previously  for  one  airport  indicate  that,  with  relatively  drastic 
worst-case  inputs  the  results  do  not  change  dramatically,  thus  increasing  our 
overall  confidence  in  the  national  results  presented  here.  In  Figure  9.4  the 
Phase  I  and  Phase  II  results  are  compared,  showing  that  the  new  Phase  II  inputs 
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FIGURE  9.1  FLOWCHART  FOR  NATIONAL  RISK  ASSESSMENT  MODEL 
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FIGURE  9.2  CONCEPTUAL  TABULAR  ARRAY  OF  DETAILED  RESULTS  OF  SINGLE  AIRPORT  SIMULATION  RUNS 


FIGURE  9.3  1993  NATIONAL  RISK  PROFILE  WITH  95%  STATISTICAL  CONFIDENCE  LIMITS 
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FIGURE  9.4  COMPARISON  OF  PHASE  I  AND  PHASE  II  NATIONAL  RISK  PROFILES 
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result  in  a  greatly  reduced  estimated  risk,  despite  the  fact 
model  is  considerably  more  likely  to  generate  results  with  a 
(extreme  values). 


that  the  Phase  II 
large  variance 
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X.  SUMMARY  AND  CONCLUSIONS 


ORI,  Inc.,  has  developed  and  demonstrated  a  general  stochastic  risk 
assessment  model.  This  model  has  been  applied  to  the  assessment  of  the  risk 
associated  with  proposed  applications  of  carbon  fiber  composite  materials  in 
commercial  aircraft  structures.  The  airport  risk  model  replicates  key  elements 
of  the  accident  scenario  (i.e.  fire,  plume,  diffusion  and  transport,  transfer, 
equipment  failure,  and  economic  impact)  in  a  logic  structure  that  supports  the 
calculation  of  statistical  measures  of  the  risk.  The  national  risk  is  designed 
to  use  the  results  from  several  individual  airport-city  calculations. 

One  of  the  principal  changes  from  Phase  I  to  Phase  II  has  been  the 
availability  of  improved  input  data,  particularly  for  key  parameters  associated 
with  the  amount  of  fiber  in  a  projected  aircraft,  the  fraction  of  carbon  fiber 
released  in  a  burn,  and  the  vulnerability  of  electronic  equipment.  Where  sound 
experimental  data  is  not  available  we  have  always  taken  a  conservative  approach; 
that  is,  the  inputs  are  selected  so  as  to  maximize  the  estimated  risk.  Experi¬ 
mental  data  indicates  that  the  exponential  law  tends  to  overestimate  the  prob¬ 
ability  of  equipment  failure  for  low  values  of  carbon  fiber  exposure.  We  have 
used  the  exponential  law,  since  it  is  relatively  simple,  and  is  conservative. 

The  model  has  been  subjected  to  a  variety  of  tests  all  showing  its 
inherent  stability.  We  have  developed  statistical  confidence  bounds  about  our 
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risk  probabilities.  The  model  outputs  have  been  subject  to  a  variety  of  sen¬ 
sitivity  tests;  these  have  shown  that  even  when  drastic  unrealistic  shifts  in 
the  input  data  are  introduced  the  resulting  annual  risk  increases  by  less  than 
a  factor  of  ten. 

The  national  model  was  run  using  nine  relatively  busy  airports  heavi 
ly  weighted  toward  those  with  large  concentrations  of  business  and  industry, 
and  thus  overestimates  the  national  risk.  Even  in  that  case,  the  chance  of  ex 
ceeding  an  annual  economic  impact  of  more  than  $100,000  due  to  aircraft-fire 
accidents  and  subsequent  carbon  fiber  release  is  less  than  one  in  a  thousand. 

A  worst-case  analysis  of  the  possible  impact  of  such  accidents  on  the  electric 
power  distribution  system  showed  that  the  expected  number  of  outages  would  be 
negligible  compared  to  outages  due  to  other  causes. 

Our  conclusion,  based  on  the  results  presented  earlier  in  this  re¬ 
port,  and  summarized  briefly  above  is  that: 

The  risk  due  to  accidental  release  of  carbon  fibers  following  an 
accident  and  fire  involving  civil  aircraft  is  quite  small. 
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THE  EFFECTS  OF  CARBON  FIBER 
EXPOSURE  ON  ELECTRIC  UTILITY  SYSTEMS 


INTRODUCTION 

Objectives 

This  study  assesses  the  expected  impacts  of  the  exposure  of  an 
electric  utility  system  to  carbon  graphite  fibers  in  response  to  the  combined 
requirements  of  Tasks  8E  and  9  of  the  scope  of  work:— 

•  Task  8E  -  "is  directed  toward  the  assessment  of  the  likelihood 
of  failures  of  the  electric  power  distribution  system  due  to 
carbon  fiber  release  accidents.  Information  is  available  that 
relates  failures  of  individual  critical  components  of  the  sys¬ 
tem  to  different  levels  of  carbon  fiber  exposure." 

•  Task  9  -  "In  order  to  properly  assess  the  possible  impact  of  car¬ 
bon  fiber  incidents  on  the  national  power  supply  system,  it  is 
necessary  to  develop  a  historical  perspective  on  previous  break¬ 
downs  of  the  system." 

Scope 

Based  on  a  review  of  literature  and  visits  to  operating  electric 
utilities,  for  Task  8E  we  have  defined  typical  electric  power  distribution 

Certain  measurements  in  this  Appendix  are  in  English  Units.  The  following 
factors  may  be  used  to  convert  these  measurements  to  SI  Units:  (1)  kilometer  = 
0.62  mile;  (2)  square  kilometer  =  0.3861  square  mile. 
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systems  in  the  2.4KV-38KV  range,  using  actual  distribution  circuits  to  the 
maximum  extent  possible.  Using  CF  exposure  failure  threshold  (r)  values  pro¬ 
vided  by  Westinghouse  Electric  Company  for  various  types  of  insulators  and 
bushings,  expected  customer  outages  due  to  CF  are  estimated  for  the  typical 
distribution  circuits. 

Using  published  information  and  visits  to  electric  utilities  and 
to  selected  industrial  facilities,  for  Task  9  we  have  estimated  the  frequency 
and  duration  of  failures  to  power  systems  and  described  the  protective  measures 
taken  by  both  the  utilities  and  by  business  and  industrial  users.  This  infor¬ 
mation  is  used  as  a  baseline  against  which  to  assess  the  expected  outages  due 
to  exposure  to  carbon  fibers. 

Background 

Tests  which  have  been  performed  by  the  U.S.  Army  Ballistics  Research 
Laboratory  and  by  the  Westinghouse  Electric  Company  indicate  that  the  effects 
of  carbon  fibers  on  electrical  systems  operating  in  excess  of  about  38KV  can 
be  neglected.  Likewise,  the  effects  of  carbon  fiber  exposure  on  the  secondary 
side  of  distribution  transformers  serving  industries,  businesses  and  residences 
are  addressed  as  a  part  of  the  general  NASA  assessment  of  carbon  fiber  risks, 
described  in  the  main  body  of  this  report.  There  is  a  need  to  analyze  the 
potential  risks  due  to  CF  exposure  primary  distribution  circuits  which  operate 
in  the  2.4KV  to  38KV  range. 

Tests  have  been  performed  by  Westinghouse  Electric  Company  on  var¬ 
ious  types  of  insulators  and  bushings  operating  at  7.5KV,  15KV,  and  34.5KV 
and  exposed  to  various  lengths  of  carbon  fiber.  These  tests  which  are  summarized 
in  the  final  section  of  this  Appendix  provide  the  basic  inputs  to  this  analysis. 

Rationale 

The  risks  are  defined  in  terms  of  outage  frequencies  with  and  without 
exposure  to  carbon  fiber.  All  effects  of  carbon  fiber  exposure  are  estimated 
on  the  basis  of  3m.m-length  fiber  segments.  Power  systems,  typical  distribution 
systems,  and  some  protective  measures  are  first  described.  General  reliability 
considerations  together  with  outage  data,  are  then  discussed.  Finally,  esti¬ 
mates  of  the  effects  of  CF  exposure  are  made  on  the  basis  of  a  typical  circuit 
in  the  7.5KV  range  and  an  actual  23KV  circuit. 


A-2 


DESCRIPTION  OF  POWER  SYSTEMS  (References  1,  2,  3,  4) 

Illustrative  Power  System 

Figure  A.l  is  an  illustrative  power  distribution  system  consisting  of 
transmission  circuits,  subtransmission  circuits,  and  primary  feeders.  Sta¬ 
tions  A,  B,  and  C  are  bulk  power  stations  interconnected  by  a  transmission  net¬ 
work.  The  distribution  system  is  that  portion  of  the  power  system  which  inter¬ 
connects  the  customer  service  connections  with  the  bulk  power  sources.  In 
Figure  A.l,  this  includes  the  subtransmission  network.  Station  D,  and  the  primary 
feeders.  In  some  systems  radial  subtransmission  circuits  may  be  used  or  a 
loop  arrangement,  in  which  the  transmission  circuit  connects  (loops  to)  several 
distribution  stations  and  returns  tc  the  bulk  station,  may  be  used.  Often 
the  primary  feeders  are  connected  directly  to  a  bulk  station  without  use  of 
intermediate  subtransmission  voltages.  It  should  also  be  noted  that  major 
power  users  such  as  large  industrial  plants  may  be  fed  directly  from  trans¬ 
mission  or  subtransmission  circuits  rather  than  from  primary  feeders.  In  this 
case,  the  industrial  plants'  substations  are  equivalent  to  the  distribution 
substation. 

Circuit  breakers  are  usually  located  as  shown  except  that  the  breakers 
on  the  high  voltage  side  of  a  transformer  bank  are  often  omitted  because  of  the 
high  reliability  of  transformers.  Bus  tie  breakers  such  as  those  at  Station  A 
are  used  where  bus  sectional izing  in  the  event  of  a  fault  is  justified.  Bulk 
power  outages  are  usually  due  to  interruptions  in  the  transmission  network 
except  for  events  such  as  tornadoes  that  affect  large  portions  of  the  distri¬ 
bution  system.  Failures  due  to  CF  exposure  are  expected  to  occur  mainly  at 
the  primary  feeders;  however,  in  some  systems  the  subtransmission  circuits 
operate  at  voltages  below  38KV  and  might  also  be  susceptible. 

Illustrative  Distribution  Networks 

As  mentioned  above,  the  oistribution  system  generally  consists  of  dis¬ 
tribution  stations  connected  to  bulk  sources  via  subtransmission  circuits  and 
connected  to  the  customers  by  primary  feeders.  Subtransmission  circuits  may 
be  radial,  loop  or  network  configurations.  Likewise,  the  primary  feeders  may 
be  radial,  loop,  or  network  configurations.  The  distribution  transformers  may 
be  3  phase  for  industrial  or  commercial  establishments  or  may  be  single  phase- 
to-ground  to  supply  residential  circuits.  Secondary  voltages  are  usually 
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230/n5V  and  are  not  considered  in  this  analysis.  Figure  A. 2  shows  a  typical 
radial  distribution  feeder  in  some  detail  (from  Reference  4). 

A  fault  due  to  CF  exposure  on  a  primary  feeder  can  be  expected  to 
open  the  feeder  at  the  nearest  protective  device  (fuse  or  breaker)  toward  the 
source  (or  sources),  thus  cutting  off  all  loads  downstream  from  the  fuse  or 
breaker. 

Figure  3  shows  a  typical  distribution  circuit  for  single  family 
residences  in  an  urban  or  suburban  area.  A  single  phase  230/115V,  3-wire 
secondary  circuit  is  fed  from  a  three-phase  primary  feeder  by  a  distribution 
transformer  at  each  end  of  an  isolated  secondary  circuit.  Each  distribution 
transformer  is  connected  to  the  primary  feeder  through  a  fused  disconnect. 

The  isolated  secondary  circuit  may  be  extended  to  include  other  distribution 
transformers  connected  in  the  same  manner.  Other  single  phase  isolated 
secondaries  will  be  connected  to  other  phases  of  the  primary  feeder  in  order 
to  balance  the  loads  between  phases. 

In  some  residential  areas,  each  isolated  secondary  is  supplied  by  a 
single  transformer  and  often  in  rural  areas  a  distribution  transformer  will 
supply  a  single  residence.  A  common  distribution  system  for  commerical  areas 
uses  three-phase  transformers  supplying  each  establishment  directly  from  the 
primary  feeder,  as  shown  in  Figure  A. 3. 

Insulators  and  Bushings 

A  flashover  due  to  carbon  fibers  on  any  insulator  or  bushing  will 
cause  a  failure  to  the  entire  circuit  associated  with  that  insulator  or 
bushing.  The  circuit  in  this  case  is  defined  as  that  portion  of  the  system 
protected  by  an  associated  fuse  or  circuit  breaker.  For  the  case  of  a  pri¬ 
mary  feeder,  a  fault  on  any  phase  will  fail  all  three  phases.  For  the  case 
of  a  single-phase  fused  lateral  circuit,  a  fault  will  only  affect  the  par¬ 
ticular  single  phase  faulted. 

A  count  of  the  number  of  distribution  transformers,  disconnects, 
etc.  in  a  distribution  circuit  is  available  from  distribution  circuit  maps; 
however,  a  count  of  the  actual  number  of  pin  insulators  is  not  usually 
directly  available  (poles  are  usually  not  shown).  The  number  of  bushings  can 
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FIGURE  A. 3.  TYPICAL  DISTRIBUTION  CIRCUIT  ARRANGEMENT 


be  obtained  from  the  number  of  transformers  and  line  capacitors,  while  the 
number  of  post  type  insulators  can  be  obtained  from  the  number  of  disconnect 
switches  and  fuses. 

The  number  of  pin  (or  line-post)  insulators  is  proportional 
to  the  number  of  poles,  which  can  be  estimated  with  sufficient  accuracy  for  this 
study  on  a  per  mile  basis.  The  following  visual  observations  have  been  made 
in  the  Washington,  D.C,  area  for  urban,  suburban,  and  rural  primary  feeders 
owned  by  PEPCO,  VEPCO,  and  Baltimore  Gas  and  Electric  Company: 

•  In  urban  areas  there  are  about  60-80  poles  per  mile  and  in 
suburban  areas  there  are  about  30-40  poles  per  mile  with  pole 
spacings  dependent  on  the  location  of  customers  and  street  inter¬ 
sections.  In  rural  areas  the  poles  average  about  20-30  per  mile, 
depending  on  design  span  widths. 

•  Each  pole  has  1  pin  insulator  (or  line-post)  per  phase  of  primary 
feeder.  Each  distribution  transformer  has  1  vertical  bushing  per 
phase  plus  2  post  insulators  for  fused  cut-outs.  Each  capacitor 
bank  has  2  vertical  bushings  per  capacitor  and  each  disconnect 
switch  has  2  post  insulators  per  switch. 

PROTECTIVE  MEASURES 

The  protective  measures  include  devices  and  procedures  used  by  both 
the  utility  and  its  customers.  These  measures  are  of  interest  to  this 
problem  to  the  extent  that  they  affect  the  probability  of  an  outage  given 
a  failure,  the  duration  of  the  outage,  and  the  impact  of  an  outage  of  given 
duration. 

Generation/Transmission  Systems 

Protective  measures  at  this  system  level  are  only  of  background 
interest  to  this  analysis.  Faults  due  to  CF  on  distribution  circuits  in 
the  2.4KV-38KV  range  will  normally  be  cleared  at  the  distribution  level 
and  will  not  affect  the  generation/transmission  system.  Note  that  high 
load/long  duration  outages  affecting  the  generation/transmission  system 
are  included  in  the  bulk  outage  data  discussed  later. 
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Protection  measures  at  the  bulk  level  include  automatic  operation  of 
breakers  to  disconnect  faults,  load  shedding  and  generator  dropping  practices, 
schemes  for  bringing  spinning  reserves  onto  line,  means  for  obtaining  emergency 
start  up  in  the  event  of  system  outages,  etc.  The  automatic  opening  and  closing 
of  breakers  is  the  "first  line  of  defense"  in  the  event  of  unexpected  faults 
and  is  of  more  interest  to  this  study  than  other  protective  measures. 

Some  automatic  generation/transmission  system  protective  measures  are 
described  below  for  background  purposes: 

t  Generators  are  protected  from  internal  faults  by  means  of  percent 
differential  (current  balance)  relays  which  insure  that  all 
secondary  0-0  +  0-G  currents  are  equal.  Generators  are  also 
protected  for  overspeed,  overload,  anti -motoring,  loss  of  field, 
and  high  temperatures. 

•  Power  Transformers  are  protected  by  0-0  +  0-G  differential,  over¬ 
current,  and  thermal  relays. 

•  Station  Buses  have  percent  differential  protection  to  protect 
against  0-0  and  0-G  faults  anywhere  on  the  bus. 

t  Transmission  lines  are  protected  by  directional  overcurrent 
relays  and  impedance  (distance)  relays  at  the  station  at  each 
end.  These  relays  also  provide  backup  in  the  event  of  a  failure 
to  clear  faults  at  adjacent  stations  and  lines.  Breakers  are 
normally  programmed  to  open  in  about  5  cycles  and  reclose  in 
about  15  cycles  in  order  to  clear  lightning  strikes  without  power 
interruption. 

Distribution  System 

Subtransmission  lines  which  carry  high  distribution  loads  at  relatively 
high  voltages  will  be  protected  by  means  similar  to  the  bulk  transmission  lines. 
When  a  single  radial  subtransmission  line  feeds  a  distribution  station,  the 
breaker  at  the  distribution  station  end  of  the  line  may  be  omitted  with  protec¬ 
tion  provided  by  the  breaker  at  the  supply  end  only.  This  breaker  will  open 
for  any  faults  on  the  upstream  side  of  the  primary  feeder  breaker  and  will  pro¬ 
vide  backup  to  the  primary  feeder  breakers.  Primary  feeder  breakers  are  con¬ 
trolled  by  overcurrent  relays  for  each  phase  and  for  phase-to-ground  currents. 
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Three  automatic  reclosures  followed  by  a  lockout  (in  the  event  of  a  sustained 
fault)  is  common. 

Fuses  in  subfeeders  an'^  laterals  are  usually  designed  to  open  prior 
to  operation  of  the  primary  feeder  breaker  while  fuses  on  the  distribution 
transformers  are  designed  to  open  prior  to  damage  to  the  lateral  fuses,  etc. 
Therefore,  insulator  or  bushing  flashovers  due  to  carbon  fibers  will  affect 
only  that  portion  of  the  feeder  protected  by  the  next  protective  device  toward 
the  source  from  the  fault. 

Customer  Response  to  Power  Outages 

Tables  A.l  and  A. 2  from  Ref.  (5)  summarize  types  of  responses  of  resi¬ 
dential,  commercial,  and  industrial  customers  to  electric  utility  power  out¬ 
ages  based  on  experience  in  the  Pacific  Northwest. 

ORI  surveys  of  businesses  and  industries  showed  a  wide  variety  of 
protective  measures.  Facilities  in  which  power  continuity  is  absolutely 
necessary  (e.g.,  hospitals,  air  traffic  control  centers,  many  chemical  plants, 
radio  stations,  etc.)  have  auxiliary  generators  for  backup.  However,  other 
industries  we  surveyed  provided  no  backup  power  even  though  power  out¬ 
ages  lasting  more  than  a  few  minutes  would  be  costly.  Examples  include  a 
large  truck  engine  plant,  puolishing  companies,  etc.  The  rationale  for  this 
is  the  extremely  high  utility  reliability  that  these  facilities  have  experienced 
in  the  past  coupled  with  the  high  cost  of  auxiliary  power.  In  one  or  two  in¬ 
stances  in  which  the  locations  permitted,  large  industries  were  supplied  by 
more  than  one  primary  feeder  circuit. 

GENERAL  RELIABILITY  CONSIDERATIONS  AND  ESTIMATES  FOR  A  TYPICAL  DISTRIBUTION 
SYSTEM 

Reliability  Measures  (Ref.  5,  6) 

Commonly  applied  measures  of  power  systems  reliability  include 
the  Loss  of  Load  Probability  (LOLP)  and  the  Frequency  and  Duration  measure 
(FAD).  The  LOLP  accounts  for  the  total  fraction  of  time  that  a  power  sys¬ 
tem  is  expected  to  have  a  deficit  without  regard  to  the  distribution  of 
outage  durations.  The  FAD,  on  the  other  hand,  accounts  for  both  the  fre¬ 
quency  and  the  duration  of  the  outages.  It  is  desirable  to  account  also 
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TABLE  A.l 


SUMMARY  OF  POTENTIAL  CONSUMER  RESPONSES 
TO  LOWER  ELECTRIC  POWER  RELIABILITY 


Consumer 

Function 

Response 

Residential 

Heating 

Firewood  stored,  oil  or  gas  heat 

Refrigeration 

Dry  Ice 

Lighting 

Candles,  flashlights 

Cooking 

Camp  stove 

Commercial 

Lighting 

Batteries,  standby  generators 

Data  processing 

UPS  (Uninterrupted  Power  Supply) 
standby  generators 

Refrigeration 

Standby  generators 

Industrial 

Electric  drive 

Standby  generators 

Lighting 

Batteries 

Space  conditioning 

Nonelectric  heating  &  cooling 

TABLE  A. 2 

SUMMARY  OF  ESTIMATED  CONSUMER  RESPONSES  TO  LOWER  RELIABILITY 
LEVELS,  AS  A  FUNCTION  OF  FREQUENCY  AND  DURATION 


Total 

Consumer  Response 

Outage 

Industrial 

Commercial 

Residential 

144  min/ 
year 

Emergency  equipment 
installed 

Emergency  equipment 

No  response 

288  min/ 
year 

Non-electric  equipment 
installed 

Non-electric  equip¬ 
ment  installed 

No  response 

1440  min/ 
year 

Standby  generator 
installed 

Standby  generator 
installed 

Emergency 

equipment 

installed 

Source;  Impacts  from  a  Decrease  in  Electric  Power  Service  Reliability 
Stanford  Research  Institute,  June  1976. 
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for  the  class  of  customers  affected  (residential,  commercial,  industrial)  and 
for  the  size  of  the  loads  lost.  However,  in  general,  data  are  not  sufficient 
for  this  level  of  detail. 

Analysis  of  Bulk  Outage  Data  (Ref.  5,  7,  8,  9) 

Bulk  power  outages  are  defined  under  FPC  Order  331-1,  as  interrup¬ 
tions  of  a  generating  unit  or  electrical  facility  operating  at  a  nominal  volt¬ 
age  of  69KV  or  higher  and  resulting  in  a  load  loss  of  15  minutes  or  longer  of 
at  least  100  megawatts.  Smaller  systems  must  report  to  DOE  if  one  half 
or  more  of  the  annual  system  peak  is  involved.  (Ref.  5,  9). 

During  CY  1978  there  were  62  bulk  outages  reported  (Ref.  9).  These 
involved  a  loss  of  about  12,155  MW  and  65,000  MWH  to  about  3.1  million  cus¬ 
tomers.  The  distribution  of  outage  duration  by  number  of  outages,  customers 
affected,  customer-hours  lost,  and  loads  are  summarized  in  Table  A. 3.  Figure  A. 4 
shows  the  percent  distribution  of  customers  affected  by  various  duration  times. 
These  reflect  outages  only,  not  counting  load  reduction  measures.  When  ranges 
of  outage  times  were  given,  a  midpoint  value  was  used. 

fi  6 

If  it  is  assumed  that  all  63.4  x  10  households  and  4.1  x  10  com¬ 
mercial  business  establishments  are  utility  customers,  then  about  4.6%  of  all 
customers  were  affected  by  bulk  outages  in  1978.  This  neglects  the  fact  that 
some  customers  were  affected  by  more  than  one  bulk  outage  during  1978. 

Reference  9  reports  the  following  distribution  of  bulk  outages  by  customer 
types : 

Residential  55% 

Commercial  30% 

Industrial  25% 

The  number  of  feeders  required  in  an  area  is  primarily  a  function 
of  the  peak  loads,  with  a  typical  feeder  handling  2-3MW.  Bulk  outages  will 
usually  drop  a  large  number  of  feeders,  and  the  fraction  of  feeders  dropped 
is  roughly  proportional  to  the  fraction  of  loads  dropped.  If  it  is  assumed 
that  bulk  outages  are  randomly  distributed  among  feeders,  then  the  annual 
likelihood  of  a  particular  feeder  being  dropped  is  approximately  equal  to 
the  fraction  of  feeders  dropped  per  year.  At  an  average  U.S.  load  of  about 
350,000MW,  the  12,157MW  load  dropped  due  to  1978  bulk  outages  represents 
about  3.4%  of  the  feeders  so  that  the  likelihood  of  a  feeder  being  dropped 
due  to  a  bulk  outage  is  approximately  .034  per  year. 
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TABLE  A. 3.  DISTRIBUTION  OF  BULK  POWER  OUTAGES  CY  1978 


OUTAGE 

NO. 

CUSTOMERS 

CUSTOMER- 

LOAD 

DURATION 

OF 

AFFECTED 

HOURS  LOST 

LOST 

(Hours-Days) 

OUTAGES 

(Thousands) 

(Thousands) 

(MW) 

0-1  Hours 

35 

1612 

895 

6879 

1-3 

13 

410 

465 

2635 

3-6 

5 

344 

2217 

1459 

6-12 

4 

145 

1722 

678 

12-24  " 

3 

72 

759 

183 

1-5  Days 

1 

200 

9600 

323 

5-10  " 

1 

300 

43,200 

Not  Available 

TOTAL 

62 

3083 

58,858 

12,157 

FIGURE  A. 4.  PERCENTAGE  OF  CUSTOMER-OUTAGES  OF  DURATION  D 
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Outages  Due  to  Normal  Equipment  Failures 

Estimates  of  outages  due  to  normal  equipment  failures  (not  in  the 
presence  of  carbon  fibers  )  can  be  made  by  applying  equipment  failure  data 
to  "typical"  distribution  systems. 

Figure  A. 5  shows  a  one-line  diagram  of  a  "typical  distribution  cir¬ 
cuit.  This  circuit  which  feeds  an  urban  residential  area  consists  of  3 
circuit-miles  of  three-phase  primary  feeder,  16  circuit  miles  of  three- 
phase  subfeeders,  and  8  circuit  miles  of  single-phase  lateral  feeders. 

These  circuits  supply  4000  homes  through  398  distribution  transformers. 
Sectional izing  fuses  are  provided  in  each  subfeeder  and  lateral  and  at 
the  load  midpoint  of  the  primary  feeder.  One  or  more  normal ly-open  manual 
tie-switches  permit  interconnection  to  other  feeders  in  the  event  of  an 
emergency. 

The  expected  outages  (Eq)  per  customer  (or  per  distribution 
transformer)  can  be  estimated  from: 

where : 

E^  =  outages  in  section  i 

F^.  =  fraction  of  customers  (or  transformers)  affected  by  an 
outage  of  section  i 

and: 

Ei  =  X,  (t) 

=  failure  rate  of  equipment  in  section  i 
t  =  time  over  which  outages  are  estimated 
The  probability  of  an  outage  in  section  i  is  given  by: 

P.  =  1-  Exp  (-E^.) 

Failure  rates  and  average  outage  durations  are  shown  in  Table  4 
from  the  sources  indicated.  It  should  be  noted  that  there  are  large 
uncertainties  and  variations  in  both  published  failure  rates  and  out¬ 
age  durations.  For  example,  the  rates  shown  in  Table  A. 4  for  open  line 
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SECTIONALIZING 
FUSE  - 


— #\#- 


► 


T 

1  MILE 


SECTION 

A 


1  MILE 


SECTION 

B 


MILE 


TO 

OTHER  -4 
CIRCUIT 


NORMALLY 

OPEN 


/ 


•  2500  KVA  PRIMARY  FEEDER  -  SINGLE  FED  RADIAL 


•  4  SO  MILE  RESIDENTIAL  AREA,  4000  HOMES 

•  16  SUBFEEDERS  EACH  1  MILE  LONG 


•  16  LATERALS  EACH  1/2  MILE  LONG 

•  15-5  KVA  DISTRIBUTION  TRANSFORMERS/MILE  ON  EACH  SUBFEEDER,  EACH  LATERAL,  AND  ON  PRIMARY 
FEEDER  PAST  FEED  POINT. 

•  SUBFEEDERS,  LATERALS,  AND  TRANSFORMERS  ALLOCATED  EQUALLY  TO  SECTION  A  &  B. 


FIGURE  A. 5.  "TYPICAL"  DISTRIBUTION  CIRCUIT 
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TABLE  A. 4 


FAILURE  RATES  OF  VARIOUS  DISTRIBUTION  EQUIPMENTS 


DISTRIBUTION 

EQUIPMENT 

NO. 

FAILED 

UNIT- 

HOURS 

FAILURE 

RATE 

OUTAGE 

DURATION 

SOURCE 

Dist.  Sta. 
Transformer 

- 

.012/Yr 

T3 

Feeder 

Cir.  Breaker 

20 

1870 

.011/Yr 

2,3 

Distribution 

Sta.  Bus 

35 

584 

.06/Yr 

5  Hr 

Distribution 

Regulator 

10 

280 

.036/Yr 

4  Hr 

2,  3 

Primary 

Feeder 

- 

- 

.07/Mi/Yr 

3  Hr 

1,  3 

Lateral 

Feeder 

- 

- 

.18/Mi/Yr 

2  Hr 

1,3 

Line 

Fuse 

- 

- 

.0007/Yr 

1  Hr 

3 

Distribution 

Transformer 

28 

290 

.097/Yr 

2  Hr 

2 

SOURCES: 

1  —  Power  Systems  Reliability  Calculation  (Ref.  10) 

2  —  Determination  and  Analysis  of  Data  for  Reliability  Studies  (Ref.  1 1 ) 

•  Field  Data  from  Texas  Electric  Service,  Co. 

3  —  Reliability  Information  for  Electric  Utility  Transmission  and  Distribution  Systems 

(Ref.  6) 
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feeders  apparently  include  little  if  any  consideration  of  severe  wind,  ice, 
or  lightning  conditions.  Table  A. 5  provides  equipment  outage  rates  for  the 
various  equipments  shown  in  Figure  A. 5. 

If  there  were  no  sectional i zing  fuses,  all  customers  would  be  affected 
by  all  outages  and  the  average  outage  per  customer  would  be  the  same  as  the 
total  for  the  circuit,  2.859/year.  The  annual  outage  probability  per  customer 
would  be: 

1  -  Exp  (-2.859)  =  0.94. 

Now  if  the  primary  sectional i zing  fuse  is  added  and  coordinated 
so  as  to  open  for  faults  in  Section  B  without  interrupting  customers  in 
Section  A,  and  if  the  subfeeder  fuses  and  lateral  fuses  are  omitted,  then 
all  customers  would  be  affected  by  faults  in  Section  A  while  only  50%  of 
the  customers  would  be  affected  by  faults  in  Section  B.  This  results  in: 

Section  Outages/Year  Outage  Probability  Customers  Affected 

A  1.50  0.78  100% 

B  1.36  0.74  50% 

The  average  annual  outages  per  customer: 

=  1 .50  X  100%  +  1  .36  X  50%  =  2.18 

The  annual  outage  probability  per  customer  is  obtained  from: 

(100%)  +  (50%)  +  P^g  (100%)  =  (.78)  (1-.74)  + 

(.74)  (1-.78)  (.50)  +  (.78)  (.74)  +  0.86. 

If  fuses  are  now  added  to  each  subfeeder  and  each  lateral  the 
following  allocation  of  failure  rates,  outage  probability,  and  customer 
outages  results: 
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TABLE  A. 5 


ANNUAL  OUTAGES  FOR  TYPICAL  DISTRIBUTION 


DISTRIBUTION  STA.  TRANSFORMER 
CIRCUIT  BREAKER 
DIST.  FEEDER  REGULATOR 
MAIN  FEED  (1  MILE  OPEN  LINE) 

SEC  A  PRI  FEED  (1  MILE  OPEN  LINE) 

SEC  A  SUBFEEDERS  (8  MILE  OPEN  LINE) 
SEC  A  LATERALS  (4  MILE  OPEN  LINE) 

SEC  A  FUSES  (12  TOTAL) 

TOTAL  SEC  A 


SEC  B  PRI  FEED  (1  MILE  OPEN  LINE) 

SEC  B  SUBFEEDERS  (8  MILE  OPEN  LINE) 
SEC  B  LATERALS  (4  MILE  OPEN  LINE) 
SEC  B  FUSES  (13  TOTAL) 


TOTAL  SEC  B 


TOTAL  FOR  CIRCUIT 


CIRCUIT 


=  .012/YR 

=  .011/YR 

=  .060/YR 

=  .070/YR 

=  .070/YR 

=  .560/YR 

=  .720/YR 

=  .008/YR 

=  1.500/YR 


=  .070/YR 

=  .560/YR 

=  .720/YR 

=  .009/YR 

=  1.359/YR 


=  2.859/YR 
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Portion  of 

Outages  per 

Outage 

Customers 

System 

Year 

Prob. 

Affected 

Equip,  prior  to 

Pri  fuse,  P 

2.12 

0.19 

100% 

Primary  feeder 
Section  B,  b 

.07 

0.07 

50% 

Each  subfeeder,  S 

.07 

0.07 

40% 

Each  lateral,  L 

.09 

0.09 

2% 

The  average  number 

of  annual  outages  per  customer  is: 

.212  (100%)  +  .07  (50%)  +  16  (.07)  (4%)  +  16  (.09)  (2%)  =  .32. 


The  annual  outage  probability  per  customer  is  approximately: 

1-  Exp  (-0.32)  =  0.27 

Addition  of  bulk  outages  to  equipment  outages  results  in: 

Circuit  Configuration  Annual  Outage  Prob.  per  Customer 

^^"^Bulk  •  ^Equip.) 


No  Fuse 

0.95 

Primary  sectionalizing  fuse 

0.87 

Subfeeder  and  Lateral  f^L>es 

0.30 

Note  that  these  values  apparently  understate  the  outages  due  to  tornadoes, 
ice  storms,  etc. 

Distribution  Outage  Data 

Outage  data  at  the  distribution  circuit  level  are  not  generally 
available  in  the  published  literature.  Discussions  with  various  utility 
companies  indicate  that  about  5-10  outages  per  year  per  distribution  cir¬ 
cuit  can  be  expected.  The  23KV  distribution  circuit  described  later  had 
8  outages  resulting  in  a  total  of  about  3700  customer-outages  in  10  calen¬ 
dar  months. 

EFFECTS  OF  CF  ON  DISTRIBUTION  CIRCUITS 

The  exposure  of  electric  utility  distribution  systems  to  high  con¬ 
centrations  of  carbon  fibers  is  expected  to  result  in  flashovers  across 
insulators  and  bushings.  Since  a  large  number  of  insulators  will  be  ex¬ 
posed  over  a  period  of  time,  a  series  of  interruptions  can  be  expected 
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to  occur.  If  the  time  between  flashovers  exceeds  the  feeder  breaker  re¬ 
closing  cycle  times  (adjustable  from  a  few  seconds  to  about  6  minutes),  a 

series  of  intermittent  interruptions  will  occur.  If  a  series  (burst)  of  flash¬ 
overs  spans  the  reclosing  time,  the  feeder  breaker  will  lock  open.  In  this 
analysis,  exposure  to  CF  is  assumed  to  result  in  an  interruption  with  a 
probability  determined  by  the  Weibull  distribution. 

The  exposure  values  expected  to  result  in  flashovers  are  obtained 
from  Westinghouse  Electric  Company  tests  which  indicate  that  the  percentage 
of  insulators  failed  versus  the  exposure  values  can  be  approximated  by  the 
Weibull  function.  Some  representative  values  from  the  Westinghouse  data 
are  contained  in  Attachment  A.  Note  that  large  differences  in  failure 
characteristics  result  from  tests  on  different  types  of  insulators  and 
different  applied  voltages. 

In  the  paragraphs  which  follow,  a  typical  distribution  circuit  at 
7.5KV  primary  voltage  and  a  selected  actual  distribution  circuit  operating 
at  23KV  are  analyzed.  The  closest  applicable  insulation  types  and  voltage 
are  selected  from  the  Westinghouse  data  using  the  tests  of  2mm  fiber 
lengths;  the  exposure  values  at  2mm  are  then  linearly  extrapolated  to 
values  at  the  3-mm  fiber  lengths  assumed  for  this  study. 

The  maximum  exposure  value  expected  to  be  encountered  are  less 
than  1  X  10^  fiber-sec  per  cubic  meter  which  yield  per-insulator  failure 

O 

probabilities  on  the  order  of  1  x  10"  .  The  insulator  types  which  occur 
in  small  numbers  in  a  circuit  (e.g.  circuit  breaker  bushings)  can  therefore 
be  neglected.  Exposure  values  are  translated  to  insulator  failure  probabili¬ 
ties  using  the  Weibull  function  with  input  constants  from  the  Westinghouse 
data. 

Estimates  of  customer  outages  are  made  based  on  accidents  at  the 
Los  Angeles  International  Airport  assuming  (a)  first  that  all  circuits 
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fit  the  7.5KV  circuit  characteristics  and  then  (b)  that  all  circuits  fit 
the  23K\/  circuit  characteristics. 


Analysis  of  7.5  KV  Distribution  System 

The  7.5KV  distribution  system  was  described  previously  and  shown 
in  Figure  5.  The  failure  probabilities  are  dominated  by  the  pin  insula¬ 
tors.  If  the  effects  of  sectional izing  are  neglected,  which  represents 
the  worst  case,  all  pin  insulators  can  be  considered  in  series.  The 
total  number  of  insulators  are  estimated  as  follows: 


Primary  feeder  -  2  miles  x  240  insulators/mi .=  480 

Secondary  feeders  -16  miles  x  240  insulators/mi .=  3840 

Lateral  feeders  -  8  miles  x  80  insulators/mi .=  640 

Total  4860 

Therefore  5000  insulators  are  assumed  to  be  in  the  circuit. 


function: 


Insulator  failure  probabilities  arey estimated  from  the  Weibull 


Ppj  =  1  -EXP 


cto 


a 


"X 


Ppj  =  Prob.  of  failure  of  a  single  insulator 


^0  =  Exposure  value  in  fiber-sec  per  cubic  meter 
ao  =  Cut  off  point  on  Kq 

a,  3 =  Weibull  constant  for  the  applicable  insulator  test 
data 

VIq  =  y  at  fiber  length  in  test,  y^  =  y  of  desired  fiber  length 


The  probability  of  failure  (outage)  of  an  entire  circuit  with  n 
insulators  in  series  is  equal  to: 


PFC  =  1  - 


g 

The  Westinghouse  data  for  the  wet  7.5KV  pin  insulator  shows:  a=.69  x  10  , 

3=7.4,  y^  at  2mm  =  5.77  x  10^.  Using  linear  extrapolation  (see  attachment  A) 

the  y  a  3mm  =  3.10  x  10^.  Applying  the  Weibull  function  with  these 
X 
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values  at  =  0  results  in  the  value  of  failure  probability  versus  exposure 
shown  in  Table  A. 6. 

Accident  Scenario.  There  are  expected  to  be  2.6  aircraft  accidents 
per  year  which  will  result  in  fiber  release.  Previous  ORI  work  has  shown  that 
a  large  release  consists  of  5  x  10^^  fibers.  At  stability  class  6  and  a  5.5 
meter/sec  wind  speed  the  exposure  footprints  shown  in  Figure  A. 6  will  result 
from  each  release  of  5  x  10^"*  fibers.  Based  on  our  earlier  analyses,  there 
conditions  represent  the  worst  case. 

Los  Angeles  International  Airport  is  chosen  as  an  accident  site 
which  represents  a  severe  case  from  the  standpoint  of  downwind  population  den¬ 
sities.  Data  for  downwind  areas  are  shown  in  Table  A. 7.  Each  household,  busi¬ 
ness,  and  industry  is  assumed  equal  to  one  utility  customer. 

The  expected  number  of  customer  outages  may  now  be  obtained  as 

follows : 

Nq  =  Ng  ^  (PC.)  (Di)  (Ai) 

where: 

=  Number  of  customer  outages  per  year 
=  Number  of  accidents  per  year 

PC^.  =  Probability  of  an  outage/customer  for  exposure  to  level  i 

Di  =  Customer  density  in  exposure  level  i 

Ai  =  Area  covered  by  exposure  level  i  (per  accident) 

The  values  of  PCi ,  Di  and  Ai  are  shown  in  Table  A. 8.  Intersected  areas  were 
obtained  by  manually  comparing  footprints  versus  demographic  areas  from  Table  7. 
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TABLE  A. 6 

OUTAGE  PROBABILITY  VS.  EXPOSURE  LEVEL  FOR  7.5KV  CIRCUIT 


EXPOSURE 

LEVEL 

FAILURE  PROBABILITY 

1  INSULATOR 

ENTIRE  CIRCUIT 

1  X  103 

7.29  X  10-” 

5.00  X  10-2 

5x  103 

3.64  X  lO-lO 

2.00  X  10-O 

1  X  10^ 

7.29  X  10-10 

3.50  X  10-0 

5x  lO'* 

3,64  X  10-O 

1.80  X  10-0 

1  X  105 

6.16  X  10-9 

3.10X  10-0 

2.5  X  10^ 

2.00  X  10-0 

1.05  X  10-^ 

5x  105 

3.64  X  10-0 

1.80  X  10-1 

1  x  106 

6.15  X  10-0 

3.03  X  10-2 

1  X  10^ 

6.15  X  10-0 

2.65  X  10-1 

2  X  10^ 

1.03  X  10-2 

0.9999  .  .  . 

3x10^ 

1.88  X  10-1 

0.9999  ,  .  . 
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FIGURE  A. 6.  EXPOSURE  FOOTPRINTS  FOR  RELEASE  OF  5x10"  SINGLE  FIBERS 
AND  METEOROLOGICAL  CONDITIONS  SHOWN 


TABLE  A. 7 


DEMOGRAPHIC  DATA  FOR  AREA  DOWNWIND  OF  LOS  ANGELES  INTERNATIONAL  AIRPORT 


COUNTY  OR  CITY 

AREA 
SO.  Ml. 

DISTANCE 
FROM  A/P 

DISTANCE 

ALONG 

DOWNWIND 

AXIS 

CUSTOMERS 
PER  SO.  Ml. 

Los  Angeles 

SMSA 

4069 

0  Mi. 

35  Mi. 

650 

Orange 

County 

782 

35  '' 

30  " 

615 

Riverside 

City 

71.3 

40  '' 

10  '' 

230 

Ontario 

City 

32.1 

40 

10  '' 

671 

San  Bernadino 

County 

20,117 

30  " 

10  " 

Urban  -  370 
Other  -  15 

Riverside 

County 

7176 

40  " 

180 

22 

TABLE  A. 8 

INTERSECTED  AREAS  (Ai),  CUSTOMER  OUTAGE  PROBABILITY  (PCi),  AND  CUSTOMER  DENSITY 
(Di)  FOR  INTERSECTION  OF  EXPOSURE  FOOTPRINTS  WITH  AREAS  FROM  TABLE  8. 


EXPOSURE 

LEVEL 

Eo 

FOOTPRINT 
SIZE  (Mi.) 

DEMOGRAPHIC 

REGION 

INTERSECTED 

INTERSECTED 

AREA 

Ai 

CUSTOMER 

OUTAGE 

PROB. 

PCi 

CUSTOMER 

DENSITY 

Di 

L 

W 

Area* 

5x  105 

30 

4 

37 

Los  Angeles,  SMSA 

37  Sq.  Mi. 

1.80  X  10-4 

650 

2.5  X  10® 

66 

6 

83 

Los  Angeles 

5Sq.  Mi. 

1.05  X  10-4 

650 

San  Bernadino 
County  -  Urban 

5  Sq.  Mi. 

1.05  X  10-4 

370 

Riverside  City 

71  Sq.  Mi. 

1.05  X  10-4 

230 

Riverside  County 

2  Sq.  Mi. 

1.05  X  10-4 

22 

5x  ICH 

192 

20 

1093 

San  Bernadino 
County  “  Other 

365  Sq.  Mi. 

1.80  X  10-5 

15 

Riverside  County 

728  Sq.  Mi. 

1.80  X  10-5 

22 

5x  103 

240 

36 

1557 

San  Bernadino 
County  -  Other 

519  Sq.  Mi. 

2.00  X  10-5 

15 

Riverside  County 

1038  Sq.  Mi. 

2.00  X  10-6 

22 

‘Areas  are  those  Portions  of  Footprints  that  are  Mutually  Exclusive, 
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Applying  the  above  relations 

results  in  the  following 

customer 

outages 

per  accident: 

Customer 

Exposure 

Level  Outage/Customer 

X  Density 

X  Area  = 

Outages 

5x10^ 

1.80x10"^ 

650 

37 

4.330 

2.5x10^ 

1.05x10'^ 

650 

5 

0.340 

2.5x10^ 

1 .05x10"^ 

370 

5 

0.190 

1.05x10'^ 

230 

71 

1 .710 

1.05x10”^ 

22 

2 

0.005 

5x10^ 

1.80x10"^ 

15 

365 

0.099 

5x10^ 

1 .80x10'^ 

22 

728 

0.288 

5x10^ 

2.00x10”^ 

15 

510 

0.016 

5x10^ 

2.00x10"^ 

22 

1038 

0.046 

Total 

per  accident 

^  9.000 

Total  per  year  =  2.6  x  9  -  23.40 

This  result  compares  with  over  3  million 

customer  outages  due  to  bulk 

outages 

alone.  Applying  the  results 

of  the  reliability  analysis  for  the 

typical 

circuit  under  the  same  fusi ng 

conditions  as 

above  yields 

at  least 

another 

20  million  outages  per  year  due  to  normal 

distribution 

system  failures 

Analysis 

of  a  Selected  23KV  Distribution  Circuit 

Figure  A. 7  shows  a  one-line  schematic  of  an  actual  23KV  distribu¬ 
tion  circuit.  This  circuit  serves  about  1800  industrial,  commercial,  and 
residential  customers  in  an  urban  area  of  about  one  square  mile.  This 
circuit  consists  of  a  3(J)  radial  primary  feeder  with  3(p  subfeeders  and 
laterals  together  with  a  few  ItJ)  laterals  as  shown.  There  is  no  automa¬ 
tic  sectional izing  capability  except  for  the  fusing  of  some  laterals  as 
shown.  There  exists  the  capability  to  cross-tie  to  adjacent  circuits  at 
several  points. 
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FIGURE  A. 7 


300'  200’ 


FROM 
SHEET  B 


SYMBOLS 

30 

-►  3-Phase  subfeeder  or  lateral 

XOO'  is  the  length  of  the  subfeeder  or  lateral  open  wire 
circuits 

U.G.  —  Underground 

X  (3)  —  Number  of  30  transformer  in  subfeeder  or  lateral 

X  (1)  ~  Number  of  10  transformer  in  subfeeder  or  lateral 

fX  (1)  —  Number  of  10  transformers  located  on  primary 
feeder 

NC  —  Normally  closed 
NO  “  Normally  open 


FIGURE  A. 7.  SHEET  C 
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Circuits  and  transformers  are  summarized  below: 


Total  3<()  open  lines  -  37,000  ft. 

Total  l(j)  open  lines  10,300  ft. 

Total  exposed  distribution  transformers 
3(t)  ^  32 
Icj)  ^  116 

Total  non-exposed  distribution  transformers 
(in  underground  vaults  and  in  buildings) 

3cJ)  =  32 

l(j)  =  40 

Normal  Reliability.  There  are  about  8-10  normal  outages  (about 
3700  customer-outages)  per  year  on  this  circuit.  This  appears  to  be  typi¬ 
cal  for  the  type  of  circuit  (based  on  discussions  with  several  utilities). 

Insulator  Failure  Probability.  The  total  number  of  insulators 
are  estimated  from: 

7  miles  of  3(1)  lines  x  240/mile  =  1680 
2  miles  of  Icf)  lines  x  80/mile  =  160 
Total  =  1840 

The  34.5KV  distribution  post  insulators  are  the  closest  applicable 
insulators,  used  as  line  posts  for  the  23  KV  open  lines.  The  Westinghouse 
data  for  this  shows  a  =  .19  x  10^,  6  =  1.99,  M=2mm  =  1.8  x  10®,  M=3mm  (by 

O 

linear  extrapolation)  =  1.0  x  10  .  Failure  probabilities  for  a  single 
insulator  and  for  all  insulators  in  series  are  shown  in  Table  A. 9. 

Accident  Scenario.  The  same  accident  scenario  is  applied  as 
used  previously,  2.6  aircraft  accidents  each  releasing  5  x  lO^^fibers  at 
Los  Angeles  International  Airport,  and  producing  the  same  footprints. 
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TABLE  A. 9 


OUTAGE  PROBABILITY  VS.  VARIOUS  EXPOSURE  LEVELS  FOR  23  KV  CIRCUIT 


EXPOSURE 

LEVEL 

FAILURE  PROBABILITY 

1  INSULATOR 

ENTIRE  CIRCUIT 

1  X  103 

I.Ox  10-10 

1.84  X  10-’ 

5x  103 

3.0  X  10-0 

4.60  X  10-8 

1  X  104 

I.Ox  10-8 

1.80  X  10-8 

5x  10^ 

2.4  X  10-’ 

4.46  X  10-4 

1  X  10^ 

9.6  X  10-’ 

1.42  X  10-3 

2.5  X  10^ 

5.9  X  10-O 

1.09  X  10-2 

5x  105 

2.4  X  10-8 

4.26  X  10-2 

1  X  106 

9.4  X  10-8 

.1589 

1  X  10’ 

9.1  X  10-3 

.9999 

2  X  10’ 

3.5  X  10-2 

.9999 

3x  10’ 

7.8  X  10-2 

.9999 
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The  Expected  Number  of  Customer  Outages  are 

estimated  in 

the  same 

way  as  was  done 

previously,  producing 

the  following 

results : 

Customer 

Exposure  Level 

Outage/Customer 

X  Density  x 

Area 

Outages 

5x10^ 

4.26x10“^ 

650 

37 

1024.53 

2.5x10^ 

1.09x10'^ 

650 

5 

35.43 

1 .09x10'^ 

370 

5 

20.17 

1 .09x10'^ 

230 

71 

178.00 

1 .09x10"^ 

22 

2 

0.48 

5.0x10^ 

4.46x10"^ 

15 

365 

2.44 

4.46x10"^ 

22 

728 

7.14 

5.0x10^ 

4.60x10"^ 

15 

519 

0.04 

4.60x10"^ 

22 

1038 

0.11 

Total  per  accident 

1268.34 

Total  per  year  =  2.6  x  1268.34  -  3298 
customer  outages  per  year 


From  the  actual  circuit  outage  data  there  were  3700  outages  for 
1800  customers  or  2.06  outages  per  year  per  customer.  For  the  97636 
customers  within  the  exposure  footprints,  normal  customer  outages  are 
expected  to  be  about  2.06  x  97636  =  2.01  x  10  customer  outages.  The 
ratio  of  CF/normal  outages  is:  3298  t  2.01  x  10^  ^  16/1000  for  those 
affected  by  CF.  Over  the  entire  U.S.,  the  3298  CF  induced  outages/year 
compares  with  over  20  million  which  are  expected  to  occur  normally. 

CONCLUSION 

The  carbon  fibers  appear  to  present  an  insignificant  problem  to 
electrical  distribution  circuits,  even  when  estimates  are  made  using 
very  pessimistic  (worst  case)  assumptions. 


A-32 


APPLICATION  OF  INSULATOR  TEST  DATA  TO  POWER  DISTRIBUTION  STUDY 


The  Westinghouse  Electic  Company  has  performed  tests  on  various 
types  of  insulators  and  bushings  at  various  appl ied  voltages  and  exposed 
to  various  fiber  lengths  .The  results  of  these  tests  are  summarized  in 
Table  A. 10  (Westinghouse  Table  3.3).  Figure  A. 8  (Westinghouse  Figure 
3.3-1)  shows  interpolation  of  some  of  these  data  for  various  fiber 
lengths. 


Figure  A. 8  has  been  used  to  scale  values  of  F  to  3mm  fiber  lengths 
for  those  insulators  shown  on  Figure  3.3-1.  Insulators  and  bushings  not 
shown  on  Figure  A, 8  are  assumed  to  scale  to  fiber  length  at  the  same 
rate  as  the  15  KV  C  Neck  Distribution  post  since  this  is  the  most  con¬ 
servative  rate  indicated  on  Figure  A. 8.  The  results  are  shown  in  Table 
A. 11 . 


The  values  of  F  at  3mm  are  selected  for  the  insulator  types  and 
voltages  closest  to  those  representing  the  distribution  systems  being 
analyzed. 


Figure  A. 9  shows  the  values  of  exposure-to-flashover  for  a  wet 
7.5KV  pin  insulator  while  Figure  A. 10  shows  the  values  of  exposure-to- 
flashover  for  34.5KV  distribution  post.  Thse  figures  and  associated 
Weibull  constants  are  used  as  the  basis  for  failure  estimates  made  in 
this  study  (after  translation  of  mean  values  for  3mm  fiber  lengths). 
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TABLE  A.l  RESULTS  OF  TESTS  FOR  SELECTED  INSULATORS 
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FIGURE  A. 8  LINEAR  INTERGEOLATION  OF  EXPOSURE  TO  FIBER  LENGTH 


TABLE  A. 11 

ESTIMATES  OF  E  AT 
3  MM  FIBER  LENGTH 


INSULATOR 

E 

7.5  KV  Pin  Insulator  (Dry) 

1.5  X  10^ 

7.5  KV  Pin  Insulator  (Wet) 

3.1  X  10^ 

15  KV  C-Neck  D.  Post  Vert. 

2.5  X  10^ 

1 5  KV  C-Neck  D.  Post  Vert  (Wet) 

2.1  X  lO'^ 

34.5  KV  C-Neck  D.  Post  Vert 

CO 

o 

X 

7.5  KV  Station  Post  Vert. 

4.4  X  10^ 

7.5  KV  Station  Post  Hor. 

1.2  X  109 

5  KV  Trans.  Bushing  Vert. 

9.5  X  10^ 

15  KV  Trans.  Bushing  Vert. 

2.5  X  lO'^ 

5  KV  Trans.  Bushing  Hor. 

5.8  X  10^ 

7.5  KV  Suspension  Vert. 

2.0  X  10^ 
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TOTAL  AIRCRAFT  ACCIDENT  COSTS 


APPENDIX  B 


TOTAL  AIRCRAFT  ACCIDENT  COSTS 


ORI  utilized  a  1978  FAA  report—'^  and  historical  records  from  the 
National  Transportation  Safety  Board  (NTSB)  in  determining  the  costs  of  major 
airplane  accidents  of  the  past  ten  years.  These  costs  were  used  as  the  basis 
for  a  comparison  of  the  additional  risks  (costs)  presented  by  the  use  of 
carbon  fibers  with  the  accepted  costs  associated  with  major  airplane  accidents 
of  the  recent  past. 

The  ORI  costing  methodology  considered  only  two  items:  aircraft 
hull  damage  and  costs  of  personal  injuries  to  crew,  passengers  and  persons 
on  the  ground.  Other  costs  such  as  investigative  costs,  property  damage  to 
ground  structures,  and  other  incidental  costs  were  not  included. 

Subsets  of  the  entire  NTSB  aircraft  accident  file  (1966-1975)  were 
drawn  in  order  to  compare  carbon  fiber  related  costs  to  the  costs  of  an 
appropriate  population  of  past  accidents.  It  was  decided  that  this  appropriate 
population  would  be  those  earlier  accidents  which  were  on  similar  scales  in 
terms  of  aircraft  type  and  accident  severity  as  the  future  accidents  could 
conceivably  result  in  carbon  fiber  release.  Only  accidents  which  involved 
U.S.  commerical  transport  jets  of  significant  size  were  included  in  the  analysis. 
These  aircraft  closely  approximate  the  types  of  planes  which  will  be  flying 
with  carbon  fibers.  All  the  recorded  accidents  involved  one  of  the  following 
aircraft: 


1_/  William  L.  Fallon,  Cost  Analysis  of  Aircraft  Accidents,  FAA  Office  of 
Aviation  System  Plans,  June  1978. 
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Small  Jets 
(150  passengers) 


Medium  Jets 
(150-260  passengers) 


Large  Jets 
(260  passengers) 


BAC-111 

B-707 

B-720 

B-727 

B-737 

DC-8  (non  60's  series) 
DC- 9 


DC-8-61 

DC-8-62 

DC-8-63 


B-747 
DC-10 
L  1011 


In  addition  to  aircraft  size  and  type,  the  severity  of  the  accidents  was 
taken  into  account.  Only  those  NTSB  file  entries  in  which  the  fuselage  had 
been  labelled  as  being  "destroyed"  or  having  "substantial"  damage  were  retained. 
This  strategy  was  required  in  order  to  exclude  minor  incidents  contained  in 
the  NTSB  file.  These  less  calamitous  accidents  may  be  as  small  as  a  person¬ 
al  injury  caused  during  strong  air  turbulence  and  were  excluded  so  as  not  to 
reduce  the  significance  of  the  individual  accidents  in  which  partial  or  total 
fuselage  destruction  occurred.  The  restrictions  these  limitations  produced 
reduced  the  original  1966-1975  file  size  from  560  accessible  commercial 
aircraft  accident  files  to  155  jet  accidents  with  the  proper  accident 
severity  levels  (16  large  jets,  0  medium  jets,  149  small  jets). 

For  hull  cost  computation,  the  selling  price  of  an  identical  air¬ 
craft  in  the  accident  year  was  culled  from  printouts  documenting  the  cited 
FAA  study.  The  original  source  of  this  information  was  publically  available 
records  of  commercial  transactions  such  as  those  contained  in  Aviation  Week. 

No  individual  consideration  was  given  to  special  avionics  or  other  equipment. 

One  hundred  percent  of  the  replacement  cost  was  assigned  to  those  hulls  which 
were  "destroyed";  one  third  of  the  ocst  was  assigned  to  the  substantially 
damaged  aircraft.  The  FAA  developed  the  one  third  factor  through  consulta¬ 
tions  with  NTSB  and  industry  experts. 

On  the  subject  of  personal  injury  costs,  $300,000  was  selected  as  the 
the  cost  of  a  fatal  injury. 

The  $300,000  was  based  on  the  projection  of  non-Warsaw  Pact  air¬ 
craft  accident  claims  settlements,  as  reported  by  the  Civil  Aeronautics  Board. 
The  figure  has  been  endorsed  by  the  Associated  Aviation  Underwriters,  and  is 
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used  by  the  agency  (FAA)  in  its  cost  benefit  work  in  facility  establishment 
criteria.  Serious  injury  is  $45,000  based  on  |ctual  settlement;  minor  injury 
cost  is  $6,000  based  on  recognized  methodology—  . 

In  all  cases  the  injury  values  and  aircraft  replacement  costs  have 
been  converted  to  1974  dollars. 

The  following  tables  (Table  B.l  and  B.2)  present  the  results  of  this 
analysis  by  jet  size.  The  word  "Significant"  in  the  titles  refers  to  the 
accident  severity  being  substantial  or  worse.  (Note:  There  are  no  accidents 
in  the  medium  jet  category  which  fit  all  the  requirements.)  A  few  cases  with 
extremely  large  personal  injury  costs  are  seen  to  dominate  the  upper  range  of 
the  total  cost  figures.  These  are  those  truly  disastrous  accidents  in  which 
a  large  number  of  deaths  takes  place. 


TABLE  B.l  SIGNIFICANT  1966-1975  SMALL  JET  ACCIDENT  COST  SUMMARY 


Cost  in 


0 

1 

1,001 

2,001 

3,001 

4,001 

5,001 

6,001 

7,001 

8,001 

9,001 

10,001 

20,001 

30,001 

40,001 


Number  of  Cases 


Thousands 


Personal 
Injury  Only 


Aircraft 
Damage  Only 


Total 

Costs 


73 

0 

-  1,000 

41 

11 

-  2,000 

4 

52 

-  3,000 

0 

43 

-  4,000 

3 

11 

-  5,000 

0 

6 

-  6,000 

3 

4 

-  7,000 

0 

7 

-  8,000 

1 

2 

-  9,000 

0 

0 

-  10,000 

0 

1 

-  20,000 

5 

2 

-  30,000 

7 

0 

-  40,000 

2 

0 

-  50,000 

0 

0 

Total 

139 

139 

0 

10 

58 

34 


1 

139 


Number  of  Cases  -  139 

Worst  Case  (Personal  Injury  Only)  - 
Worst  Case  (Aircraft  Damage  Only)  - 
Worst  Case  (Total)  -  $40,390,000 
Average  Case  (Personal  Injury  Only) 
Average  Case  (Aircraft  Damage  Only) 
Average  Case  (Total)  -  $5,096,000 


$34,140,000 

$12,000,000 

-  $2,546,000 

-  $2,550,000 
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TABLE  B.2  SIGNIFICANT  1966-1975  LARGE  JET  ACCIDENT  COST  SUMMARY 


Number  of  Cases 


Cost  in  Thousands 


Personal  Aircraft  Total 

Injury  Only  Damage  Only  Costs 


0  9 

1  -  1,000  6 

1,001  -  2,000  0 

2,001  -  3,000  0 

3,001  -  4,000  0 

4,001  -  5,000  0 

5,001  -  6,000  0 

6,001  -  7,000  0 

7,001  -  8,000  0 

8,001  -  9,000  0 

9,001  -  10,000  0 

10,001  -  20,000  0 

20,001  -  30,000  0 

30,001  -  40,000  •! 

40,001  -  50,000  __0 


0 

0 

0 

0 

0 

0 

5 

9 

0 

0 

0 

1 

1 

0 

0 


Total  16  16 


Number  of  Cases  -  16 

Worst  Case  (Personal  Injury  Only)  -  $32,502,000 
Worst  Case  (Aircraft  Damage  Only)  -  $27,000,000 
Worst  Case  (Total)  -  $49,002,000 
Average  Case  (Personal  Injury  Only)  $2,149,000 
Average  Case  (Aircraft  Damage  Only)  $8,180,000 
Average  Case  (Total)  -  $10,329,000 


1 

0 

J_ 

16 
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